初心者のための R および RjpWiki に関する質問コーナー

新規投稿はできません



丸め方向の指定

abc (2004-11-30 (火) 19:25:46)

windows版のRを使用しているのですが、丸めモードの変換を行い、上向けに丸めたときの結果と下向きに求めた時の結果の比較を行いたいのですが、
どのようにすればよろしいのでしょうか?教えてください。
よろしくお願いします。

エクセルへ値コピー

ABA (2004-11-29 (月) 23:47:56)

Rで数値データを一列に出力させたとき、値の横に一行目には[1,]、
2行目には[2,]・・・と表示されるのですが、出力されたデータの値のみを
エクセルにコピー&ペーストするには、どうすればよいのでしょうか。
色々と試行錯誤したのですが、分からないので教えてください。
宜しくお願いします。

標準化回帰係数

B (2004-11-26 (金) 13:58:42)

重回帰をした際に標準化回帰係数を知りたいのですが、summaryでは表示されません。手計算で出すことができると思うのですが、出力する関数があれば教えてください

rimage インストール時のエラー

ささかま (2004-11-26 (金) 00:10:45)

R2.0.1 mac版を使っています。OS10.3.6 Xcode 1.1を入れてあります。 rimageをインストールするため、Finkを使ってfftwとlibjpegを入れました。その後R上からrimageをインストールしようとしましたが、

checking jpeglib.h usability... no
checking jpeglib.h presence... no
checking for jpeglib.h... no
configure: error: Sorry, can't find jpeglib header

と出てしまい、インストールできませんでした。fftwをcheckするときにも同様のメッセージが出ましたが、以前教えて頂いたように(「macでソースファイルからパッケージをインストールする方法」という題です)/usr/local/include/にシンボリックリンクを作ることでエラーは回避できました。しかし、jpeglib.hについては、シンボリックリンクを作ってもエラーが出てしまいます。別の場所にリンクを作らなければならないのでしょうか?どなたかご存知の方ご教授願います。 この掲示板の使い方がよくわからず上手く書き込めませんでしたが、どなたか直して頂いたようで、ありがとうございました。

統計的有意の記号

y.a. (2004-11-20 (土) 22:37:54)

例えば,回帰分析の結果に対して summary.lm()を signif.stars=T
オプション付きで実行するとt値の各有意水準に応じて星をつけてくれます.
この星の水準を変更するにはどうすればよいでしょうか?
具体的には,10%, 5%, 1% 水準でそれぞれ,*, **, ***
をつけたいのです(経済学ではこれがデフォルトです).

summary.lm をみると,係数の行列は ans$coefficients ですが,
これに星をつけている箇所はよく分かりませんでした.

ラベルの縦書き

学生 (2004-11-19 (金) 19:30:47)

日本語版R 2.00(Win版)を使用しています.

グラフを描くに,x軸やy軸のラベルを縦書きにしたいです.
以前,掲示板に少し書かれていますが,具体的に,どこをどう設定すれば,実現するのかが分かりません.
どなたか,ご存知の方,ご助言頂ければ幸いです.

Decisiontreeはどうすればできるのでしょうか?

MI (2004-11-19 (金) 17:51:39)

RでDecisiontreeができると人づてに聞きました。
調べるとtreeのような感じですが、よくわかりません。
どなたか、ご存知ないでしょか?

Scatterplot3Dについて

みち (2004-11-15 (月) 16:12:57)

色々調べてみても分からなかったので、こちらに投稿させて頂きました。

1.Scatterlot3D で Type=h として、プロットした点から底面に垂直線を引かせた後、点を滑らかな線で繋ぐにはどうしたら良いのでしょうか。

2.Scatterlot3Dで軸ラベルを回転させたり移動させたりするにはどうすればよいのでしょうか。

宜しくお願いします。

複数条件下でのベクトルの取り出し

P (2004-11-13 (土) 18:14:19)

初歩的な質問と思いますが,あれこれと調べてどうしても分からなかったので,質問させていただきます.

例えば,dataという配列もしくはデータフレームがあったとして,A列が1,B列が2を満たすC列を取り出したい場合,
data[,3][data$A==1][data$B==2]
では上手く行かないのですが,どのようにすればよいのでしょうか.
ご教示いただければ幸いです.

macでソースファイルからパッケージをインストールする方法

ささかま (2004-11-12 (金) 20:33:39)

R2.0.0 mac版を使用しています。OSは 10.3.6です。xcode1.0を入れてあります。
lpSolveというパッケージを導入したいのですが、CRANにはmac用のバイナリはありません。Package Installer を使ってソースからインストールしようとしましたが、
In file included from commonlib.c:17:
commonlib.h:6:20: malloc.h: No such file or directory
make: *** [commonlib.o] Error 1
ERROR: compilation failed for package 'lpSolve'
というエラーが出てしまい、出来ませんでした。
ソースを手動でダウンロードし、ターミナルから
R CMD INSTALL lpsolve
と実行しても同様のエラーが出てしまいました。ルートになってもだめでした。
Rを使い始めて数週間ですし、UNIXは素人です。色々と検索しては見たのですが、解決できませんでした。どなたか解決法をご存知の方、教えていただけないでしょうか?

RExcelのインストール

NK (2004-11-10 (水) 14:44:56)

R(D)COM Server V1.2を今まで使っていて、V1.2をアンインストールして残っていた関係したファイルも全て削除してからV1.35をインストールしましたが,インストールの途中で”C:WINNT?System32?msvcirt.dll Restart Replace failed Move FileEx Failed code 5"というエラーが出てしまい正常にインストールできません。どのように対応すれば良いでしょうか?宜しくお願いします。

出力画像のサイズを変更したい

fujim (2004-11-04 (木) 22:15:36)

初めまして。近頃Rをはじめようと思い、R-tips,データサイエンス入門
などを読みつついろいろと試しているものです。
環境:Windows2000,R2.0.0jp 実行方法:test.Rファイルを作成し.batファイルより R.exe CMD BATCH --slave --vanilla test.Rで実行
分からない点:例えばグラフィックス参考実例集barplotの項のbarplot03.R等において、出力画像.pngのサイズをどのように変更すればよいのでしょうか?
調べた点:parパラメータを以下のように変更してみた。
1) par(pin=c(5,10))
2) par(fin=c(5,10))
いずれも.outファイルにエラーが出て画像が出力できなくなります。
エラーは以下の通りです。
Error in plot.new() : Figure region too large
Execution halted
Warning message:
calling par(new=) with no plot
宜しくお願いします。

関数単位default width & heightdefault pointsize
bitmap()inch6‐(適当)
bmp()pixel48012
dev2bitmap()inch6
jpg()pixel48012
pdf()inch6
pictex()inchw=5, h=4なし
png()pixel48012
postscript()inch
quartz()inch512
windows()inch712
win.graph()inch712
win.metafile()inch712
win.print()inch712
x11()inch712
X11()inch712
xfig()inch

本に載せるいいネタを戴きました.ありがとうございます! -- 舟尾 2004-11-06 (土) 13:53:27

スクリプト実行

taku (2004-11-04 (木) 14:19:56)

perlやrubyの様に、一行目に次のようなRのパスを書いて、
一連のコマンドをスクリプトとして実行したいのですが、(ファイル名hello)

#!/usr/bin/R

ARGUMENT './hello' __ignored__
というメッセージが出て、普通にRのプロンプトが出てしまいます。

Rプログラムをスクリプトとして実行するにはどうしたらいいのでしょうか?

関数の中身を見る関数

akosuke (2004-10-29 (金) 10:45:01)

自分で作った関数で名前はわかるけど中身を忘れてしまった、、、というような場合に関数の中身を見るための関数はありますでしょうか?
ここでも調べてみたのですが調べ方が悪いのか見つかりませんでした。
くだらない質問で申し訳ありませんがどうぞよろしくお願いします。

AICの計算

AIC (2004-10-26 (火) 21:50:37)

The R Bookの第5章を読んでいる者です。p.109の上から2行目で「ARIMA分析で行ったようにプログラムを組む・・・」と書いてあるのですが、このARIMAモデルではq=0としてpのみ動かすという形にしていますよね。ところがGARCHでパラメタp,qを2つ動かしてAICを最小にするp,qを求めるにはどのようにプログラムを組んだら良いのかが分からないので、教えていただけると幸いです。
また、AICのパラメータ数はこの場合npram<-2と定義していますが、GARCHの
場合、式(5.14)にあるようにω、α、βの3と定義するべきではないのでしょうか?
それから(長くてすみません)GARCHの派生モデルを"fSeries"というパッケジで分析できると聞き、garchと同じ要領で書けばよい(例えばGJRなら
gjr(data,order(p,q,γ))??)と思ってやってみたのですができません。
ご存知でしたら教えていただけないでしょうか。

0を埋める

(2004-10-25 (月) 18:52:25)

とってもFAQだと思うのですが…
1 -> 001,42 -> 042
のように有効桁数(この例だと3ケタ)を決めて,
それに満たない場合は0で埋めるにはどうすればいいでしょうか?

変数名に別の変数の内容を含めたい

質問者 (2004-10-23 (土) 15:43:00)

大変初歩的な質問と思いますが、うまくいかないのでよろしくお願いします。

month <- c("Jan","Feb","Mar","Apr","May","Jun",
           "Jul","Aug","Sep","Oct","Nov","Dec")
for (i in 1:12) paste("Hoge.",month[i],".total",sep="") <- 0

として、例えばHoge.Jan.totalに0が格納されるようにしたいのですが、paste()の結果がクオートされているせいか、Error: Target of assignment expands to non-language objectとエラーになります。ご指導をお願いします。

データフレームの分割処理について

R-Fan (2004-10-18 (月) 14:38:21)

このサイトでRを知りました。いつも有益な情報をありがとうございます。
表題についてですが、大きなデータをキー項目(下記例ではdate)ごとに分割して処理するプログラムを作っています。subset()によるデータ抽出は大変便利ですがデータが大きくなると処理にかかる時間も増えてしまいます。この処理を高速化したいと考えておりますが、データフレームをキー項目で分割するような組み込み関数はないでしょうか。それともデータが大きくなる場合は、前処理段階でデータファイルを分割しておくのが良いでしょうか。全体のデータサイズは20〜30万件で、取り出すデータは1000〜2000件です。

(案1)
for (key in sort(unique(maindata$date))){
 subdata <- subset(maindata, date == key)
 #subdataの処理
}

自己レスになりますが、subset()を使用するかわりにキー項目に対応する行の索引を作成することで処理時間を4割程度削減することができました。

(案2)
idx <- tapply(1:nrow(maindata), maindata$date,
              function(subidx){return (subidx)}) 
for (key in sort(unique(maindata$date))){
  subdata <- maindata[idx[[toString(key)]],]
  #subdataの処理
}

もっと根本的に改善する方法がありましたらご教授いただければ思います。 どうぞよろしくお願いいたします。

関数の読み込み

恐縮です (2004-10-17 (日) 00:24:05)

Rを使い始めて間もない者です。
Rで色々な関数をinstall.packages()でインストールするのですが、ワークスペースを閉じるたびにまたインストールし直しになります。
一度インストールした関数をそのまま使い続けるにはどう保存すればよいのでしょうか?宜しくお願いします。

積集合(INTERSECT)と和集合(UNION)の中間的なものはありますでしょうか?

撃墜王 (2004-10-14 (木) 15:39:55)

R言語には、INTERSECTとUNIONといったコマンドが存在します。
ただ、入力されるリストが非常に多く存在した場合に、「それら全てのINTERSECT」を取ると、最終的なアウトプットがほぼ0になってしまいます。
また、「それら全てのUNION」を取ると、当然ながらアウトプットが膨大なものになってしまいます。。。どちらを採択するにしても大変悩ましく感じています。
そこで、質問させて頂きたいのは、INTERSECTとUNIONの中間的のようなもの、もしくは考え方などがあったら是非アドバイスいただければと思います。

なお、自分の中のイメージとしては

入力されるリスト 100個
1リストあたりに含まれる変数名(Aomori,Miyagi,Hokkaido...etc) 100個
全てのINTERSECTを取った場合 最終的に残る変数名 3個(少なすぎる。。。)
全てのUNIONを取った場合、 最終的に残る変数名 1024個(多すぎる。。。)

そこで、単純なUNIONではなくて、1つのリストにしか含まれないような変数名は除くようにしたいと思っています。
具体的な例をさらにあげますと

list1 Kanagawa,Tokyo,Saitama
list2 Chiba,Kanagawa,Tokyo
list3 Tokyo,Irabaki,Gunma

Intersect = Tokyo,Union=Kanagawa,Tokyo,Saitama,Chiba,Ibaraki,Gunma
となってしまいますが、上のようなリストの場合Kanagawaはlist1,2に共通して含まれ、唯一list3には含まれないような形なので、これは採用するといったようなUnionの考え方とIntersectの考え方をいっしょくたにしたような方法と言うのは実現可能でしょうか。色々調べたのですが、よい情報がみつからなかったので、諸先輩方アドバイスの程どうぞよろしくお願いいたします。

CSVファイルのインポートで errorが。

竹内 (2004-10-10 (日) 19:42:29)

こんにちは。Rを使い始めてまだ5日目の初心者です。使用環境は Mac OS X(10.3)で Rの versionは 1.91(英語版)です。 R-Bookの p59にある read.csv("ファイル名 ")で CSVファイルを読み込ませようとしても次のようなメッセージが出てしまいます。

Error in file(file, "r") : unable to open connection
In addition: Warning message:
cannot open file `test.csv'

R-Bookや googleで一応は調べた(つもり)なのですが、 errorが何を意味しているかわかりません。どなたかご教授をお願いします。

frequencyについて

R初心者 (2004-10-07 (木) 23:49:37)

Rで株価の時系列分析をしようとしている者ですが、frequencyについて、年次のデータであれば1を、月次のデータであれば12とすればよいと分かるのですが、株価のように日次データでしかも日曜日の値がない場合、frequency=365とするとずれが生じて横軸の値が正確でなくなるのですが、どのように解決すればよいのでしょうか?

GARCHモデルのパラメータ推定

ピエロ (2004-10-07 (木) 22:18:58)

RでGARCHモデル分析する際、パラメータ推定は何の方法を用いているのでしょうか?
擬似最尤法とかいうやつでしょうか?

あと、Rでの時系列分析の細かな手法を扱っている本あるいはサイト等ございましたら、教えて頂けると幸いです。
宜しくお願い致します。

pairsのプロット範囲

ぷーすか (2004-10-07 (木) 20:17:03)

散布図行列で、プロット点をドットに変え、y=xの直線を書きもうとして、

pairs(x,panel=function(x,y){points(x,y,pch=".");abline(0,1,col=2)})

でうまくできました。
ここでさらに各図のプロット範囲を同じにしようとxlim、ylimの指定を追加しようと思いましたが、

pairs(x,xlim=c(-3,3), ylim=c(-3,3), panel=function(x,y){points(x,y,pch=".");abline(0,1,col=2)})

ではエラーが出てしまいます。
プロット範囲以外のカスタマイズはあきらめ、

pairs(x,xlim=c(-3,3), ylim=c(-3,3))

とすれば、warningは出ますが、一応プロット範囲は制御できました。

プロット内容とプロット範囲のカスタマイズを両立させる方法はないのでしょうか?

eclipseの使い方

(2004-10-04 (月) 19:18:27)

eclipse (eclipse-SDK-3.0-linux-gtk.zip)をダウンロードし、3.0 用 Plug-inを入れました。eclipseを起動させたのですが、はて、使い方が分かりません。Window->Open Perspective->Other->StatET perspectiveを選択してみましたが、新規プロジェクトでRを選ぶこともできず、最初から躓いています。Eclipseを用いたRの使用について分かりやすく解説したサイトはないでしょうか?

エクセルファイルのデータをRに取り込む方法は?

初心者 (2004-09-27 (月) 11:51:31)

Rの初心者ですが、どうやってエクセルファイルのデータをRに取り込むんでしょうか?よろしくお願いいしマース〜〜



ヒストグラムX軸目盛について。

tofu (2004-09-21 (火) 11:11:26)

度数分布をヒストグラムにしたときののX軸の目盛の数字、目盛り幅を変えるのはどのようにすればいいのでしょうか?区間の境界値をいれたいのですが…教えてください。

plotでドットだけを打ちたい

KabuTaro (2004-09-20 (月) 04:02:13)

plot(c(1,2,3),c(4,5,6),type="p")
とすると(1,4),(2,5),(3,6)に丸印でプロットされますが
これを単にドットだけにするにはどうしたらよいでしょうか?
pchを弄くっても形が変わるだけでうまくいきません。
よろしくお願いいたします。

package(gregmisc)のheatmap.2の引数設定について

ちょろべぇ (2004-09-14 (火) 18:57:31)

gregmiscパッケージのheatmap.2を用いてヒートマップを作図しております。
heatmap.2は様々なパラメータ設定が可能であり、自分としては、distfun,hclustfunを下記のように設定してHeatmapを書きたいと思っております。
(希望)
?distfun→距離関数だと思われるので、ユークリッド距離を計算したい
?hclustfun→average linkageで結合したい
?できればdivisiveではなくagglomerativeでクラスタしたい

ただ、

heatmap.2(data,distfun=dist(method="euclidian"),
          hclustfun=hclust(dist(method="euclidian"),method="average"),
          dendrogram="both")

のように設定しても、hclustfun部分でエラーが帰ってきてしまいます。(T_T)
かなり愚かな質問だとは思いますが、アドバイスありましたらどうぞよろしくお願いいたします。

Hclust <- function(d) hclust(d, method = "average", members=NULL)
heatmap.2(data, hclustfun=Hclust, dendrogram="both")
# もしくは一度に済ませるには
heatmap.2(data, hclustfun=function(d) hclust(d, method = "average", members=NULL), 
          dendrogram="both")

エラーバー

Kei (2004-09-11 (土) 16:54:52)

実験のデータをエラーバーで表現したいですが、R言語でどのふうすればいいのかよく分かりません。一応、自分で作ってみました。内容は以下のようです:

mybar <- function(X) {
  n <- nrow(X)
  df <- n-1
  X.mean <- apply(X,2,mean)
  X.var <- apply(X,2,var)
  d.error<-qt(0.975,df)*sqrt(X.var/n)
  d.error.x <- (0.2 + 1) * 1:5 - 1/2   # エラーバーを表示する位置(x座標)
  barplot(X.mean, names.arg=colnames(X), 
              ylim=c(0,ceiling(max(X.mean+(d.error)))),
              col=c("blue", "green", "red", "yellow"))
  arrows(d.error.x, X.mean, d.error.x, X.mean+ d.error, angle=90)
}
mybar(分析するデータ[,列:列])

軸の調整方についてちょっと質問したいですが、教えてください。

棒グラフの重ね合わせ

TETSU (2004-09-10 (金) 11:51:50)

2つの棒グラフをbarplot関数を用いて描画したのですが、これらのグラフを比較するために、積み上げではなく重ね合わせで描画したいと思っています。色々と調べてみたのですがうまい方法が見つかりませんでした。どなたか描画法をご存知でしょうか?

層別の相関値

HOJO (2004-09-10 (金) 11:39:45)

 欠損値を含む二つの変数について,相関値をcor(x,y,use=pairwise complete observations)で計算しました.そうしたところ,全体では0.919だったのに,上位30件では,0.910,次の70件では,0.395, 残りの100件のデータでは,0.045となってしまいました.このようなことはどうして起こるのでしょうか.層別の方法に何か問題があるのでしょうか.

単純主効果

(2004-09-10 (金) 02:07:41)

2元配置の分散分析(完全無作為)をする場合に
交互作用がある場合,
1.「単純主効果の検定」,
2.「単純主効果の多重比較」
をしたいのですが,Rではどのよなプログラムを書けば良いのでしょうか.

参考までに, SASでの単純主効果についての情報と,そのデータを以下にリンクします.
例:http://koko15.hus.osaka-u.ac.jp/~kano/lecture/faq/interaction.doc
データ:http://mat.isc.chubu.ac.jp/R/tech.html#CRFpq

ご存知の方,教えて頂ければ幸いです

Rcmdrでグラフの変数選択欄

HOJO (2004-08-31 (火) 16:06:53)

たびたびですみません.RcmdrでIndexPlotなどを用いてグラフを作画しようとしても,変数選択欄に必要な変数名称が全ては現れず,作画することができません.またこの内容は作画すべきグラフの種類によっても異なりますがこれはどうしてでしょうか.なおPCは,W2Kを利用しています.

CSVファイルのインポート

HOJO (2004-08-31 (火) 11:26:52)

エクセルで作成した200行,20列程度の数値データのみのCSVファイルをRcmdrで読み込んだ時に次のエラーメッセージが出ます。RはV.1.9.0です.
1. Error: unable to open connection
2. data is not a data frame and cannot be attached
これはどのような点に気をつければよいのでしょうか。

非心パラメータ

(2004-08-24 (火) 19:12:06)

正確な用語がわからないのですが
SASではTNONCT関数で算出できる,
「t分布の非心パラメータ(引数は,x, df, prob)を求める関数」は,
Rにないでしょうか?
http://support.sas.com/91doc/getDoc/lrdict.hlp/a000245960.htm

例えば,
x <- 3.1; df<- 58; prob <- .975;
?(x,df,prob)

1.04844
と求めたいのです.

行列のグラフィカル表示image()の色に任意の値を割り当てるには

川村健介 (2004-08-18 (水) 13:14:01)

お世話になります。
行列をグラフィカルに表示するimage()関数の色調を任意に割り当て,さらにその色を凡例(棒)で表示するにはどうすればよいのでしょうか。
例えば配列の値が0〜0.88の間にあった場合,0.0〜0.5に白を割り当て,0.5〜0.9の間を0.1間隔でGreyの色調を与えたいのですが。

plot()で軸の数値をなくしたい

(2004-08-18 (水) 00:22:26)

10ほど並べた散布図で、最初の1つだけ軸の数値を表示し、残りの9個は軸の数値を省略したいとのですが、ラベルはxlab="",ylab=""でOKですが、数値の方はどうすればよいか分かりません。よろしくお願いします。

par(mfrow=c(2,5))
plot(sin)
plot(sin, xaxt = "n", yaxt = "n", xlab="", ylab="")
plot(sin, xaxt = "n", yaxt = "n", xlab="", ylab="")
# あと 7 回繰り返し

上のような感じで如何でしょうか?ご参考になれば幸いです.-- 舟尾 2004-08-18 (水) 00:40:00

  • ありがとうございます。目的の図を得ることができました。 -- 2004-08-18 (水) 11:21:11

Rで扱える統計手法

石川誠 (2004-08-16 (月) 01:26:14)

お世話になります。
Rで扱える統計手法を知りたいのですが、何処かに日本語で一覧のようなものはないでしょうか?

行列の値上位5つの要素の位置を知りたい

川村健介 (2004-08-14 (土) 09:05:51)

初めて投稿させていただきます。
630*630の配列の値の中から,上位5位までの要素の位置(x,y)をしるいい方法はありませんか?
which.maxだと最も高い値の位置,例えば39726と帰ってきます。
できれば,上位から[326,52],[555,156]...といった感じで位置を知りたいのですが。

statsライブラリ内のheatmap関数について

○田五郎 (2004-08-10 (火) 01:45:09)

こんばんわ。○田です。
この度はstasライブラリ内のheatmap関数について相談させていただきたく。
statsライブラリをロードして、heatmapがかけることはわかったのですが、
?heatmap内のデモスクリプトを書くなどしてみましたが、
heatmap関数の引数であるcol引数の種類が良くわかりません。。。
自分としては正の値が赤系、マイナスの値を青系で表示したいのですが
デモのcol = cm.colors(256)はレインボー系、、、
col = topo.colors(16)は緑と薄いこげ茶系と、自分の望んでいる色が出せないんですよね。。ググってもあまり良い情報も無く、もしみなさんどこか素敵なソース(情報源)等お知りでしたら、ご教授ください。
どうぞよろしくお願いしますm(_ _)m。

mvtnorm内で帰ってくる値が小数点4桁以下がその都度異なる?仕様?バグ?

ちょろちゃん (2004-08-08 (日) 15:34:16)

お世話になっております。
Rのmultcompライブラリのsimtest、method=dunnettを
利用しております。
同じデータセットを入力した際に、simtestのダネットの結果中、
T値は常に一定の値が帰ってくるのですが、P値の小数点4桁よりも小さな値は
計算するたびに異なっております。
どこで値が異なってしまったか追跡したところ、mulcompから呼び出されている
mvtnorm内でP値を計算するときにどうやら異なった値が帰ってきているようでした。
これは仕様なのでしょうか、それともmvtnormのバグなのでしょうか。
もしご存知の方がいらっしゃいましたらアドバイスいただければ幸いです。
どうぞよろしくお願いいたします。

繰り返しのt検定の結果を表示させるには

はまちゃん (2004-08-03 (火) 18:12:58)

お世話になっております。

a<-matrix(rnorm(96),c(16,6))
for (i in 2:6) t.test(a[,1],a[,i])

のように、繰り返しでt検定を行う場合に、結果が表示されません。
なにか方法があったら、教えてお願いできますか?
ほかの検定でも同じことが起こるのでしょうか?宜しくお願いします。


添付ファイル: filehm2.png 1627件 [詳細] filetest.png 2387件 [詳細] filecor.png 2433件 [詳細] filebp.png 2496件 [詳細] filescatterplot3d.png 1458件 [詳細] filetopoheat.png 2481件 [詳細] filehist3.png 2532件 [詳細] fileoutput.png 1352件 [詳細] filehist2.png 1354件 [詳細] filemultiplot.png 1231件 [詳細] filepairspanel.png 2435件 [詳細] filepchtest.png 2478件 [詳細] fileplot.png 2876件 [詳細] filescplot1.jpeg 2576件 [詳細]

トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2015-03-01 (日) 01:15:59