K010 appstat201201
- 2. ガイダンス 2
? vB初
? 表寄僥 h廠僥冩梢親 縮娩
換邦 慌岻┐燭襪漾，箸發罎
e-mail tarumi@ems.okayama-u.ac.jp
? 娩Iの古勣
? y僥秘T
? yソフトウェア R を旋喘したシミュレ`ションgYを宥して
デ`タ盾裂のgHに乾れるとともにy尖を_かめる
? デ`タ盾裂に乾れる
? y尖を_かめる
- 3. ガイダンス 3
? 廣吭並
? 朕の念に PC があるが尖のh苧も謹い。医しているときは PC
に乾れずにh苧をくこと。
? テキスト
? 換邦?V广
仝 R/S-PLUS によるy盾裂秘T々
慌羨竃井 , 2500
? 歌深
? R によるやさしいy僥 表弥 雰 票 冷拭 翫小 匯隻 慌广
オ`ム芙
? 弥拐雰 (2004) 仝 The R BookD デ`タ盾裂h廠 R の試喘並箭鹿々
湘爺芙 .
? 巍硫劃仂 (2005) 仝 The R TipsD デ`タ盾裂h廠 R の児云室?グラ
フィックス試喘鹿々湘爺芙 .
? その麿 R S にvする云やy僥にvする云
- 4. 讐吶の坪否 4
? y喘 S 冱Z
? S, S-PLUS 斌喘ソフトウェア
? R フリ`ソフトウェア
? Web
? RjpWiki http://www.okada.jp.org/RWiki/
? R 塘下圷 http://www.r-project.org/
? vxサポ`トペ`ジ
? http://webclass.el.okayama-u.ac.jp/
? ID 寄僥原嚥?パスワ`ドでアクセス
? 鞠hキ` Rstat2012
? 歌深ペ`ジ
? http://case.f7.ems.okayama-u.ac.jp
- 5. R のインスト`ル (Window の栽 ) 5
? ダウンロ`ド
? http://cran.md.tsukuba.ac.jp/bin/windows/base/ CRAN mirror
? Download R 2. １ 5.0 for Windows (?? megabytes)
? 嘔のウィンドウが燕幣されたら
仝g佩々を
? この少御にも
仝g佩する々を
? 參瘁テキスト１７０を
歌孚のこと
- 6. 扮 6
? vxの扮
? マイドキュメント に appstat というフォルダを恬撹しよう
? この讐吶の坪否はこのフォルダで恬Iを佩う
? R のアイコン┘轡腥`トカットの俐屎
? R のアイコンを嘔クリック
? 仝プロパティ々をxk
? 恬Iフォルダを
＾ c:Documents and Settingst2My Documents"
を
＾ c:Documents and Settingst2My Documentsappstat＾
に筝して仝 Ok 々を兀す
- 7. 軟強嶮阻 7
? R の軟
? R のアイコンをダブルクリック
? R のK阻 (quit)
K阻の栽はワ`クスペ`スのイメ`ジを隠贋するかど
q() うかいてくる
軟咾靴討い謇g恬った篳v方などを隠贋したけれ
ば
y
ただgに弯としてまたは_Jとして聞っただけで
隠贋までしたくない栽は
n
- 8. demo 8
? R はグラフィックが掲械にれています。
? お返Xにyグラフが恬撹できます。
? これからのY創でコマンドの枠^に召秘ることがありま
すがこれはプロンプトを吭龍しています。タイプする駅勣
はありません。
? デモのg佩
>demo()
たとえば
? >demo(graphics)
- 9. 窮弯として旋喘 9
? 弯としても聞える
?
?
?
?/
?2 \３\は x^2, x^3
ル`トは
? sqrt(2)
など
? ┌などもつかえる
? ル`ルは宥械の麻と揖じ
? ┌枠
- 10. 旗秘 10
? x に 10 を旗秘する
x <- 10
x に採がはいっているか_Jするときは
x とだけ秘薦しよう。
? その瘁
x <- x+1
をしたらどうなるか深えてみよう
? 肝の麻をしよう
y <- 10
x <- 20
y+x y*x, y-x, y/x, x^2, x^3, sqrt(x), log(x)
をやってみよう
- 11. ベクトル柴麻 11
? ベクトル
? ベクトルを恬る児云議なv方は c( )
? x <- c(10,12,14,15,15)
? x には採が秘っているか。 x と秘薦してエンタ`を兀し_Jしてみよ
う。
? ベクトルのLさ
? length(x)
? 栽をだしたかったら
? sum(x) v方侏冱Z
v方兆 ( 哈方のリスト )
? x の２\を麻してみよう v方兆 () ? 哈方がoくても凄察
? x^2 は駅勣
? と秘薦する
? x の２\の才を麻してみよう
- 12. 佩双處麻 12
? > X<-matrix(c(1,2,3,4),2,2)
? >X
? とすれば佩双が恬られる .
寄猟忖と弌猟忖は曝eされていることに廣吭。
? X^2 としたときどういう､卦されるか_Jしよう
? Y<-matrix(c(4,5,6,1),2,2) としたとき
? X+Y, X-Y, X*Y, X/Y を_Jしよう
? XY を麻したいrは
? X%*%Y
? X*Y では勣殆揖平の譴泳磴砲覆
- 13. 13
? 剃佩双は
? solve(X)
? X%*%solve(X) を麻してみてg了佩双になっていることを_かめよう
。
? 佩双塀 determinant
? det(X)
? 崔佩双 transpose
? t(X)
? ここで匯業K阻してみよう。K阻の碧圭は
- 14. 14
? デ`タの扮
> height <- c(148, 160, 159, 153, 151, 140, 156, 137, 149, 160, 151, 157, 157,
144)
http://mikawaya.to/pukiwiki/
? _J
? ls()
? }u
? height2<-height
? 茅
- 15. デ`タの俐屎 15
? 恷瘁に弖紗
> x <- c(23, 22, 52, 13, 4)
>x
[1] 23 22 52 13 4
> newx <- c(x, c(10, 11, 12, 13))
> newx
[1] 23 22 52 13 4 10 11 12 13
> newxx <- append(x, c(10, 11, 12, 13))
> newxx
[1] 23 22 52 13 4 10 11 12 13
- 16. 余嶄に弖紗 16
> x <- c(1, 2, 3, 4, 5)
>x
[1] 1 2 3 4 5
> newx <- c(x[1:2], c(10, 11, 12, 13), x[3:length(x)])
> newx
[1] 1 2 10 11 12 13 3 4 5
> newxx <-append(x, c(10, 11, 12, 13), after = 2)
> newxx
[1] 1 2 10 11 12 13 3 4 5
- 17. 茅崔算俐屎 17
? m 桑朕のデ`タを茅
> x <- c(3, 3, 2, 3, 3)
>x
[1] 3 3 2 3 3
> newx <- x[-3] #x の 3 桑朕の､任△ 2 を茅し newx
に旗秘
> newx
[1] 3 3 3 3
? m 桑朕のデ`タの､ z に筝
> x <- 1:5
>x
[1] 1 2 3 4 5
> x[3] <- 100
>x
[1] 1 2 100 4 5
- 18. 勣殆の渇竃 18
? ベクトルの栽
? x<-c(10,12,14,15,15)
? x の２桑朕の､鰔たい
? x[2]
? で勣殆を函り竃せる
? 佩双の栽
? x<- matrix(c(1,2,3,4),2,2)
? x の２佩１双朕をたい
? x[2, 1]
? で勣殆をとりだせる
? ２佩朕畠何をとりだしたいときは
? x[2, ]
? 揖に２双朕
? x[, 2]
? x の３桑朕の､ y に旗秘したかったら
? y <- x[3]
- 19. g佩潤惚のファイル竃薦 19
? 參和の鮫中に竃る燕幣 ( グラフを茅く ) をファイルへ竃薦
>sink(＾ ファイル兆￣ )
> ??? 鮫中には採も燕幣されない
>sink()
> sink("result.txt")
result.txt の嶄附を_Jしよう
> height
> mean(height)
> max(height)
> sink()
- 20. デ`タの隠贋隠贋したデ`タ?v方のiみzみ 20
? 屡協x篳の隠贋
workspace を兆念をつけて隠贋する。
? save.image(＾ ファイル兆￣ )
? save.image(＾intro.Rdata￣)
? 隠贋したデ`タv方のiみzみ
? load("intro.RData")
? ls()
? なにも峺協せず隠贋した栽┐燭箸┐丕 q() でｙとした
とき .RData というファイルに隠贋されます。
? 肝指軟rは採も峺協がなければそのディレクトリにあ
.Rdata をiみzみます。なければ採も秘ってない彜Bになり
ます。
- 21. ベクトル柴麻 21
? 處
? x の２\の才を麻してみよう
? x+x はどういう､砲覆辰討い襪
x <- c(10,12,14,15,15)
y <- c(1, 2, 3 4, 5)
? としたときに
? x+y
がどのような､砲覆辰討い襪。
- 23. 竃薦デバイス 23
? TeX に函りzみやすい PS ファイルで竃薦するには
? 2F のg片の栽は
? postscript(＾***.eps", horizontal = FALSE, onefile = FALSE, paper =
"special", height = 9, width = 9, family=＾Courier￣) )
? plot(height, weight)
? dev.off()
? Windows の栽 Windows を軟咾靴憧_Jしよう。
![デ`タの俐屎 15
? 恷瘁に弖紗
> x <- c(23, 22, 52, 13, 4)
>x
[1] 23 22 52 13 4
> newx <- c(x, c(10, 11, 12, 13))
> newx
[1] 23 22 52 13 4 10 11 12 13
> newxx <- append(x, c(10, 11, 12, 13))
> newxx
[1] 23 22 52 13 4 10 11 12 13](https://image.slidesharecdn.com/k010appstat201201-130416203036-phpapp02/85/K010-appstat201201-15-320.jpg)
![余嶄に弖紗 16
> x <- c(1, 2, 3, 4, 5)
>x
[1] 1 2 3 4 5
> newx <- c(x[1:2], c(10, 11, 12, 13), x[3:length(x)])
> newx
[1] 1 2 10 11 12 13 3 4 5
> newxx <-append(x, c(10, 11, 12, 13), after = 2)
> newxx
[1] 1 2 10 11 12 13 3 4 5](https://image.slidesharecdn.com/k010appstat201201-130416203036-phpapp02/85/K010-appstat201201-16-320.jpg)
![茅崔算俐屎 17
? m 桑朕のデ`タを茅
> x <- c(3, 3, 2, 3, 3)
>x
[1] 3 3 2 3 3
> newx <- x[-3] #x の 3 桑朕の､任△ 2 を茅し newx
に旗秘
> newx
[1] 3 3 3 3
? m 桑朕のデ`タの､ z に筝
> x <- 1:5
>x
[1] 1 2 3 4 5
> x[3] <- 100
>x
[1] 1 2 100 4 5](https://image.slidesharecdn.com/k010appstat201201-130416203036-phpapp02/85/K010-appstat201201-17-320.jpg)
![勣殆の渇竃 18
? ベクトルの栽
? x<-c(10,12,14,15,15)
? x の２桑朕の､鰔たい
? x[2]
? で勣殆を函り竃せる
? 佩双の栽
? x<- matrix(c(1,2,3,4),2,2)
? x の２佩１双朕をたい
? x[2, 1]
? で勣殆をとりだせる
? ２佩朕畠何をとりだしたいときは
? x[2, ]
? 揖に２双朕
? x[, 2]
? x の３桑朕の､ y に旗秘したかったら
? y <- x[3]](https://image.slidesharecdn.com/k010appstat201201-130416203036-phpapp02/85/K010-appstat201201-18-320.jpg)
