DAA03
|
|
- ちえこ うるしはた
- 5 years ago
- Views:
Transcription
1
2 par(mfrow=c(1,2)) # figure Dist. of Height for Female Participants Dist. of Height for Male Participants Density Density Height Height
3 HISTOGRAM hist(dat[dat$gender=='f',]$h, main="dist. of Height for Female Participants", xlab="height", xlim=c(140,190), probability=t) # dat[dat$gender=='f',]$h # dat gender F dens.f = density(dat[dat$gender=='f',]$h) lines(dens.f, col='blue',lwd=2) hist(dat[dat$gender== M,]$h, main= Dist. of Height for Male Participants, xlab= Height, xlim=c(140,190), probability=t,ylim=c(0,0.08)) dens.m = density(dat[dat$gender=='m',]$h) lines(dens.m, col='green', lwd=2)
4 HISTOGRAM Dist. of Height for Female Participants Dist. of Height for Male Participants Density Density Height Height
5 par(mfrow=c(1,1)) plot(dens.f,col='blue',lwd=2, ylab='density', xlim=c(140,190), main="dist. of Height by gender",xlab='height') lines(dens.m,col='green',lwd=2) legend("topleft", c('female','male'), col=c('blue','green'), cex=1.5,lwd=2)
6 # text(x, y, TEXT ) text(x = 157.5, y = 0.04, 'Female', col='blue', cex=2) text(x = 170, y = 0.04,'Male', col='green', cex=2) Dist. of Height by gender density Female Male Height
7
8 # plot(x,y, ) plot(dat$shoesize, dat$h, main="relationship b/w shoesize and height, xlab = 'shoesize, ylab='height, pch=19, col='red ) Relationship b/w shoesize and height height shoesize
9 txt = paste("r =",round(cor(dat$shoesize,dat$h), 4)) > txt[1] "r = 0.875" # paste a=1:3 paste("s",a) => [1] "s 1" "s 2" "s 3" paste( s, a, sep= ) => [1] "s1" "s2" "s3 # round > cor(dat$shoesize,dat$h) [1] > round(cor(dat$shoesize,dat$h),4) [1] 0.875
10 txt = paste("r =", round(cor(dat$shoesize,dat$h), 4)) text(22, 175, txt, cex = 1.5) Relationship b/w shoesize and height height r = shoesize
11 abline: abline(h = mean(dat$h), col='blue'); abline(v = mean(dat$shoesize), col='green'); H V height r = mean height Relationship b/w shoesize and height mean shoesize shoesize
12 abline(lm(dat$h~dat$shoesize), lty=2, lwd=2) Relationship b/w shoesize and height height r = shoesize
13 plot(dat[dat$gender=='f',]$shoesize, dat[dat$gender=='f',]$h, main="relationship b/w shoesize and height", xlab='shoesize', ylab='height', cex.lab=1.5, pch=19, col='blue', xlim=c(20,29), ylim=c(140,190)) Relationship b/w shoesize and height height Female Male shoesize
14 plot(dat[dat$gender=='m',]$shoesize, dat[dat$gender=='m',]$h, main="relationship b/w shoesize and height", xlab='shoesize', ylab='height', cex.lab=1.5, pch=15, col='green', xlim=c(20,29), ylim=c(140,190)) Relationship b/w shoesize and height height Female Male shoesize
15 Relationship b/w shoesize and height height Female Male shoesize
16 plot(dat[dat$gender=='f',]$shoesize, dat[dat$gender=='f',]$h, main="relationship b/w shoesize and height", xlab='shoesize', ylab='height', cex.lab=1.5, pch=19, col='blue', xlim=c(20,29), ylim=c(140,190)) points(dat[dat$gender=='m',]$shoesize,dat[dat$gender=='m',]$h, pch = 15, col = 'green') legend("topleft", c('female','male'), pch =c(19,15), col = c('blue','green'), cex = 1.5) points: lines height Female Male Relationship b/w shoesize and height shoesize
17 dat<-read.csv(" plot(dat, pch=20, col='blue') material price design sales
18 plot(dat.pca, pch = rownames(dat.pca), cex = 1.7, col = 'blue ) > dat.pca writing thesis interview a b c d e f g h i d d writing c b g h i e b c g h i e a a f f e d d e i i a h g b thesis b g ha c c f f e d e d i i a h a h interview c g b g b f f c
19
20
21
22 ( "% " = $ %&' )
23 ! = # $ = $ % & % + $ ( & ( + + $ * & * = + $, &, *,-%! = # $ = $ % & % + $ ( & ( + + $ * & * =. 0 /0 $1 $ 2$
24
25 ! " # = %&' ( = ) ( + # = ) ( ) ( # %&' ( = (, + # -, + ( # + # - # + + ( 0 + # - 0 = 1 = 1 2 = 1 2 ( 2 # 2( # - 2 = 1 2 ( 2 # = ) ( # # 5 1 = ) ( # ) ( # 2 ( # , ( 2 # - 2 2( # - 2 ( 2 + # - 2
26
27 ! " = $%&&'( = ()*
28 ! $ "# = &'( ), + =, ), ), +, + =, )+, ) +, + ) +, ), + =, )+, ), +, ), + +, ), + =, )+, ), +! $ "" = &'( ), ) =, )), ), ) =, ) $, ) $ =! $ " = (/0 )
29 : cov(x,y) = cov(y,x) cov(x,x) = var(x) cov(ax, by) = ab*cov(x,y) cov(a+x, b+y) = cov(x,y) cov(x, Y+Z) = cov(x,y)+cov(x,z) cov(σx, ΣY)=ΣΣcov(X,Y)
30 " #$ = &'(( ), + = &',(), +),/( ),/((+)
31
32 ( "% " = $ %&' ( * +, = $ %&' ) ( /01 = $ %&' " % ", ) 1 /01(", 2) 4 = * + * 8 " % " 2 % 32 ) 1 9 = : " ; +, = : " 9, = : " : ", ; +8 = : " : " : 2 : 2 /01(", 2) < = 1=4 " 1=4(2)
33 ! " # = % ) *&,- # &'(. 1 ). 1! " # = % &'( ) = % &'( ) = % &'( ) = % &'( ) * &,- # = % &'( ) * & 1 # + % &'( ) * & + 1 1,- # = % &'( ),- 1 # 2,- 1 % &'( * & 1 # +.,- 1 # 2,- 1.,- 1 * & 1 #.,- 1 #. 1 3! # " = % 3 * & 1 #.3,- 1 # 3! " # = 4 " # ) &'( =.4 " #.4, 5 # =.4 # ". 4 " #. =. 1 4 # " * & + 1,- 1 # * & 1
34
35 ( ) = E ( X 1 ) + E ( X 2 ) +!+ E ( X n ) = µ + µ +!+ µ = nµ ( ) = var( X 1 ) + var( X 2 ) +!+ var( X n ) = σ 2 +σ 2 +!+σ 2 = nσ 2 E X 1 + X 2 +!+ X n var X 1 + X 2 +!+ X n ( ) = E X 1 + X 2 +!+ X # n E X var X! " ( ) = var! # " n X 1 + X 2 +!+ X n n $! & = E X $! 1 # &+!+ E X $ n # & = nµ % " n % " n % n = µ $ & = 1 % n var X 2 1 ( ) +!+ 1 n var ( X 2 n) = nσ 2 = σ 2 n 2 n
36 dat<-read.table(" > head(dat) shoesize h gender affil club M psy tetsudo M cs tennis M psy karate M psy tennis F cs tetsudo F psy tennis > summary(dat) shoesize h gender affil club Min. :21.00 Min. :145.0 F:36 cs :35 karate :24 1st Qu.: st Qu.:157.2 M:34 psy:35 tennis :23 Median :24.50 Median :164.0 tetsudo:23 Mean :24.82 Mean : rd Qu.: rd Qu.:170.9 Max. :28.00 Max. :182.3
37 > mean(dat$shoesize) [1] > mean(dat$h) [1] > var(dat$shoesize) [1] > var(dat$h) [1] > cov(dat[,1:2]) shoesize h shoesize h > cor(dat[,1:2]) shoesize h shoesize h
38 )! "# = % &'( ) "* & + & = "* ( + ( + "* "* ) + ) = " % * & + & = "! "# &'( > dat.meter = dat[,1:2]*0.01 > mean(dat.meter$h) [1] > mean(dat$h)*0.01 [1] > mean(dat$h) [1]
39 *! " + $ = & '() * + ', ' + $ &, ' = & + ', ' + $ - 1 =! " + $ '() * '() > dat.h.ext5 = dat$h+5 [1] > mean(dat.h.ext5) [1] > mean(dat$h) [1] > mean(dat$h)+5 [1]
40 *! " + $ = & + ', " = + ' + & - ', $ = - ' =! " +! $ '() * '() > handshoe = dat$h+dat$shoesize > mean(handshoe) [1] > mean(dat$h)+mean(dat$shoesize) [1]
41 ! "# + % =! "# +! % = "! # + % > dat.h.meter.ext5 = dat$h* > mean(dat.h.meter.ext5) [1] > 0.01*mean(dat$h)+0.05 [1]
42 !"# $% = ' $% ' $% ) > var(dat.meter$h) [1] = ' $% $* ) = ' $ ) % * ) = $ ) ' % * ) = $ )!"# % > var(dat$h)*(0.01^2) [1]
43 !"# $ + & = ( $ + & ( $ + & = ( $ + & + + & * = ( $ + * =!"# $ * > var(dat.h.ext5) [1] > var(dat$h) [1]
44 !"# $% + ' = ) $% + ' ) $% + ' + > var(dat.h.meter.ext5) [1] > var(dat$h)*(0.01^2) [1] = ) $% + ' $, + ' + = ) $ % +, + = $ + ) %, + = $ +!"# %
45 !"# $ + & = ( $ + & * + + *, - > var(handshoe) [1] = ( $ + & - 2 $ + & * + + *, + * + + *, - = ( $ - + 2$& + & - + 2$* + 2$*, 2&* + 2&*, + * * + *, + *, - = ( $ - 2$* + + * & - 2&*, + *, - + 2$& 2$*, 2&* + + 2* + *, = ( $ * ( & *, - + 2/01 $, & > var(dat$h)+var(dat$shoesize)+2*cov(dat$h,dat$shoesize) [1]
46 !"# $, $ = ' $$ ' $ ' $ = ' $ ) ' $ ) = *+, $ > cov(dat$h,dat$h) [1] > var(dat$h) [1]
47 !"# $%, '( = * $% * $% '( * '( = * $'%( * $% '( * '( $% + * $% * '( = $'* %( $'* % * ( $'* % * ( + $'* % * ( = $' * %( * % * ( = $' -!"# %, ( > cov(dat.meter$h,dat.meter$shoesize) [1] > (0.01*0.01)*cov(dat$h,dat$shoesize) [1]
48 !"# $ + &, ( + ) = + $ + & + $ + & ( + ) + ( + ) = + $ + & ( + ) + $ + & ( + ) + ( + ) $ + & + + $ + & + ( + ). = + $ + & ( + ) + $ + & + ( + ) + $ + & + ( + ) + + $ + & + ( + ) = + $( + (& + $) + &) $ + + & ( + + ) = $( + (+ & + $+ ) + + &) $( (+ & $+ ) + & + ) = + &) + & + ) =!"# &, ) > cov(dat.h.ext5,dat.ss.ext1) [1] > cov(dat$h,dat$shoesize) [1]
49 !"# $, & + ( = * $ * $ & + ( * & + ( = * $ & + ( * $ & + ( * & + ( $ + * $ * & + (. = * $ & + ( * $ * & + ( * $ * & + ( + * $ * & + ( = * $ & + ( * $ * & + ( = * $& + $( * $ * & + * $ * ( = * $& * $ * & + * $( * $ * ( =!"# $, & +!"# $, ( > cov(dat$h,handshoe) [1] > cov(dat$h,dat$shoesize)+cov(dat$h,dat$h) [1]
50 !"# $1 + $2, )1 + )2 = + $1 + $2 + $1 + $2 )1 + )2 + )1 + )2 = + $1 + $2 )1 + )2 + $1 + $2 + )1 + )2 = + $1)2 + $1)2 + $2)1 + $2)2 + $1 + )1 + + $1 + )2 + + $ $2 + )2 = + $1)1 + $1 + )1 + + $1)2 + $1 + )2 + + $2)1 + $2 + )1 + + $2)2 + $2 + )2 =!"# $1, )1 +!"# $1, )2 +!"# $2, )1 +!"# $2, )2 =..!"# $ / + ) 0 / $ / =..!"# $ / + ) 0 =. 123 $ / +..!"# $ / + $ 0 / / 0 / / 04/ =. 123 $ / +..!"# $ / + $ 0 / / 05/
DAA04
# plot(x,y, ) plot(dat$shoesize, dat$h, main="relationship b/w shoesize and height, xlab = 'shoesize, ylab='height, pch=19, col='red ) Relationship b/w shoesize and height height 150 160 170 180 21 22
More informationDAA02
c(var1,var2,...,varn) > x x [1] 1 2 3 4 > x2 x2 [1] 1 2 3 4 5 6 7 8 c(var1,var2,...,varn) > y=c('a0','a1','b0','b1') > y [1] "a0" "a1" "b0" "b1 > z=c(x,y) > z [1] "1" "2"
More informationDAA12
Observed Data (Total variance) Predicted Data (prediction variance) Errors in Prediction (error variance) Shoesize 23 24 25 26 27 male female male mean female mean overall mean Shoesize 23 24 25 26 27
More informationDAA09
> summary(dat.lm1) Call: lm(formula = sales ~ price, data = dat) Residuals: Min 1Q Median 3Q Max -55.719-19.270 4.212 16.143 73.454 Coefficients: Estimate Std. Error t value Pr(> t ) (Intercept) 237.1326
More information1 R Windows R 1.1 R The R project web R web Download [CRAN] CRAN Mirrors Japan Download and Install R [Windows 9
1 R 2007 8 19 1 Windows R 1.1 R The R project web http://www.r-project.org/ R web Download [CRAN] CRAN Mirrors Japan Download and Install R [Windows 95 and later ] [base] 2.5.1 R - 2.5.1 for Windows R
More information講義のーと : データ解析のための統計モデリング. 第2回
Title 講義のーと : データ解析のための統計モデリング Author(s) 久保, 拓弥 Issue Date 2008 Doc URL http://hdl.handle.net/2115/49477 Type learningobject Note この講義資料は, 著者のホームページ http://hosho.ees.hokudai.ac.jp/~kub ードできます Note(URL)http://hosho.ees.hokudai.ac.jp/~kubo/ce/EesLecture20
More information「統 計 数 学 3」
関数の使い方 1 関数と引数 関数の構造 関数名 ( 引数 1, 引数 2, 引数 3, ) 例 : マハラノビス距離を求める関数 mahalanobis(data,m,v) 引数名を指定して記述する場合 mahalanobis(x=data, center=m, cov=v) 2 関数についてのヘルプ 基本的な関数のヘルプの呼び出し? 関数名 例 :?mean 例 :?mahalanobis 指定できる引数を確認する関数
More informationplot type type= n text plot type= n text(x,y) iris 5 iris iris.label >iris.label<-rep(c(,, ),rep(50,3)) 2 13 >plot(iris[,1],iris
23 2 2 R iris 1 3 plot >plot(iris[,1],iris[,3]) 11 iris[, 3] 1 2 3 4 5 6 7 iris[, 3] 1 2 3 4 5 6 7 119 106 118 123 108 132 101 110 103 104 105 126 131 136 144 109 125 145121 130 135137 129 133 141 138
More information講義のーと : データ解析のための統計モデリング. 第5回
Title 講義のーと : データ解析のための統計モデリング Author(s) 久保, 拓弥 Issue Date 2008 Doc URL http://hdl.handle.net/2115/49477 Type learningobject Note この講義資料は, 著者のホームページ http://hosho.ees.hokudai.ac.jp/~kub ードできます Note(URL)http://hosho.ees.hokudai.ac.jp/~kubo/ce/EesLecture20
More informationRとExcelを用いた分布推定の実践例
R Excel 1 2 1 2 2011/11/09 ( IMI) R Excel 2011/11/09 1 / 12 (1) R Excel (2) ( IMI) R Excel 2011/11/09 2 / 12 R Excel R R > library(fitdistrplus) > x fitdist(x,"norm","mle")
More information- 1-128 - 2 -
127 - 1-128 - 2 - - 3-129 - 4 - 2-5 - 130-6 - - 7-131 - 8 - - 9-132 - 10 - 6041 3 () 1 ( ) () 6041 (1010) 1041 (192) 1941 () 2 (1) (2) (3) () 3 1 1 () 4 2 () 5 1 2 3 4 () 6 () 7-11 - 133-12 - 134 135 136
More informationスライド 1
担当 : 田中冬彦 016 年 4 月 19 日 @ 統計モデリング 統計モデリング 第二回配布資料 文献 : A. J. Dobson and A. G. Barnett: An Introduction to Generalized Linear Models 3rd ed., CRC Press. 配布資料の PDF は以下からも DL できます. 短縮 URL http://tinyurl.com/lxb7kb8
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<4D F736F F F696E74202D BD95CF97CA89F090CD F6489F18B4195AA90CD816A>
主な多変量解析 9. 多変量解析 1 ( 重回帰分析 ) 目的変数 量的 説明変数 質的 あり量的 重回帰分析 数量化 Ⅰ 類 質的 判別分析 数量化 Ⅱ 類 なし 主成分分析因子分析多次元尺度構成法 数量化 Ⅲ 類数量化 Ⅳ 類 その他 クラスタ分析共分散構造分析 説明変数 : 独立変数 予測変数 目的変数 : 従属変数 基準変数 3 1. 単回帰分析各データの構造 y b ax a α: 1,,,
More informationバイオインフォマティクス特論12
藤 博幸 事後予測分布 パラメータの事後分布に従って モデルがどんなデータを期待するかを予測する 予測分布が観測されたデータと 致するかを確認することで モデルの適切さを確認できる 前回と同じ問題で事後予測を う 3-1-1. 個 差を考えない場合 3-1-2. 完全な個 差を考える場合 3-1-3. 構造化された個 差を考える場合 ベイズ統計で実践モデリング 10.1 個 差を考えない場合 第 10
More informationZNR ( ) 8/20 µs 8/20 µs (A) (V) ACrms (V) C (V) max.(v) Ip (A) /0 µs 2 ms ERZV ERZV
ZNR ( ) ZNR V E R Z V φ φ φ φ φ φ V φ E R Z V A E H φ φ φ φ φ V 00 Sep. 20 ZNR ( ) ERZV05 @ 8/20 µs (J) @ 8/20 µs (A) (V) ACrms (V) C (V) max.(v) Ip (A) /0 µs 2 ms 40 0.6 0.4 2 25 ERZV07 36 2.5. 0.9 0
More information講義のーと : データ解析のための統計モデリング. 第3回
Title 講義のーと : データ解析のための統計モデリング Author(s) 久保, 拓弥 Issue Date 2008 Doc URL http://hdl.handle.net/2115/49477 Type learningobject Note この講義資料は, 著者のホームページ http://hosho.ees.hokudai.ac.jp/~kub ードできます Note(URL)http://hosho.ees.hokudai.ac.jp/~kubo/ce/EesLecture20
More information掲示用ヒート表 第34回 藤沢市長杯 2017
34 8 4 2 Round 1 Round 2 SEMI FINAL 30 16 8 H1 H5 H1 H1 Red 12401821 2 Red 12601360 2 1-1 Red 12501915 1 1-1 Red 12501915 4 White 12900051 4 White 12600138 3 3-1 White 12802412 2 3-1 White 12801091 1 Yellow
More information5 LATEX 2ε 2010
2010-11-27 15:30 16:00 TEX 5 LATEX 2ε 2010 1986 Lisp-Stat 1996 ptex 1987 ASCII TEX 1990 ptex 1993 JIS X 4051 1994 ptex JIS 1995 ptex 3.0 platex 2ε 2000 jsarticle 2008 ε-ptex e-ptex 2010 TEX Live 2010
More informationPackageSoft/R-033U.tex (2018/March) R:
................................................................................ R: 2018 3 29................................................................................ R AI R https://cran.r-project.org/doc/contrib/manuals-jp/r-intro-170.jp.pdf
More information2009 SEP. No.664 9 18 100460345 209710798 %0.11 3350955 9750799 12.1 0.93 100350 19 100600453 209700800 %0.03 3200065 9740873 11.0 0.98 90490 20 100750364 209680396 %0.81 321 844 978 591 11.6
More information宿題の解答
1 2006 1 17 prob=0.4 size=700 x 8.1 > x plot(x, pbinom(x, size=700, prob=0.4)) > abline(h=c(.025,.975), lty=2) 0.025 L = 254 0.025 U = 305 > pbinom(254, size=700, prob=0.4) [1] 0.02410222
More informationMH MH 9.50 8.50 9.40 8.40 9.30 9.20 8.30 9.10 9.00 8.20 8.90 8.80 8.10 8.70 8.60 8.00 8.50 7.90 8.40 8.30 7.80 8.20 8.10 7.70 8.00 7.60 7.90 7.80 7.50 7.70 7.60 7.40 1 7.50 7.40 7.30 2 7.30 7.20 7.20
More information2 Part A B C A > B > C (0) 90, 69, 61, 68, 6, 77, 75, 20, 41, 34 (1) 8, 56, 16, 50, 43, 66, 44, 77, 55, 48 (2) 92, 74, 56, 81, 84, 86, 1, 27,
/ (1) (2) (3) ysawano@tmu.ac.jp (4) (0) (10) 11 (10) (a) (b) (c) (5) - - 11160939-11160939- - 1 2 Part 1. 1. 1. A B C A > B > C (0) 90, 69, 61, 68, 6, 77, 75, 20, 41, 34 (1) 8, 56, 16, 50, 43, 66, 44,
More informationMicrosoft PowerPoint - 統計科学研究所_R_重回帰分析_変数選択_2.ppt
重回帰分析 残差分析 変数選択 1 内容 重回帰分析 残差分析 歯の咬耗度データの分析 R で変数選択 ~ step 関数 ~ 2 重回帰分析と単回帰分析 体重を予測する問題 分析 1 身長 のみから体重を予測 分析 2 身長 と ウエスト の両方を用いて体重を予測 分析 1 と比べて大きな改善 体重 に関する推測では 身長 だけでは不十分 重回帰分析における問題 ~ モデルの構築 ~ 適切なモデルで分析しているか?
More informationkubostat2017c p (c) Poisson regression, a generalized linear model (GLM) : :
kubostat2017c p.1 2017 (c), a generalized linear model (GLM) : kubo@ees.hokudai.ac.jp http://goo.gl/76c4i 2017 11 14 : 2017 11 07 15:43 kubostat2017c (http://goo.gl/76c4i) 2017 (c) 2017 11 14 1 / 47 agenda
More information★結果★ 藤沢市長杯 掲示用ヒート表
AA 35 Round 1 8 4 Round 2 28 16 SEMI FINAL H1 H5 H1 H1 Red 12802015 1 Red 12802109 1 1-1 Red 12802015 2 1-1 Red 12702346 White 12800232 2 White 12702406 3 3-1 White 12702346 1 3-1 White 12802109 Yellow
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- 111 -
- 110 - - 111 - - 112 - - 113 - - 114 - - 115 - - 116 - - 117 - - 118 - - 119 - - 120 - - 121 - - 122 - - 123 - - 124 - - 125 - - 126 - - 127 - - 128 - - 129 - - 130 - - 131 - - 132 - - 133 - - 134 - -
More information2 3 2
1 2 3 2 ha ha ha 3 4 ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 5 5 ( ) ( ) 3 ( ) ( ),,, 5 6 7 ( ) ( ) ( ) ( ) ( ) ( ) 6 36.3% ( ) 8 ( ) ( ) ( ) ( ) ( ) ( ) 3 5 ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )
More information300 10
300 10 16 1695 4 2 44 7 3 1 2 2 1869 11 192210 506.75ha 3 (1594) (1596) 16 17 7 1630 2 14 1637 16 17 4 1664 5 6 1678 6 1678 6 3 8 16954 8 (1758) 7 (1770) 17 2 (1745) 2 (1765) 4 11 11 10 12 1 12 12 12 2
More informationuntitled
2 23 3 28 10 12 1 2 3 18 6 1 1 2 9 9 5 22 1 22 10 2 100 9 9 9 1 11 2,300 22 10 8 90 5 80 90 100 5,000 100 90 5 90 8 90 13 1 22 20 / 14 / 10 / 1 13 80 8 80 80 7 8 13 80 7 90 2 200 300 3 4 10 1 1 90 5 90
More information14 12 ( 17 4 2 ) 2004
14 12 ( 17 4 2 ) 2004 1 4,700 1980 1,757 1,150 25,000ha 6,000ha 10 10 10 2 1981 56 90 4,700 1980 1,757 1,150 2003 187,614 96,940 51.7 17 102,128 54.7 2004 545 76 14 2004 4 13 49.1 1981 84,638 38,272 17,570
More information0
... 1... 2... 3... 5... 10... 13... 15... 19... 22... 24... 26... 27 0 0 1 2 H13 12 3 4 H21 15.7 5 29.8 10 1 6 4 1 10 1 50,284 48,834 46,185 41,056 33,892 26,726 7 8 9 1012 10 11 12 136ha 80.7ha 13 14
More informationuntitled
1,985ha 41.5 2,799ha 58.5 4,784ha 100.0 3 14,210 0 0.0 27 57,590 37,140 64.5 3 1,800 930 51.7 6 1,150 1,150 31.0 39 77,310 39,220 50.7 34 8-3 34 8-4 52 1-6 52 2-3 53 1-6 56 1-2 56 1 18 12 17 30 50
More information9 5 1 9 30 17 10 1 21 30 18 10 1 22 30 18 10 1 20 30 16 10 1 20 30 31 16 60 9 63 9 65 9 3 14 1 2 3 4 19 10 1 5,600 20,000 21,500 10 10 1,540 1,540 950 950 178 15 5,329 19 1,900 200 19 12 1 200 3,429
More information0.表紙
60 10 11 11 11,600 36,500 10,500 7,900 +1,500+1,500+9,000+18,0003,70041,600 12 3.84ha 3,100 3,100 4,400 2,785 5,000 50 70 3,500 3,500 2,600 30 2,600 45 3,375 40 4,500 20 2,000 12 1,500 12 2,300 7,500
More information1 4 8 19 65 87 108 132 142 156 - 1 - - 2 - - 3 - - 4 - - 5 - - 6 - - 7 - - 8 - - 9 - - 10 - - 11 - - 12 - - 13 - - 14 - - 15 - - 16 - - 17 - - 18 - 30-19 - - 20 - ( ) - 21 - - 22 - - 23 - - 24 -
More information,
34 29 5 29 6 29 7 30 28 26 1 1-71 - , - 72 - 27 3 1, - 73 - 5811 5821 5822 5831 58 2 5841 5851 5861 5862 5863 5864 5892 5896 5897 5898 5891 5895 7711 622 6311 6312 6314 6324 6032 6033 7611 7621 7623 7624
More information-2-
-1- -2- -3-24 32ha 13-4- 30-5- -6- -7- 100-8- 20 20 100-9- 100-10- 100 100-11- -12-30 -13- 10-14- 52 71 59-15- 30-16- -17- -18- 31 31-19- -20- 47 52 71 29 41 59 306 88 348-21- -22- -23- -24- 12-25- 1,300
More information総務委員会会議録
29 28 29 29 () 121 1 111 11928 3 () () () () 11928 3 1,320 32,7897,000 () 1 23725,000 9 1,320 313 2,7897,000 2 2 2 5.5 16 15200 31 29 2931 3110 21 21 5 6 8,2195,000 2134,1097,000 1,320 32,7897,000 31 ()
More information100 26331 26331 7,082ha(6.9% 9,027ha(8.8) 903ha(0.9%) 1,2292.5% 5,45711.2% 2,916(6.0%) 80,622ha79.0 4,507ha4.4% 18,277(37.4) 2949 (42.9% ( ) !! 1008940 (03)3581-4111 http://www.mof.go.jp/ 0608579
More information- 8 -. 1 (1) (2) (3) (4) (5) (6) (7) (8)
- 7 - (1) (2).. ha 43 ( (1) ( ) (2).. 18 11 30 43 (1).. (2).. 0.5. () () A) B) C) - 8 -. 1 (1) (2) (3) (4) (5) (6) (7) (8) (1) (2) 43-9 - - 10 - (1).. 43.. 43... 43 - 11 - (2) (1).. A B A A B A - 12 -.
More informationPowerPoint Presentation
R で自動作図 どんな絵でも, 何十回でも, 何百枚でも 竹中明夫 ( 国立環境研 ) プログラミング環境 + 実行環境 統計解析関数群 描画関数群 プログラム作図のすすめまずは描いてみる高水準作図関数と低水準作図関数繰り返し作業をプログラムで画像の保存 この話の目指すところ R のプログラムで絵が描けるよう になる R のプログラムで絵が描けそうな 気になる R のプログラムで絵を描きたくなる やる気にさえなれば,
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 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²ó
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 informationhttps://sites.google.com/site/prj311/ 2011-03- 11 14:46:18 ##### Ruby ##### h = Hash.new(0) while line = gets() if line =~ /^.*?\x01(.*?)\x01/ h[$1] += 1 end end h.sort_by { k,v - v }.each { k,v print
More information! & # # w w w w w w w w l & w_ # w_ w # w w w # w w # w w # w w w w bw w bw w bw w w bw w b w w_ l !!!!!! 6!!!! 6 ' ' ' ' ' ' ' ' ' ' ' '! ' ' ' ' ' ' ' ' ' ' ' ' ' '! ' ' ' ' ' ' ' ' ' ' '
More informationMicrosoft PowerPoint - R-intro-02.ppt
R で学ぶデータ解析とシミュレーション 2 ~ グラフ作成入門 ~ 2 時間目のメニュー グラフの作成方法 グラフ作成の第一歩 高水準作図関数 高水準作図関数の種類 関数 plot() を用いた作図例 低水準作図関数 低水準作図関数の種類 低水準作図関数を用いた作図例 数学関数のプロット 数学関数の定義方法 数学関数の作図例 参考 ( 重ねた図の描き方,R の画像の編集方法 ) 2 グラフィックスは
More information一般化線形 (混合) モデル (2) - ロジスティック回帰と GLMM
.. ( ) (2) GLMM kubo@ees.hokudai.ac.jp I http://goo.gl/rrhzey 2013 08 27 : 2013 08 27 08:29 kubostat2013ou2 (http://goo.gl/rrhzey) ( ) (2) 2013 08 27 1 / 74 I.1 N k.2 binomial distribution logit link function.3.4!
More information第32回新春波乗り大会2018
AA 32 Round 1 4 SEMI FINAL 2 20 8 FINAL H1 H1 H1 Red 12701793 1 1-1 Red 12701793 2 1-1 Red 11800623 White 12900058 4 3-1 White 12402115 4 2-1 White 12402209 Yellow 11603976 3 2-2 Yellow 12301534 3 1-2
More information( )/2 hara/lectures/lectures-j.html 2, {H} {T } S = {H, T } {(H, H), (H, T )} {(H, T ), (T, T )} {(H, H), (T, T )} {1
( )/2 http://www2.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html 1 2011 ( )/2 2 2011 4 1 2 1.1 1 2 1 2 3 4 5 1.1.1 sample space S S = {H, T } H T T H S = {(H, H), (H, T ), (T, H), (T, T )} (T, H) S
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 informationk3 ( :07 ) 2 (A) k = 1 (B) k = 7 y x x 1 (k2)?? x y (A) GLM (k
2012 11 01 k3 (2012-10-24 14:07 ) 1 6 3 (2012 11 01 k3) kubo@ees.hokudai.ac.jp web http://goo.gl/wijx2 web http://goo.gl/ufq2 1 3 2 : 4 3 AIC 6 4 7 5 8 6 : 9 7 11 8 12 8.1 (1)........ 13 8.2 (2) χ 2....................
More information3 M=8.4 M=3 M=.8 M=4.7 M=5.6 M=3 M=5. M=4.6 M=7 M=3 M= (interaction) 4 - A - B (main effect) - A B (interaction)
1 (two-way ANOVA) - - A B 1 3 M=8.4 M=3 M=.8 M=4.7 M=5.6 M=3 M=5. M=4.6 M=7 M=3 M=4 - - 1 (interaction) 4 - A - B (main effect) - A B (interaction) two-way ANOVA 5 1 A - H0: µ A 0 = µ A 1 = = µ A n - H1:
More informationバイオインフォマティクス特論4
藤 博幸 1-3-1. ピアソン相関係数 1-3-2. 致性のカッパ係数 1-3-3. 時系列データにおける変化検出 ベイズ統計で実践モデリング 5.1 ピアソン係数 第 5 章データ解析の例 データは n ペアの独 な観測値の対例 : 特定の薬剤の投与量と投与から t 時間後の注 する遺伝 の発現量 2 つの変数間の線形の関係性はピアソンの積率相関係数 r で表現される t 時間後の注 する遺伝
More information31 gh gw
30 31 gh gw 32 33 1406 1421 640 0 (mm) (mm) MAX1513 MIN349 MIN280 MAX900 gh gw 34 gh gh gw gw gh gh gw gw gh gh gw gw 35 175 176 177 178 179 180 181 195 196 197 198 202 203 2 1 L L L2 L2 L2 L 2 2 1 L L
More informationyamadaiR(cEFA).pdf
R 2012/10/05 Kosugi,E.Koji (Yamadai.R) Categorical Factor Analysis by using R 2012/10/05 1 / 9 Why we use... 3 5 Kosugi,E.Koji (Yamadai.R) Categorical Factor Analysis by using R 2012/10/05 2 / 9 FA vs
More information統計研修R分散分析(追加).indd
http://cse.niaes.affrc.go.jp/minaka/r/r-top.html > mm mm TRT DATA 1 DM1 2537 2 DM1 2069 3 DM1 2104 4 DM1 1797 5 DM2 3366 6 DM2 2591 7 DM2 2211 8
More information1 2 S-PLUS python MeCab csv S-PLUS csv S- PLUS S-PLUS csv python word1,word2,word3 2 python MeCab MeCab ipadic html 2010/3/26
2 1 2 2 2 3 3 4 4 4.1 2........................... 4 4.2........................ 5 4.3........................... 6 5 7 6 8 A 9 B 10 C 10 D 11 E 12 1 1 2 http://app.xrea.jp/dat/ S-PLUS python MeCab csv
More informationスライド 1
担当 : 田中冬彦 2017 年 4 月 18 日 @ 統計モデリング 統計モデリング 第二回配布資料 文献 : A. J. Dobson and A. G. Barnett: An Introduction to Generalized Linear Models 3rd ed., CRC Press. 配布資料の PDF は以下からも DL できます. 短縮 URL http://tinyurl.com/lxb7kb8
More information第122号.indd
-1- -2- -3- 0852-36-5150 0852-36-5163-4- -5- -6- -7- 1st 1-1 1-2 1-3 1-4 1-5 -8- 2nd M2 E2 D2 J2 C2-9- 3rd M3 E3 D3 J3 C3-10- 4th M4 E4 D4 J4 C4-11- -12- M5 E5 J5 D5 C5 5th -13- -14- NEWS NEWS -15- NEWS
More informationR-introduction.R
による統計解析 三中信宏 minaka@affrc.go.jp http://leeswijzer.org 305-8604 茨城県つくば市観音台 3-1-3 国立研究開発法人農業 食品産業技術総合研究機構農業環境変動研究センター統計モデル解析ユニット専門員 租界 R の門前にて : 統計言語 R との極私的格闘記録 http://leeswijzer.org/r/r-top.html 教科書と参考書
More information1.2 R R Windows, Macintosh, Linux(Unix) Windows Mac R Linux redhat, debian, vinelinux ( ) RjpWiki ( RjpWiki Wiki
R 2005 9 12 ( ) 1 R 1.1 R R R S-PLUS( ) S version 4( ) S (AT&T Richard A. Becker, John M. Chambers, and Allan R. Wilks ) S S R R S ( ) S GUI( ) ( ) R R R R http://stat.sm.u-tokai.ac.jp/ yama/r/ R yamamoto@sm.u-tokai.ac.jp
More informationA4“û”q
Contents 02 Q Q Q Q 04 Interview A A A A Q Q Q Q 05 Interview A A A A Q Q Q Q Interview A A A A 06 1 2 08 09 3 10 4 11 12 5 6 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 14 15 2 16 2
More informationHA TO MARK
HA TO MARK 02 03 04 2015 75 50 25-3.1pt +0.7pt 24.4% 24.2% 28.6% 18.3% 15.2% 23.0% 21.0% 17.9% 22.6% 23.3% +2.4pt 59.1% 61.5% 52.2% 54.8% 53.5% 0 2009 2010 2011 2012 2013 2009 2010 2011 2012 2013 2009
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 information¸×ÌÞÆ�°½-No.94-’ÓŠ¹
St. Paul s Kanko Club NEWS (2006) No.94 St. Paul s Kanko Club NEWS (2006) No.94 St. Paul s Kanko Club NEWS (2006) No.94 St. Paul s Kanko Club NEWS (2006) No.94 St. Paul s Kanko Club NEWS (2006) No.94 St.
More informationスライド 1
担当 : 田中冬彦 2018 年 4 月 17 日 @ 統計モデリング 統計モデリング 第二回配布資料 文献 : A. J. Dobson and A. G. Barnett: An Introduction to Generalized Linear Models 3rd ed., CRC Press. 配布資料の PDF は以下からも DL できます. 短縮 URL http://tinyurl.com/lxb7kb8
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 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 information_FH28_J.qxd
23.9 + + ø ø ø ø ø ø ø ø ø ++ + ++ + + + ± FH 28 D - 5 (25) S B -.5 SH (5) 2 3 4 5 6 7 8 9 2 6 7 3 8 4 9 5 ** ** ** ** ** ** ** ** ** ** ** ** ** CL586-86-4-** CL586-868-3-** CL586-823-5-** CL586-835-4-**
More informationIsogai, T., Building a dynamic correlation network for fat-tailed financial asset returns, Applied Network Science (7):-24, 206,
H28. (TMU) 206 8 29 / 34 2 3 4 5 6 Isogai, T., Building a dynamic correlation network for fat-tailed financial asset returns, Applied Network Science (7):-24, 206, http://link.springer.com/article/0.007/s409-06-0008-x
More information.....Z...^.[.......\..
15 10 16 42 55 55 56 60 62 199310 1995 134 10 8 15 1 13 1311 a s d f 141412 2 g h j 376104 3 104102 232 4 5 51 30 53 27 36 6 Y 7 8 9 10 8686 86 11 1310 15 12 Z 13 14 15 16 102193 23 1712 60 27 17 18 Z
More information