R/S.5.72 (LongTerm Strage 1965) NASA (?. 2? (-:2)> 2?.2 NB. -: is half
|
|
- ゆき すみだ
- 5 years ago
- Views:
Transcription
1 SHIMURA Masato R/S References 22 C.Reiter 5 1 R/S 1.1 Harold Edwin Hurst England) Leicester, Oxford Sir Henry Lyons ( Hurst Black Simaika [Nail Basin] Jonglei 88 1
2 R/S.5.72 (LongTerm Strage 1965) NASA (?. 2? (-:2)> 2?.2 NB. -: is half _1 ( ((-:2)> 2?.2){1 _1 _1 _ _1 _ _1 1 _1 _1 _1 1 _1 1 _1 7{. <\((-:2)> 2?.2){1 _1 NB _1 _1 _1 _1 _1 1 _1 _1 1 1 _1 _ _1 _ _1 _ _ (box ( ) +/\((-:2)> 2?.2){1 _1 _1 _2 _ rw=: 3 : +/\((-:y)>y?y){1 _1.5 2
3 RW plot rwp 1 rwp (3#.5) rwp 1.5 rwp=:4 : NB. Random walk with many probability NB. Usage: plot rwp 1 NB. or.5 ;("1),. +/\ L:({(x * y)> y? (#x)#y){(l:)_1 1 ) 1.3 dt =.1, dx = ±.1 dt, dx = σ dt P(X n = dx) = µ 2σ dx 2 P(X n = dx) = 1 2 µ 2σ dx 2 X n ( 1 E(X n ) = dx 2 + µ ) ( 1 2σ dx dx µ V(X n ) = E(Xn) 2 (E(X n )) 2 = σ 2 dt µ 2 (dt) 2 X n Y n Y =, Y n = Y 1 + Y Y n ) 2σ dx 2 = dx µ dx = µdt 2σ2 3
4 script N.Thomson rho plot.1 count rno,.5, 1,1 +/\ ] a=. rn.5 1 a,. +/\ a _ _ _ _ _ _ _ _ _.7598 _ _1.886 plot((i.1)%1); +/\ rno.1 1 pd eps /temp/bwn_1.eps 4
5 Brown µ = σ = R/S long memory) R/S R t,n = max 1 k N ( k i=1 (x t+i x t,n )) min 1 k N ( k i=1 (x t+i x t,n )) S t,n R t,n S t,n Q t,n cn H MAX MIN R randsn ] a=. randsn 5 _
6 a,.((- mean)a),. +/\ (- mean) a deviation ccd _ _ _ _ _ _ _ e_17 ccd=:[: +/\ (- mean) r a r=:(>./ - <./)@ccd q a q=:r % sd NB. sd is standard deviation * ) C.Reiter S IDI 18 (n = 58) t (time span) * 1 span NB. span=: NB. n=58 span *1 wavelet 2 n 2 n 6
7 ave ave1 ave2 ave3 ave4 4. span,. span rs_test sp NB. 1/ NB. 1/ NB. 1/ span x ( (ˆ. span rs sp1) %. 1,. ˆ. span _ f = lnx.73 calc hurst x H =.769 calc_hurst sp1 _
8 1.4.2 rs_test=:([: ;._3])" 1 NB. same calc_hurst=: 3 : NB. calc_hurst randsn 1 N=: 2+i.<.-:#y RS=. N rs_test y NB. R/S static (ˆ. RS)%. 1,. ˆ. N ) n N=. 2+i.<.-:#y rs test cut 3 5 * 2 (5,:5)<;.(_3) 1? box(<) (5,:5)q ;.(_3) 1? ,: 5 ;._3]) 1? mean=. +/%#) *2 3 cut 8
9 5 ([: ;._3]) 1? rs_test=: ([: ;._3]) 5 rs_test 1? rs_test"( 1) 1? ( ) calc_hurst randsn 1 _ H sunspots C.Reiter S IDC(Solar Influence Data Analsys Center,Bergium)
10 2.2 Autocorreration ACF(k) = Σn t=k 1 (Y t Ȳ)(Y t k Ȳ) Σ n t=1 (Y t Ȳ) 2 EXAMPLE i.1 ( 9) 1 copy * 3 ({@>1+i.8),. }.}: ( (<\.a),..@(<\) a=. i.1) tmp,. +/ L: tmp=. */ L:,. <"2 > ((<\. a),.. <\ a) - L: mean a _ *3 1
11 _.75 _ _1.25 _2.25 _ _1.75 _3.75 _3.75 _ _ _2.25 _5.25 _6.25 _5.25 _2.25 _ _6.75 _8.75 _8.75 _6.75 _ _11.25 _12.25 _11.25 _ _15.75 _15.75 _ _2.25 _ (+/@:ˆ&2@(- +/%#) ) a=. i ac f acfx a _ _ _ _ _ _ Script acfx=:3 : ((+/@:*/)"2>((<\.y),..<\y) - L: mean y) % +/@:(ˆ&2)@dev y 11
12 (15-18) 15 15{.(>:i. # a),.a=.acf sp _ _ _ _ NB. runner-up NB. max _ H hwalk ( 1 2 2H h ) k+1 C.Reiter 1 2 2H 2 ( ) k h interp=: (}. + }:)@:(2: # -:) osz=: %:@-.@(2&ˆ)@+:@<:@[ sz=: osz * %@+:@<:@#@] ˆ [ 12
13 randunif=: : NB. uniform rand randsn=: * NB. normal rand randadd=: ] + sz * randsn@$@] hwalk=: 4 : x ([ randadd interp@]) ˆ: y,(x osz 1)* randsn 1 ) hwalk,h= Hwalk,H=.5 13
14 hwalk,h= R/S a=. ; }. : sp=. > jread.. plot calc_rs_sub a R/S Y R/S X N= ln2.5 ln2.5 ˆ lnr/s vs ln N.77 calc_hurst a _ ln calc_rs_masked a _
15 3.5 rs_sunspots reg lnr/s masked ln Script mean=:+/%# diff=: }. - }: NB. 2 stage Sabun ccd=: [: +/\ (- mean) NB. same avobe sd=: [: %: # % (+/@: *:@ (- mean)) r=:((>./)-<./ )@ccd NB. max - min q=: r % sd rs=: ([: q ;._3 ]) " 1 NB. R/S static rs_test=:([: ;._3])" 1 NB. same calc_hurst=: 3 : NB. calc_hurst sunspots N=: 2+i.<.-:#y RS=. N rs_test y NB. R/S static (ˆ. RS)%. 1,. ˆ. N ) NB calc_rs_sub=: 3 : X=. diff diff y N=. 2+i. <. -: # X NB. komado RS=. N rs_test X NB. rs (ˆ. N);ˆ. RS 15
16 ) calc_rs_masked=: 4 : X=. diff diff y N=. 2+i. <. -: # X NB. komado RS=. N rs_test X NB. rs MASK=. N <: ˆ x NB (ˆ. MASK # RS)%. 1,. ˆ. MASK # N ) plot_rs=: 4 : NB. find 2.5 by own eye at x axis// N RS =: calc_rs_sub y pd reset pd key rs_sunspots reg pd N; RS,: (1,. N) +/. * f=. x calc_rs_masked y pd show ) 2.5 C.Reiter C.Reiter m histories,k m-histries hen=:3 : y,1+(.3*_2{y)+_1.4**: _1{y plot hen ˆ:(4)] ref 2]\ i
17 henon dist dist %:@(+/)@:*:@:-"1 (a=.1 2 3,: 3 4 2), (-/ a),:ˆ&2 -/ a NB. A NB B _2 _2 1 NB. -/ A,B NB. ˆ& / is 9 NB. sum %: 9 is 3 NB. square 17
18 data=. hen ˆ:(4) ] m=2 k=5 size=: 35 ref=2 ]\ size{. data size ref ] x=. (35+1+i.m) {data {. ref dist x ]j=. (k=.5){. /: ref dist x x,y (m+j),. (m+j){data adress data ] coef=. ((m+j){data) %. 1,. j{ref _ f = x x 2 1,. j{ref
19 1 (>:size+i.m){data coef +/. * 1, (>:size+i.m){data 2 non-update type) param= {. a,. (a=.param,"1 ]>:i.2) fracpred data m k siz T result _ update type (c f ps1 m histry m fps1 data;(>:35+i.m){data Script fps1=:4 : ("1) NB. Usage: fps1 data m k sz =. x y x =. y NB. x=. (sz+1+i.m){y ref=. m ]\ sz {. y j=. k{./: ref dist x coef=. ((m+j){y) %. 1,.j{ref coef +/. * 1,x ) 19
20 f ps2 data y fps2=:3 : param fps1 data;y param= (2step) a=. (}.@, fps2) ˆ:(>:i.2) (>:size+i.m){data plot }.{. : a time henon,m=2 param= time henon,m=3 2
21 2.5.2 sp=. ;{: : sp NB. 58 years sunspots data (jread...) fps1 data=: sp (data ) m (m-histry) AR k m m 2 sz n-g-m G umber of steps in our goal forecast a b m c m a,b,:c Script pred_frac=: 4 : NB. Usage: pred_frac sunspots NB. x is m k sz G NB. or/ y is ; }. : sp M K SZ G =: x data=: y NB. should global definition NB calc a=. (>: SZ+M+i.G){ y b=. ((M,K,SZ),("1 )] >:i.g) fracpred y c=.(<:m)}.{. :(}.@, fps2)ˆ:(>:i.g+<:m)(>:sz+i.m){y a,b,:c ) Example pred_frac ; }. : sp plot_frac ; }. : sp m=8 k=2 sz=46 G=35 ( 21
22 a b c fractal prediction m=8,k=2 3 References C.Reiter [Fracral Visualization and J] Jsoftware 2 SIDC Solar Influence Data Analsys Center(Bergium) 22
ATM M.Shimura JCD02773@nifty.ne.jp 2003 12 13 JAPLA2003 1 queue ATM ATM queue 1.1 ATM No (Sec (Sec 1 13 37 60 26 28 99 1 25 40 39 143 202 14 88 190 27 1 184 2 170 37 40 130 317 15 121 72 28 48 115 3 101
More information1 SHIMURA Masato polynomial irr.xirr EXCEL irr
1 SHIMURA Masato 2009 12 8 1 2 1.1................................... 2 1.2 polynomial......................... 4 2 irr.xirr EXCEL 5 2.1 irr............................................. 5 2.2 d f, pv...........................................
More information1 1.1 p(x n+1 x n, x n 1, x n 2, ) = p(x n+1 x n ) (x n ) (x n+1 ) * (I Q) 1 ( 1 Q 1 Q n 0(n ) I + Q + Q 2 + = (I Q) ] q q +/. * q
Masato Shimura JCD02773@nifty.ne.jp 2008 7 23 1 2 1.1....................................... 2 1.2..................................... 2 2 3 2.1 Example...................................... 3 2.2 Script...........................................
More informationI L01( Wed) : Time-stamp: Wed 07:38 JST hig e, ( ) L01 I(2017) 1 / 19
I L01(2017-09-20 Wed) : Time-stamp: 2017-09-20 Wed 07:38 JST hig e, http://hig3.net ( ) L01 I(2017) 1 / 19 ? 1? 2? ( ) L01 I(2017) 2 / 19 ?,,.,., 1..,. 1,2,.,.,. ( ) L01 I(2017) 3 / 19 ? I. M (3 ) II,
More information0 2 SHIMURA Masato
0 2 SHIMURA Masato jcd02773@nifty.com 2009 12 8 1 1 1.1................................... 2 1.2.......................................... 3 2 2 3 2.1............................... 3 2.2.......................................
More information6.1 OOP Multi Sub a
/ WIN [ ] Masato SHIMURA JCD2773@nifty.ne.jp Last update 25 4 23 1 J 2 1.1....................................... 2 2 D 2 2.1 numeric trig................................... 6 3 6 3.1 X;Y....................................
More information80 X 1, X 2,, X n ( λ ) λ P(X = x) = f (x; λ) = λx e λ, x = 0, 1, 2, x! l(λ) = n f (x i ; λ) = i=1 i=1 n λ x i e λ i=1 x i! = λ n i=1 x i e nλ n i=1 x
80 X 1, X 2,, X n ( λ ) λ P(X = x) = f (x; λ) = λx e λ, x = 0, 1, 2, x! l(λ) = n f (x i ; λ) = n λ x i e λ x i! = λ n x i e nλ n x i! n n log l(λ) = log(λ) x i nλ log( x i!) log l(λ) λ = 1 λ n x i n =
More informationB.2 EXCEL B.3 tara B.4 CSV
Numeric Recipes for Econometrics(0) SHIMURA Masato jcd02773@nifty.com 2010 12 2 1 3-2 1.1...................................... 2 1.2 3........................... 8 1.3.......................................
More information15 P3 Pm C.Reiter dwin C.Reiter Fractal Visualization and J 4th edition fvj4 J 2D gl2 J addon Appendix (hokusai olympic0.ijs dwin * 1 coinsert *
SHIMURA Masato JCD02773@nifty.ne.jp 2017 2 23 1 2 2 6 3 9 4 15 A J 21 2 3 45 1 15 P3 Pm 1 1.1 C.Reiter dwin C.Reiter Fractal Visualization and J 4th edition fvj4 J 2D gl2 J addon Appendix (hokusai olympic0.ijs
More informationii
i 2013 5 143 5.1...................................... 143 5.2.................................. 144 5.3....................................... 148 5.4.................................. 153 5.5...................................
More information9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L) du (L) = f (9.3) dx (9.) P
9 (Finite Element Method; FEM) 9. 9. P(0) P(x) u(x) (a) P(L) f P(0) P(x) (b) 9. P(L) 9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L)
More information( 28 ) ( ) ( ) 0 This note is c 2016, 2017 by Setsuo Taniguchi. It may be used for personal or classroom purposes, but not for commercial purp
( 28) ( ) ( 28 9 22 ) 0 This ote is c 2016, 2017 by Setsuo Taiguchi. It may be used for persoal or classroom purposes, but ot for commercial purposes. i (http://www.stat.go.jp/teacher/c2epi1.htm ) = statistics
More informationii 3.,. 4. F. ( ), ,,. 8.,. 1. (75% ) (25% ) =7 24, =7 25, =7 26 (. ). 1.,, ( ). 3.,...,.,.,.,.,. ( ) (1 2 )., ( ), 0., 1., 0,.
(1 C205) 4 10 (2 C206) 4 11 (2 B202) 4 12 25(2013) http://www.math.is.tohoku.ac.jp/~obata,.,,,..,,. 1. 2. 3. 4. 5. 6. 7. 8. 1., 2007 ( ).,. 2. P. G., 1995. 3. J. C., 1988. 1... 2.,,. ii 3.,. 4. F. ( ),..
More information1 1 ( ) ( 1.1 1.1.1 60% mm 100 100 60 60% 1.1.2 A B A B A 1
1 21 10 5 1 E-mail: qliu@res.otaru-uc.ac.jp 1 1 ( ) ( 1.1 1.1.1 60% mm 100 100 60 60% 1.1.2 A B A B A 1 B 1.1.3 boy W ID 1 2 3 DI DII DIII OL OL 1.1.4 2 1.1.5 1.1.6 1.1.7 1.1.8 1.2 1.2.1 1. 2. 3 1.2.2
More information1 bmp gif,png,jpg bmp gif,png jpg BPG 2014 jpg *3 RAW TIFF RAW CCD CMOS R,G,B TIFF net *4 1.1 JPEG HP JPEG 3 1 4, 1 8, 1 16 JPEG SD jpeg JPEG RGB YCrC
Viewmat SHIMURA Masato 2015 6 12 viewnmat viewmat QT J8x 400 RGB CMYK *1 *2 RGB CMYK *1 CMYK,, *2 1 1 bmp gif,png,jpg bmp gif,png jpg BPG 2014 jpg *3 RAW TIFF RAW CCD CMOS R,G,B TIFF net *4 1.1 JPEG HP
More information分布
(normal distribution) 30 2 Skewed graph 1 2 (variance) s 2 = 1/(n-1) (xi x) 2 x = mean, s = variance (variance) (standard deviation) SD = SQR (var) or 8 8 0.3 0.2 0.1 0.0 0 1 2 3 4 5 6 7 8 8 0 1 8 (probability
More informationii 3.,. 4. F. (), ,,. 8.,. 1. (75%) (25%) =7 20, =7 21 (. ). 1.,, (). 3.,. 1. ().,.,.,.,.,. () (12 )., (), 0. 2., 1., 0,.
24(2012) (1 C106) 4 11 (2 C206) 4 12 http://www.math.is.tohoku.ac.jp/~obata,.,,,.. 1. 2. 3. 4. 5. 6. 7.,,. 1., 2007 (). 2. P. G. Hoel, 1995. 3... 1... 2.,,. ii 3.,. 4. F. (),.. 5... 6.. 7.,,. 8.,. 1. (75%)
More informationJ6 M.Shimura (1) 1 2 (2) (1824) ( (1842) 1 (1) 1.1 C.Reiter dwin require ad
J6 M.Shimura JCD02773@nifty.ne.jp 2011 12 13 1 (1) 1 2 (2)- 12 3 14 4-20 5 24 6 33 60 (1824) ( 100 1760 1849 (1842) 1 (1) 1.1 C.Reiter dwin require addons/graphics/fvj3/dwin2.ijs 1 xy find_maxmin 4 5 calc_each_poly
More informationこんにちは由美子です
1 2 . sum Variable Obs Mean Std. Dev. Min Max ---------+----------------------------------------------------- var1 13.4923077.3545926.05 1.1 3 3 3 0.71 3 x 3 C 3 = 0.3579 2 1 0.71 2 x 0.29 x 3 C 2 = 0.4386
More information平成13年度日本分析センター年報
200 150 70 234 Bq m 3 1 148 Bq m -3 100 0 550 0 11/1 0:00 am 11/2 0:00 am 11/3 0:00 am 25 20 15 10 11/1 0:00 am 11/2 0:00 am 11/3 0:00 am 39.2 Bq m -3 11/4 0:00 am 30 990 19.3 Bq m -3 60 15.8 Bq m -3 14.1
More information1 (1) vs. (2) (2) (a)(c) (a) (b) (c) 31 2 (a) (b) (c) LENCHAR
() 601 1 () 265 OK 36.11.16 20 604 266 601 30.4.5 (1) 91621 3037 (2) 20-12.2 20-13 (3) ex. 2540-64 - LENCHAR 1 (1) vs. (2) (2) 605 50.2.13 41.4.27 10 10 40.3.17 (a)(c) 2 1 10 (a) (b) (c) 31 2 (a) (b) (c)
More informationuntitled
yoshi@image.med.osaka-u.ac.jp http://www.image.med.osaka-u.ac.jp/member/yoshi/ II Excel, Mathematica Mathematica Osaka Electro-Communication University (2007 Apr) 09849-31503-64015-30704-18799-390 http://www.image.med.osaka-u.ac.jp/member/yoshi/
More information5 Armitage x 1,, x n y i = 10x i + 3 y i = log x i {x i } {y i } 1.2 n i i x ij i j y ij, z ij i j 2 1 y = a x + b ( cm) x ij (i j )
5 Armitage. x,, x n y i = 0x i + 3 y i = log x i x i y i.2 n i i x ij i j y ij, z ij i j 2 y = a x + b 2 2. ( cm) x ij (i j ) (i) x, x 2 σ 2 x,, σ 2 x,2 σ x,, σ x,2 t t x * (ii) (i) m y ij = x ij /00 y
More informationa b GE(General Erectrics) 9 4 irr (JAPLA 2009/12) Example1 120 P = C r + C 2 (1 + r) C t 1 (1 + r) t 1 + C t + F (1 + r) t 10
1 SHIMURA Masato 2010 9 27 1 1 2 CF 6 3 10 *1 irr irr irr(inner rate of return)function is able to written only few lines,and it is very powerful and useful for simulate unprofitable business model. 1
More information応力とひずみ.ppt
in yukawa@numse.nagoya-u.ac.jp 2 3 4 5 x 2 6 Continuum) 7 8 9 F F 10 F L L F L 1 L F L F L F 11 F L F F L F L L L 1 L 2 12 F L F! A A! S! = F S 13 F L L F F n = F " cos# F t = F " sin# S $ = S cos# S S
More information5 n P j j (P i,, P k, j 1) 1 n n ) φ(n) = n (1 1Pj [ ] φ φ P j j P j j = = = = = n = φ(p j j ) (P j j P j 1 j ) P j j ( 1 1 P j ) P j j ) (1 1Pj (1 1P
p P 1 n n n 1 φ(n) φ φ(1) = 1 1 n φ(n), n φ(n) = φ()φ(n) [ ] n 1 n 1 1 n 1 φ(n) φ() φ(n) 1 3 4 5 6 7 8 9 1 3 4 5 6 7 8 9 1 4 5 7 8 1 4 5 7 8 10 11 1 13 14 15 16 17 18 19 0 1 3 4 5 6 7 19 0 1 3 4 5 6 7
More informationxyr x y r x y r u u
xyr x y r x y r u u y a b u a b a b c d e f g u a b c d e g u u e e f yx a b a b a b c a b c a b a b c a b a b c a b c a b c a u xy a b u a b c d a b c d u ar ar a xy u a b c a b c a b p a b a b c a
More informationH22 BioS (i) I treat1 II treat2 data d1; input group patno treat1 treat2; cards; ; run; I
H BioS (i) I treat II treat data d; input group patno treat treat; cards; 8 7 4 8 8 5 5 6 ; run; I II sum data d; set d; sum treat + treat; run; sum proc gplot data d; plot sum * group ; symbol c black
More information最小2乗法
2 2012 4 ( ) 2 2012 4 1 / 42 X Y Y = f (X ; Z) linear regression model X Y slope X 1 Y (X, Y ) 1 (X, Y ) ( ) 2 2012 4 2 / 42 1 β = β = β (4.2) = β 0 + β (4.3) ( ) 2 2012 4 3 / 42 = β 0 + β + (4.4) ( )
More information10
z c j = N 1 N t= j1 [ ( z t z ) ( )] z t j z q 2 1 2 r j /N j=1 1/ N J Q = N(N 2) 1 N j j=1 r j 2 2 χ J B d z t = z t d (1 B) 2 z t = (z t z t 1 ) (z t 1 z t 2 ) (1 B s )z t = z t z t s _ARIMA CONSUME
More informationtokei01.dvi
2. :,,,. :.... Apr. - Jul., 26FY Dept. of Mechanical Engineering, Saga Univ., JAPAN 4 3. (probability),, 1. : : n, α A, A a/n. :, p, p Apr. - Jul., 26FY Dept. of Mechanical Engineering, Saga Univ., JAPAN
More information68 A mm 1/10 A. (a) (b) A.: (a) A.3 A.4 1 1
67 A Section A.1 0 1 0 1 Balmer 7 9 1 0.1 0.01 1 9 3 10:09 6 A.1: A.1 1 10 9 68 A 10 9 10 9 1 10 9 10 1 mm 1/10 A. (a) (b) A.: (a) A.3 A.4 1 1 A.1. 69 5 1 10 15 3 40 0 0 ¾ ¾ É f Á ½ j 30 A.3: A.4: 1/10
More information.3. (x, x = (, u = = 4 (, x x = 4 x, x 0 x = 0 x = 4 x.4. ( z + z = 8 z, z 0 (z, z = (0, 8, (,, (8, 0 3 (0, 8, (,, (8, 0 z = z 4 z (g f(x = g(
06 5.. ( y = x x y 5 y 5 = (x y = x + ( y = x + y = x y.. ( Y = C + I = 50 + 0.5Y + 50 r r = 00 0.5Y ( L = M Y r = 00 r = 0.5Y 50 (3 00 0.5Y = 0.5Y 50 Y = 50, r = 5 .3. (x, x = (, u = = 4 (, x x = 4 x,
More informationMicrosoft Word - 倫理 第40,43,45,46講 テキスト.docx
6 538 ( 552 ) (1) () (2) () ( )( ) 1 vs () (1) (2) () () () ) ()() (3) () ( () 2 () () () ()( ) () (7) (8) () 3 4 5 abc b c 6 a (a) b b ()() 7 c (c) ()() 8 9 10 () 1 ()()() 2 () 3 1 1052 1051 () 1053 11
More informationk2 ( :35 ) ( k2) (GLM) web web 1 :
2012 11 01 k2 (2012-10-26 16:35 ) 1 6 2 (2012 11 01 k2) (GLM) kubo@ees.hokudai.ac.jp web http://goo.gl/wijx2 web http://goo.gl/ufq2 1 : 2 2 4 3 7 4 9 5 : 11 5.1................... 13 6 14 6.1......................
More information分散分析・2次元正規分布
2 II L10(2016-06-30 Thu) : Time-stamp: 2016-06-30 Thu 13:55 JST hig F 2.. http://hig3.net ( ) L10 2 II(2016) 1 / 24 F 2 F L09-Q1 Quiz :F 1 α = 0.05, 2 F 3 H 0, : σ 2 1 /σ2 2 = 1., H 1, σ 2 1 /σ2 2 1. 4
More information9 2 1 f(x, y) = xy sin x cos y x y cos y y x sin x d (x, y) = y cos y (x sin x) = y cos y(sin x + x cos x) x dx d (x, y) = x sin x (y cos y) = x sin x
2009 9 6 16 7 1 7.1 1 1 1 9 2 1 f(x, y) = xy sin x cos y x y cos y y x sin x d (x, y) = y cos y (x sin x) = y cos y(sin x + x cos x) x dx d (x, y) = x sin x (y cos y) = x sin x(cos y y sin y) y dy 1 sin
More information1 1 1 1-1 1 1-9 1-3 1-1 13-17 -3 6-4 6 3 3-1 35 3-37 3-3 38 4 4-1 39 4- Fe C TEM 41 4-3 C TEM 44 4-4 Fe TEM 46 4-5 5 4-6 5 5 51 6 5 1 1-1 1991 1,1 multiwall nanotube 1993 singlewall nanotube ( 1,) sp 7.4eV
More informationこんにちは由美子です
Analysis of Variance 2 two sample t test analysis of variance (ANOVA) CO 3 3 1 EFV1 µ 1 µ 2 µ 3 H 0 H 0 : µ 1 = µ 2 = µ 3 H A : Group 1 Group 2.. Group k population mean µ 1 µ µ κ SD σ 1 σ σ κ sample mean
More informationkubostat2017e p.1 I 2017 (e) GLM logistic regression : : :02 1 N y count data or
kubostat207e p. I 207 (e) GLM kubo@ees.hokudai.ac.jp https://goo.gl/z9ycjy 207 4 207 6:02 N y 2 binomial distribution logit link function 3 4! offset kubostat207e (https://goo.gl/z9ycjy) 207 (e) 207 4
More information00 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.... 0........ 0 0 0 0 0 0 0 0 0 0..0..........0 0 0 0 0 0 0 0 0 0 0.... 0........ 0 0 0 0 0 0 0 0 0 0... 0...... 0... 0 0 0 0 0 0..0 0... 0 0 0 0 0.0.....0.
More information浜松医科大学紀要
On the Statistical Bias Found in the Horse Racing Data (1) Akio NODA Mathematics Abstract: The purpose of the present paper is to report what type of statistical bias the author has found in the horse
More information28
y i = Z i δ i +ε i ε i δ X y i = X Z i δ i + X ε i [ ] 1 δ ˆ i = Z i X( X X) 1 X Z i [ ] 1 σ ˆ 2 Z i X( X X) 1 X Z i Z i X( X X) 1 X y i σ ˆ 2 ˆ σ 2 = [ ] y i Z ˆ [ i δ i ] 1 y N p i Z i δ ˆ i i RSTAT
More information統計学のポイント整理
.. September 17, 2012 1 / 55 n! = n (n 1) (n 2) 1 0! = 1 10! = 10 9 8 1 = 3628800 n k np k np k = n! (n k)! (1) 5 3 5 P 3 = 5! = 5 4 3 = 60 (5 3)! n k n C k nc k = npk k! = n! k!(n k)! (2) 5 3 5C 3 = 5!
More information自由集会時系列part2web.key
spurious correlation spurious regression xt=xt-1+n(0,σ^2) yt=yt-1+n(0,σ^2) n=20 type1error(5%)=0.4703 no trend 0 1000 2000 3000 4000 p for r xt=xt-1+n(0,σ^2) random walk random walk variable -5 0 5 variable
More informationuntitled
23 59 13 23 24 0101 0001 0101 0002 0101 0001 0101 0002 0101 0007 0101 0009 0101 0012 0101 0026 0101 0031 0101 0033 0101 0056 0101 0059 0101 0075 0101 0076 0101 5001 0101 0002 0101 0003 0101 0008 0101 0010
More informationsec13.dvi
13 13.1 O r F R = m d 2 r dt 2 m r m = F = m r M M d2 R dt 2 = m d 2 r dt 2 = F = F (13.1) F O L = r p = m r ṙ dl dt = m ṙ ṙ + m r r = r (m r ) = r F N. (13.2) N N = R F 13.2 O ˆn ω L O r u u = ω r 1 1:
More informationuntitled
2011/6/22 M2 1*1+2*2 79 2F Y YY 0.0 0.2 0.4 0.6 0.8 0.000 0.002 0.004 0.006 0.008 0.010 0.012 1.0 1.5 2.0 2.5 3.0 3.5 4.0 Y 0 50 100 150 200 250 YY A (Y = X + e A ) B (YY = X + e B ) X 0.00 0.05 0.10
More informationuntitled
Ÿ Ÿ ( œ ) 120,000 60,000 120,000 120,000 80,000 72,000 100,000 180,000 60,000 100,000 60,000 120,000 100,000 240,000 120,000 240,000 1,150,000 100,000 120,000 72,000 300,000 72,000 100,000 100,000 60,000
More information統計的データ解析
ds45 xspec qdp guplot oocalc (Error) gg (Radom Error)(Systematc Error) x, x,, x ( x, x,..., x x = s x x µ = lm = σ µ x x = lm ( x ) = σ ( ) = - x = js j ( ) = j= ( j) x x + xj x + xj j x + xj = ( x x
More informationN88 BASIC 0.3 C: My Documents 0.6: 0.3: (R) (G) : enterreturn : (F) BA- SIC.bas 0.8: (V) 0.9: 0.5:
BASIC 20 4 10 0 N88 Basic 1 0.0 N88 Basic..................................... 1 0.1............................................... 3 1 4 2 5 3 6 4 7 5 10 6 13 7 14 0 N88 Basic 0.0 N88 Basic 0.1: N88Basic
More informationexample2_time.eps
Google (20/08/2 ) ( ) Random Walk & Google Page Rank Agora on Aug. 20 / 67 Introduction ( ) Random Walk & Google Page Rank Agora on Aug. 20 2 / 67 Introduction Google ( ) Random Walk & Google Page Rank
More informationuntitled
2 : n =1, 2,, 10000 0.5125 0.51 0.5075 0.505 0.5025 0.5 0.4975 0.495 0 2000 4000 6000 8000 10000 2 weak law of large numbers 1. X 1,X 2,,X n 2. µ = E(X i ),i=1, 2,,n 3. σi 2 = V (X i ) σ 2,i=1, 2,,n ɛ>0
More information2 A I / 58
2 A 2018.07.12 I 2 2018.07.12 1 / 58 I 2 2018.07.12 2 / 58 π-computer gnuplot 5/31 1 π-computer -X ssh π-computer gnuplot I 2 2018.07.12 3 / 58 gnuplot> gnuplot> plot sin(x) I 2 2018.07.12 4 / 58 cp -r
More informationpower.tex
Contents ii 1... 1... 1... 7... 7 3 (DFFT).................................... 8 4 (CIFT) DFFT................................ 10 5... 13 6... 16 3... 0 4... 0 5... 0 6... 0 i 1987 SN1987A 0.5 X SN1987A
More informationx = a 1 f (a r, a + r) f(a) r a f f(a) 2 2. (a, b) 2 f (a, b) r f(a, b) r (a, b) f f(a, b)
2011 I 2 II III 17, 18, 19 7 7 1 2 2 2 1 2 1 1 1.1.............................. 2 1.2 : 1.................... 4 1.2.1 2............................... 5 1.3 : 2.................... 5 1.3.1 2.....................................
More information: , 2.0, 3.0, 2.0, (%) ( 2.
2017 1 2 1.1...................................... 2 1.2......................................... 4 1.3........................................... 10 1.4................................. 14 1.5..........................................
More information(lm) lm AIC 2 / 1
W707 s-taiji@is.titech.ac.jp 1 / 1 (lm) lm AIC 2 / 1 : y = β 1 x 1 + β 2 x 2 + + β d x d + β d+1 + ϵ (ϵ N(0, σ 2 )) y R: x R d : β i (i = 1,..., d):, β d+1 : ( ) (d = 1) y = β 1 x 1 + β 2 + ϵ (d > 1) y
More information読めば必ずわかる 分散分析の基礎 第2版
2 2003 12 5 ( ) ( ) 2 I 3 1 3 2 2? 6 3 11 4? 12 II 14 5 15 6 16 7 17 8 19 9 21 10 22 11 F 25 12 : 1 26 3 I 1 17 11 x 1, x 2,, x n x( ) x = 1 n n i=1 x i 12 (SD ) x 1, x 2,, x n s 2 s 2 = 1 n n (x i x)
More information¥¤¥ó¥¿¡¼¥Í¥Ã¥È·×¬¤È¥Ç¡¼¥¿²òÀÏ Âè2²ó
2 2015 4 20 1 (4/13) : ruby 2 / 49 2 ( ) : gnuplot 3 / 49 1 1 2014 6 IIJ / 4 / 49 1 ( ) / 5 / 49 ( ) 6 / 49 (summary statistics) : (mean) (median) (mode) : (range) (variance) (standard deviation) 7 / 49
More informationJAB RL504:201 JAB NOTE 4 2 January 25, A B
JAB RL504:201 JAB NOTE 4 JAB RL504:201 2201125 1200025 200-0-25-1/- 2 201-01-25 JAB RL504:201 JAB NOTE 4 2 January 25, 201............... 6 5.1... 6 5.2... 6 5. A... 7 5.4 B... 7 5.5... 8 5.6... 9 5.7...
More informationuntitled
SPring-8 RFgun JASRI/SPring-8 6..7 Contents.. 3.. 5. 6. 7. 8. . 3 cavity γ E A = er 3 πε γ vb r B = v E c r c A B A ( ) F = e E + v B A A A A B dp e( v B+ E) = = m d dt dt ( γ v) dv e ( ) dt v B E v E
More information卒業論文
Y = ax 1 b1 X 2 b2...x k bk e u InY = Ina + b 1 InX 1 + b 2 InX 2 +...+ b k InX k + u X 1 Y b = ab 1 X 1 1 b 1 X 2 2...X bk k e u = b 1 (ax b1 1 X b2 2...X bk k e u ) / X 1 = b 1 Y / X 1 X 1 X 1 q YX1
More informationkubostat2017b p.1 agenda I 2017 (b) probability distribution and maximum likelihood estimation :
kubostat2017b p.1 agenda I 2017 (b) probabilit distribution and maimum likelihood estimation kubo@ees.hokudai.ac.jp http://goo.gl/76c4i 2017 11 14 : 2017 11 07 15:43 1 : 2 3? 4 kubostat2017b (http://goo.gl/76c4i)
More informationi
009 I 1 8 5 i 0 1 0.1..................................... 1 0.................................................. 1 0.3................................. 0.4........................................... 3
More information2 I I / 61
2 I 2017.07.13 I 2 2017.07.13 1 / 61 I 2 2017.07.13 2 / 61 I 2 2017.07.13 3 / 61 7/13 2 7/20 I 7/27 II I 2 2017.07.13 4 / 61 π-computer gnuplot MobaXterm Wiki PC X11 DISPLAY I 2 2017.07.13 5 / 61 Mac 1.
More information( ) 1.1 Polychoric Correlation Polyserial Correlation Graded Response Model Partial Credit Model Tetrachoric Correlation ( ) 2 x y x y s r 1 x 2
1 (,2007) SPSSver8 1997 (2002) 1. 2. polychoric correlation coefficient (polyserial correlation coefficient) 3. (1999) M-plus R 1 ( ) 1.1 Polychoric Correlation Polyserial Correlation Graded Response Model
More informationII A A441 : October 02, 2014 Version : Kawahira, Tomoki TA (Kondo, Hirotaka )
II 214-1 : October 2, 214 Version : 1.1 Kawahira, Tomoki TA (Kondo, Hirotaka ) http://www.math.nagoya-u.ac.jp/~kawahira/courses/14w-biseki.html pdf 1 2 1 9 1 16 1 23 1 3 11 6 11 13 11 2 11 27 12 4 12 11
More informationEPSON LP-8900ユーザーズガイド
3 4 5 6 7 8 abc ade w p s 9 10 s s 11 p 12 p 13 14 p s 15 p s A B 16 w 17 C p 18 D E F 19 p w G H 20 A B 21 C s p D 22 E s p w 23 w w s 24 p w s 25 w 26 p p 27 w p s 28 w p 29 w p s 30 p s 31 A s B 32
More information数値計算:常微分方程式
( ) 1 / 82 1 2 3 4 5 6 ( ) 2 / 82 ( ) 3 / 82 C θ l y m O x mg λ ( ) 4 / 82 θ t C J = ml 2 C mgl sin θ θ C J θ = mgl sin θ = θ ( ) 5 / 82 ω = θ J ω = mgl sin θ ω J = ml 2 θ = ω, ω = g l sin θ = θ ω ( )
More information1 J 2 tasu =: + (Tacit definition) (Explicit definition) 1.1 (&) x u&v y Fork Bond & Bond(&) 0&{ u u v v v y x y 1&{ ( p) ( q) x v&
1 J SHIMURA Masato jcd02773@nifty.ne.jp 2008 12 8 1 J 1 2 J 4 3 5 4 8 5 /de Morgan law 11 6 16 7 19 8 Reference 21 A 21 J 5 1 J J Atom ) APL J 1 J 2 tasu =: + (Tacit definition) (Explicit definition) 1.1
More informationランダムウォークの境界条件・偏微分方程式の数値計算
B L06(2018-05-22 Tue) : Time-stamp: 2018-05-22 Tue 21:53 JST hig,, 2, multiply transf http://hig3.net L06 B(2018) 1 / 38 L05-Q1 Quiz : 1 M λ 1 = 1 u 1 ( ). M u 1 = u 1, u 1 = ( 3 4 ) s (s 0)., u 1 = 1
More informationII 2 II
II 2 II 2005 yugami@cc.utsunomiya-u.ac.jp 2005 4 1 1 2 5 2.1.................................... 5 2.2................................. 6 2.3............................. 6 2.4.................................
More informationSample function Re random process Flutter, Galloping, etc. ensemble (mean value) N 1 µ = lim xk( t1) N k = 1 N autocorrelation function N 1 R( t1, t1
Sample function Re random process Flutter, Galloping, etc. ensemble (mean value) µ = lim xk( k = autocorrelation function R( t, t + τ) = lim ( ) ( + τ) xk t xk t k = V p o o R p o, o V S M R realization
More information,. Black-Scholes u t t, x c u 0 t, x x u t t, x c u t, x x u t t, x + σ x u t, x + rx ut, x rux, t 0 x x,,.,. Step 3, 7,,, Step 6., Step 4,. Step 5,,.
9 α ν β Ξ ξ Γ γ o δ Π π ε ρ ζ Σ σ η τ Θ θ Υ υ ι Φ φ κ χ Λ λ Ψ ψ µ Ω ω Def, Prop, Th, Lem, Note, Remark, Ex,, Proof, R, N, Q, C [a, b {x R : a x b} : a, b {x R : a < x < b} : [a, b {x R : a x < b} : a,
More informationCMP Technical Report No. 4 Department of Computational Nanomaterials Design ISIR, Osaka University 2 2................................. 2.2......................... 2 3 3 3................................
More information¥¤¥ó¥¿¡¼¥Í¥Ã¥È·×¬¤È¥Ç¡¼¥¿²òÀÏ Âè2²ó
2 212 4 13 1 (4/6) : ruby 2 / 35 ( ) : gnuplot 3 / 35 ( ) 4 / 35 (summary statistics) : (mean) (median) (mode) : (range) (variance) (standard deviation) 5 / 35 (mean): x = 1 n (median): { xr+1 m, m = 2r
More informationstatstcs statstcum (EBM) 2 () : ( )GDP () : () : POS STEP 1: STEP 2: STEP 3: STEP 4: 3 2
(descrptve statstcs) 2010 9 3 1 1 2 2 3 2 3.1............................................. 3 3.2............................................. 3 4 4 4.1........................................ 5 5 6 6 -pvot
More informationuntitled
š ( ) 300,000 180,000 100,000 120,000 60,000 120,000 240,000 120,000 170,000 240,000 100,000 99,000 120,000 72,000 100,000 450,000 72,000 60,000 100,000 100,000 60,000 60,000 100,000 200,000 60,000 124,000
More information時系列解析と自己回帰モデル
B L11(2017-07-03 Mon) : Time-stamp: 2017-07-03 Mon 11:04 JST hig,,,.,. http://hig3.net ( ) L11 B(2017) 1 / 28 L10-Q1 Quiz : 1 6 6., x[]={1,1,3,3,3,8}; (. ) 2 x = 0, 1, 2,..., 9 10, 10. u[]={0,2,0,3,0,0,0,0,1,0};
More informationi 18 2H 2 + O 2 2H 2 + ( ) 3K
i 18 2H 2 + O 2 2H 2 + ( ) 3K ii 1 1 1.1.................................. 1 1.2........................................ 3 1.3......................................... 3 1.4....................................
More informationOn model selection problems in terms of prediction mean squared error and interpretaion of AIC (slides)
Applications in Econometrics and Finance by Long Memory Processes 2007 11 6 13:30-16:10 Table of Contents PART1 PART2 PART3 PART1 1 {y t } ρ(h) = ( h ) = Cov[y t y t+h ]/ Var[y t ] (yt y)(y t+h y) ρ(h)
More informationU(a, b) U(a, b) ( a, b ) X U(0, 1) Y U(a, b) Y = (b a)x + a S 1 = x 2 dx y = 1 x 2 y x N(0, 1) randn 0 1 p(x) = 1 e x2 2 2π x [ 0, 1 ],
uniform.m 1 2 3 4 >> uniform 0.81472 0.90579 0.12699 0.91338 0.63236 0.09754 0.2785 0.54688 0.95751 0.96489 >> uniform 0.15761 0.97059 0.95717 0.48538 dice.m 0.80028 0.14189 0.42176 0.91574 0.79221 function
More informationAR(1) y t = φy t 1 + ɛ t, ɛ t N(0, σ 2 ) 1. Mean of y t given y t 1, y t 2, E(y t y t 1, y t 2, ) = φy t 1 2. Variance of y t given y t 1, y t
87 6.1 AR(1) y t = φy t 1 + ɛ t, ɛ t N(0, σ 2 ) 1. Mean of y t given y t 1, y t 2, E(y t y t 1, y t 2, ) = φy t 1 2. Variance of y t given y t 1, y t 2, V(y t y t 1, y t 2, ) = σ 2 3. Thus, y t y t 1,
More informationfx-3650P_fx-3950P_J
SA1109-E J fx-3650p fx-3950p http://edu.casio.jp RCA500002-001V04 AB2 Mode
More informationrenshumondai-kaito.dvi
3 1 13 14 1.1 1 44.5 39.5 49.5 2 0.10 2 0.10 54.5 49.5 59.5 5 0.25 7 0.35 64.5 59.5 69.5 8 0.40 15 0.75 74.5 69.5 79.5 3 0.15 18 0.90 84.5 79.5 89.5 2 0.10 20 1.00 20 1.00 2 1.2 1 16.5 20.5 12.5 2 0.10
More information24 201170068 1 4 2 6 2.1....................... 6 2.1.1................... 6 2.1.2................... 7 2.1.3................... 8 2.2..................... 8 2.3................. 9 2.3.1........... 12
More informationsolutionJIS.dvi
May 0, 006 6 morimune@econ.kyoto-u.ac.jp /9/005 (7 0/5/006 1 1.1 (a) (b) (c) c + c + + c = nc (x 1 x)+(x x)+ +(x n x) =(x 1 + x + + x n ) nx = nx nx =0 c(x 1 x)+c(x x)+ + c(x n x) =c (x i x) =0 y i (x
More information1 (1) () (3) I 0 3 I I d θ = L () dt θ L L θ I d θ = L = κθ (3) dt κ T I T = π κ (4) T I κ κ κ L l a θ L r δr δl L θ ϕ ϕ = rθ (5) l
1 1 ϕ ϕ ϕ S F F = ϕ (1) S 1: F 1 1 (1) () (3) I 0 3 I I d θ = L () dt θ L L θ I d θ = L = κθ (3) dt κ T I T = π κ (4) T I κ κ κ L l a θ L r δr δl L θ ϕ ϕ = rθ (5) l : l r δr θ πrδr δf (1) (5) δf = ϕ πrδr
More informationkubostat2015e p.2 how to specify Poisson regression model, a GLM GLM how to specify model, a GLM GLM logistic probability distribution Poisson distrib
kubostat2015e p.1 I 2015 (e) GLM kubo@ees.hokudai.ac.jp http://goo.gl/76c4i 2015 07 22 2015 07 21 16:26 kubostat2015e (http://goo.gl/76c4i) 2015 (e) 2015 07 22 1 / 42 1 N k 2 binomial distribution logit
More informationα β *2 α α β β α = α 1 β = 1 β 2.2 α 0 β *3 2.3 * *2 *3 *4 (µ A ) (µ P ) (µ A > µ P ) 10 (µ A = µ P + 10) 15 (µ A = µ P +
Armitage 1 1.1 2 t *1 α β 1.2 µ x µ 2 2 2 α β 2.1 1 α β α ( ) β *1 t t 1 α β *2 α α β β α = α 1 β = 1 β 2.2 α 0 β 1 0 0 1 1 5 2.5 *3 2.3 *4 3 3.1 1 1 1 *2 *3 *4 (µ A ) (µ P ) (µ A > µ P ) 10 (µ A = µ P
More informationgenron-3
" ( K p( pasals! ( kg / m 3 " ( K! v M V! M / V v V / M! 3 ( kg / m v ( v "! v p v # v v pd v ( J / kg p ( $ 3! % S $ ( pv" 3 ( ( 5 pv" pv R" p R!" R " ( K ( 6 ( 7 " pv pv % p % w ' p% S & $ p% v ( J /
More information26.fx95MS_Etype_J-cover_SA0311D
P fx-95ms fx-100ms fx-570ms fx-912ms (fx-115ms) fx-991ms English Manual Readers! Please be sure to read the important notice on the inside of the front cover of this manual. J http://www.casio.co.jp/edu/
More information/02/18
3 09/0/8 i III,,,, III,?,,,,,,,,,,,,,,,,,,,,?,?,,,,,,,,,,,,,,!!!,? 3,,,, ii,,,!,,,, OK! :!,,,, :!,,,,,, 3:!,, 4:!,,,, 5:!,,! 7:!,,,,, 8:!,! 9:!,,,,,,,,, ( ),, :, ( ), ( ), 6:!,,, :... : 3 ( )... iii,,
More information