![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
Rを使ったプログラミング †【トピックス】帰無仮説・有意水準・有意差・仮説検定。。。が分からない †前回授業でやった仮説検定のところがよくわからなかったという意見が多かったです。そこで、正確さは二の次にして、仮説検定で行っている作業をイメージで理解してもらえるようまとめてみました。それぞれどこの部分が、帰無仮説、有意水準、帰無仮説の棄却にあたるのか考えてみてください。また、なんとなく理解できたら、統計の教科書を読んでみることを強くお勧めします。
仮説検定：用語の整理： †
前回やった２項検定の補足：どうしてp値が0.04139になるか？ †> binom.test(15,20)
Exact binomial test
data: 15 and 20
number of successes = 15, number of trials = 20, p-value = 0.04139
alternative hypothesis: true probability of success is not equal to 0.5
95 percent confidence interval:
0.5089541 0.9134285
sample estimates:
probability of success
0.75
Rのおさらい †前回やったRを使った計算のおさらいを、サンプルデータをつかってやってみましょう。 サンプルデータ： 仲岡先生の実習から †下のデータは仲岡先生の実習でやった、調査コドラートの面積と、出現種数です コドラート1 10*10 2種 コドラート2 15*15 4種 コドラート3 20*20 7種 コドラート4 25*25 7種 コドラート5 30*30 7種 コドラート6 35*35 9種 コドラート7 40*40 8種 コドラート8 45*45 12種 コドラート9 50*50 10種 コドラート10 55*55 10種 コドラート11 60*60 11種 コドラート12 65*65 11種 コドラート13 70*70 9種 コドラート14 80*80 10種 コドラート15 100*100 12種 （データ提供、してくれた方に感謝！） Rを使ってこのデータの、面積と種数の相関関係を知りたいとき、必要なのは、コドラートの１辺の長さと種数のベクトルです。 今から何をするか、もうお分かりですよね？ 予想通り、K2Editorを立ち上げて、要らないところを消し、 10 2 15 4 . . というデータに整形して、エクセルにコピペして、１列ごとのデータを scan() 関数を使ってオブジェクトに読み込みます。では、やってみましょう。（そろそろミミタコだろうから、説明は簡単に書いておきます） 1. K2Editorでの操作 †
2. エクセルでの操作 †
3. Rでの操作 †
この授業では、それぞれの分析の意味については説明しません。仲岡先生に聞いて下さい。（グラフ用紙で提出しなくてもいいと仰っていたような、いなかったような。。。） ＊大事なこと＊ 「こんな解析、自分でやるのは絶対無理！」と思うかもしれませんが、 そんなことは無いです！ 上のような解析の方法、私だってRでできるとは知りませんでした。 では今回、私はどうやったのでしょう？「回帰分析行いなさい」なんて言われたら、どうすればいいのでしょうか？もしかしてRでできるかなと思ったら、 Adobe Readerを使って、マニュアルを検索します
仮説検定：補足 †前回アンケートで、仮説検定について「帰無仮説や有意水準などが分からない。」、「Rの操作は分かったけど、その意味が分からない」という意見があります。でも、私は、今の段階では、 全員がRを使って仮説検定の操作を行うことができただけで十分！ と考えています。生物学では収集したデータをもとに何らかの判断を下す場合、ほとんどの場合仮説検定が必要になります。これは別に系統学や生態学に限ったことではありません。最新号のNatureに載っている染色体不安定性を示すマウスの腫瘍に関する論文でもt検定が使われています。
Rを使ったプログラミング演習 †さて、今回で３週目になるRという計算ソフト。コマンドを入力するのはちょと難しいけど、エクセルみたいにマウスを使って面倒な操作をしなくても、一瞬にしてグラフが描けてしまう優れものですね。これなら、良いレポートを仕上げるために、何度も解析をやりなおして、グラフを描き直すっていう作業も面倒ではなさそう。。。 初めて使ったソフトなのに、皆さん、すでにRの操作には慣れてきているし（全員、課題ができていた）、興味も感じてくれたよう。そこで、今回は要望の多かったプログラミングをRを使って勉強します。 プログラミングって何？ †ところで、プログラミングって何をすることでしょうか？IT用語辞典で調べてみると、、、 載っていません 唯一近いところでは、「プログラミング言語」というのがあります。説明文からプログラミングに関わるところを抜き出して、「プログラミング」の説明文を作ってみると、
今日の授業でコンピュータにあれこれ命令を与えるために使うRは、計算ソフトですが、インタプリタ型のプログラミング言語ということもできます。 つまり、プログラミングとは、簡単に言うと、「言葉を使ってコンピュータが理解できる命令を作ること」 です。 コンピュータにどんな命令ができるのか？ or どんな命令がしたいか？ †「Wordを立ち上げて文章を作成し、印刷する」とか、「Excelで家計簿を管理する」とかは、市販のソフトウェアを使った作業ですが、いずれも、非常に高度で複雑な命令が、使用者からコンピュータに向かって発せられています。マウスを動かしたらポインタが動くなんていうのは、けっこう複雑な命令ですよ。でも、ほとんどの場合、一般の人がコンピュータでやりたいと思うほとんどの作業は、わざわざ「命令」なんていうものを意識する必要は無いはずです。利用者が理解しやすく・使いやすいように、専用のソフトウェアが開発され、販売されているでしょう。 生物学の研究者でも、プログラミングなどを経験したことの無い人は大勢います。だって、今や、DNAシーケンスの決定や、整列や、系統樹作成にだって、専用のソフトウェアがあるのですから。。。それぞれのソフトウェアの使い方を覚えれば、それですむわけです。私たちも新しい解析方法が開発されるたびに、そのソフトウェアの使い方を覚えるのに苦労しています。 でも、それだけで、本当にいいのでしょうか？良くないという点が、まず、２つ思いつきます。 1. 自分の目的にあったソフトウェアがいつも存在するというわけではない 2. ソフトが存在する場合、目的の数だけ、ソフトの使い方を覚えなければならない ということです。特に１番目は致命的ですね。ソフトウェアが開発されてなかったら、目的が達成できないのですから。例えば、集団遺伝学で学ぶ、遺伝的浮動の解析ソフトなんていうのは、パソコンショップに行っても売ってないです（でも、じつは、このくらいメジャーなものなら、ネットで探せばフリーのものがあります）。 そこで、この授業でプログラミングを勉強する理由は次の３点です。
初めて(?)のプログラミング："Hello World!" †「Hello World!」だなんて、なんか変なタイトルだなーと思いますよね。誰が使い始めたのかは知りませんが、様々なプログラミング言語の教科書で、最初に出てくるのがこのプログラミングです（プログラムを勉強したことのある人なら、たいてい誰でも知っています）。では、早速やってみましょう。Rを立ち上げて、次のように入力してください。 print("Hello World!")
うまくご挨拶できましたか？ ここでやったのは、print()という関数を使って"Hello World!"と画面に表示させただけですが、これも 画面にHello Worldと表示させる ということを目的としたプログラミングです。 えっ？それなら、Rの最初でやったオブジェクトの内容を画面に出すのもプログラミング？ と聞かれるかもしれませんが、その通りです。次の囲みの中を全てコピーして、Rのコンソールにペーストしてみましょう。 x=c(1,2,3,4,5,6,7,8,9,10) sum(x) 上でやったプログラミングは、 画面に1から10までの整数の合計を表示させる というプログラミングです。皆さんがこれまでにやってきたRを使った操作も、何らかの処理の結果を画面に表示させるプログラミングだったわけです。なんとなく、プログラミングが身近になったでしょ？ 複数行の命令文 †上でやった２つのプログラミングですが、違いがありますね。"Hello World!"の方は１行だけの命令だし、合計値の方は2行になっています。プログラムは通常、いろんな処理を組み合わせて作るので、複数行になることが多いです。次に３行からなるプログラミングをやってみましょう。 q=c(10,15,20,25,30,35,40,45,50,55,60,65,70,80,100) y=c(2,3,4,5,6,7,8,9,10,11,10,11,9,10,12) plot(log10(q^2), log10(y)) このプログラムは、 ２つの数値ベクトルについて計算を行って、対数グラフを描く というものでした。さっき、上でやったのと同じことを、複数行にまとめて書いて、まとめてRに命令したわけです。 Rで使う命令 †授業では一つ一つ解説しますが、復習課題をやるときとか、自分で新たな解析に挑戦するときは、どういう命令があるかが書かれた説明書が必要になるでしょう。
プログラミングの基礎: 繰り返しと代入と条件分岐 †上で複数の行で様々な命令を一度にＲに与えられるということを解説しました。でもやったことは、それぞれ個々の計算を順番に並べて、一度にペーストしただけですよね。そこで、次にもっとプログラムらしい命令をRに与えてみましょう。それは、繰り返しと代入と条件分岐というもので、プログラミングの基本中の基本です。 繰り返し †人が不得意でコンピュータが得意なことの一つは、単純な繰り返し作業を際限なく（文句も言わずに）行うことです。要するに、プログラミングの一つの目的は、そういう面倒な繰り返し計算をコンピュータにやらせることです。では、何回繰り返すかというのをどうやってコンピュータに指示するかというと、forという命令を使います。
代入 †上の for(i in 1:10) { print(i) } というプログラムでは、iの値は繰り返しの度に違う値になります。 1回目 iの値は1 2回目 iの値は2 3回目 iの値は3 .... ........ 10回目 iの値は10 では、一つ前の繰り返しの時のiの値を使いたい場合はどうすれば良いでしょうか？例えば、1から10までの合計をfor命令を使ってプログラミングすることを考えると、 今の繰り返し回数のiの値に１つ前回数までの合計を加えて表示させたい 1回目 iの値は1 １つ前までの合計は0 1 を表示 2回目 iの値は2 １つ前までの合計は1 3 を表示 3回目 iの値は3 １つ前までの合計は3 6 を表示 .... ........ 10回目 iの値は10 １つ前までの合計は45 55 を表示 こういう場合に代入という命令を使います。その命令は = です。通常の計算では'='というのは「3+5=」というように、計算を行わせる意味で用いられますが、プログラミングにおいては、 左辺の変数に右辺の計算結果を代入する という意味で用いられます。上の計算のプログラムは次のようになります。 goukei=0 #goukeiという変数に初期値0を代入
for(i in 1:10) { #以下を10回繰り返す
goukei=goukei+i #goukeiに前回までのgoukeiの値にiの値を足したもの代入
print(goukei) #goukeiの内容を画面に表示
}
goukei=goukei+iという命令文では、右辺の「goukei+i」を先に計算して左辺の「goukei」に代入しているところが、分かりにくいかもしれません。 練習問題： 1から10までの数字を１つずつ順々にかけ合わせたた結果を全て表示させるプログラムを作りなさい
1x1=1, 1x2=2, 2x3=6, 6x4=24, 24x5=...... と計算してゆくということ
条件分岐 †プログラミングの基本技の最後は１つは条件分岐です。if()を使って表現します i=3 #iに3という数値を代入
if(i==3){print("三")} #iの値が3ならば画面に三と表示
この例の２つめのif命令で、iの値を評価しています。評価に使われるのは比較演算子というもので、 == 等しければ != 等しく無ければ >, >= 左辺が右辺よりも大きければ、左辺が右辺以上ならば <, <= 左辺が右辺よりも小さければ、左辺が右辺以下ならば ということを意味しています
繰り返しと条件分岐 †では、上の繰り返しと条件分岐を組み合わせて見ましょう 練習問題：for命令を使って1から10までの数字を表示させなさい、5だけは"five"と英単語で表示させなさい
第10回授業の課題 †
課題１．アンケート調査 †
課題１．意見調査 †下の囲みの中にあるアンケートをコピー・ペーストして、「回答：」の後に答えを書き込むこと。 *第10回授業アンケート **氏名： **課題への回答 -今日（6月22日）の授業の進み方は？（はやい、丁度いい、おそい） --回答： -今日の授業の難しさはどう感じましたか（簡単すぎ 簡単 丁度いい 難しい 難しすぎ）： --回答： -難しいと答えた人は、特にどの点が難しかったですか？： --回答： -今日の授業は（よく分かった 分かった 分からなかった）： --回答： -分からないと答えた人は、特にどの点が分からなかったですか？： --回答： -今日の講義で理解できなかった用語があったら挙げてください： --回答： -Ｒで何ができるのか、だんだん分かってきたと思います。今後自分の勉強にＲを使いたいですか？ --回答： -Ｒでも、Ｒ以外の言語でも、何かやってみたいプログラミングはありますか？ --回答： -今回やったプログラミングはおもしろかったですか？ --回答： 復習課題：Rによるプログラミング：繰り返しと条件分岐 †次の問題文の指定にあるプログラムを作成しなさい。作成したプログラムは課題提出ページに貼り付けなさい。 行頭には必ず半角の空白を入れること。
予習課題：Rを使った関数の定義 †「関数の定義」なんていうと難しそうですが、ようするに、作ったプログラムに名前をつけて、いろいろと数値を変えて解析できるようにすることです。例えば、 1から入力した数値までの全てを横一列に表示させるプログラムを作り displayという名前の関数として定義する どうやるかというと、 display=function(a){ #関数定義の始まり
kekka=c() #kekkaに空ベクトルを代入して初期化
for(i in 1:a){ #a回（iの値を1からaまで変化させる間）繰り返し
kekka=append(kekka,i) #kekkaというベクトルにiを要素として代入
}
print(kekka) #kekkaの内容を表示
} #関数定義の終わり
上の方で作ったプログラムの上と下を関数をfunction(){}で囲んで、オブジェクトに代入するだけですから、簡単です。この関数を実行するには、最初aの値を()の中に入力して、 display(10) とします。 ↑を使って何回も実行してみると、結果がいろいろ変わるのがわかりますよね。 そこで問題です。上の関数の定義方法に従って、入力した数までの合計値を計算するsumupという名前の関数を作成します。下の囲みの中の_の部分（１文字に対応するとは限らない）を埋めて、プログラムを完成して提出しなさい。 sumup = ______ #関数sumpuを定義
kotae = _ #kotaeを初期化
_ (i in ___){ #1からaまで繰り返し
kotae = _____ #kotaeにiの値を足したものとkotaeに代入
} #繰り返し終了
_____ #kotaeを表示
} #関数定義の終了
|
||||||||||||||
期待値.xls 1709件
[詳細]
snmpenwl.gif 804件
[詳細]
snmpengi.gif 781件
[詳細]
snmpensl.gif 785件
[詳細]
snmpendf.gif 800件
[詳細]
snmpenas.gif 786件
[詳細]
WS000002.JPG 773件
[詳細]
pengraph.gif 1786件
[詳細]
Untitled-8.gif 1498件
[詳細]
Untitled-7.gif 1601件
[詳細]
![[smile]](image/face/smile.png)
)
