-1

些細なことのように聞こえますが、いくつかの研究では洗練された解決策が得られませんでした: カテゴリ変数 (GROUP) と連続読み取り変数 (血圧) で構造化されたデータフレームがあります。各グループの平均と標準偏差を示す簡単な箱ひげ図を作成するにはどうすればよいですか? 複数のグループがあります: A、B、C、D データフレーム内で ANOVA 事後分析を実行するにはどうすればよいですか。Mann-Whitney-U-Test ではどのように機能しますか? 棒グラフで有意水準をマークできますか? この操作を複数の連続変数 (dia_bloodpressure、sys_bloodpressure、mean_bloodpressure) に効率化し、異なるファイルの出力を (変数の名前で) シンク() するにはどうすればよいですか?

4

3 に答える 3

0

R でさらにプログラミングを行った後、医学研究に適した別の優れたパッケージ、psych を見つけました。上記の質問を考慮して、データフレームの統計的概要describe()をdescribeBy()取得し、グループ化変数で並べ替えます。この関数error.bars.by()は、平均値 +/- SD の高度なプロット関数です。このパッケージは共変量分析に関する多くの機能を提供します。これらは心理学的研究に役立ちますが、医療およびマーケティング研究にも役立つ可能性があります。

于 2013-01-05T16:42:23.960 に答える
0

いくつかの調査の後、私は agricolae パッケージを思いつきました。これは、複数のグループ比較を提供します。結果のオブジェクトは、グループごとの棒グラフ +/- SD または SEM の適切なプロット関数にパイプライン化できます。残念ながら、プロット内のグループ間で重要なマーカーを使用する方法はありません。

于 2012-09-17T17:40:40.693 に答える
0

考えられるコード スニペット:

library(psych)

x<-c(1,2,3,4,5,6,7,8,9,NA)
y<-c(2,3,NA,3,4,NA,2,3,NA,2)
group<-rep((factor(LETTERS[1:2])),5)
df<-data.frame(x,y,group)
df

by(df$x,df$group,summary)
by(df$x,df$group,mean)

sd(df$x) #result: NA
sd(df$x, na.rm=TRUE) #result: 2.738613

v = c("x", "y")#or
v = colnames(df)[1:2]
sapply(v, function(i) tapply(df[[i]], df$group, sd, na.rm=TRUE))

describeBy(df$x, df$group)

error.bars.by(df$x, df$group, bars=TRUE)
于 2013-01-05T17:08:42.723 に答える