私は次のようなデータフレームを持っています:
case simulation temp plank oxygen
1 1 1 8 7 11
2 2 1 16 10 15
...
17 17 2 26 12 17
18 18 2 15 8 12
19 19 2 28 11 21
20 20 2 24 6 14
シミュレーション変数のレベルで要約を分割したいのですが。たとえばtemp
、simulation == 1とSimulation == 2 の平均と、標準偏差も同じです。
現在、私は次のコードを使用していますが、これはかなり恐ろしいことです。
df <- read.csv("data.csv")
attach(df)
# Create subset variables
temp1 = subset(temp, simulation==1)
temp2 = subset(temp, simulation==2)
plank1 = subset(plank, simulation==1)
plank2 = subset(plank, simulation==2)
oxygen1 = subset(oxygen, simulation==1)
oxygen2 = subset(oxygen, simulation==2)
print(sd(temp1))
print(sd(temp2))
print(sd(plank1))
print(sd(plank2))
これをRで自動的に行う方法があるはずですが、見つかりません。私はを使用してみましたsummary(df ~ simulation)
が、それは有用なものを何も生成しません。