次のスクリプトがあります。
myfunctionSD2 <- function(mydata) { return(sd(mydata,na.rm=TRUE))}
SDmeas <- tapply(datIn$Measurement,list(as.factor(datIn$Measurement.location),as.factor(datIn$Tube.number)),myfunctionSD2)
SDMeasurement <- SDmeas[order(as.numeric(rownames(SDmeas))), ]
dput(SDMeasurement)
structure(c(0.6555264303958, 0.605687762634627, NA, 0.683435292388886,
NA, 0.615207645374612, NA, 0.739018717912613, NA, 0.79341715750565,
NA, 0.769473075819677, NA, NA, NA, NA, NA, 0.6555264303958, 0.605687762634627,
0.576875822240553, NA, 0.502849025908516, NA, 0.516028792109233,
NA, 0.486023134629369, NA, 0.489281948774758, NA, 0.464851913610455,
NA, 0.499546896146173, NA, NA, 0.6555264303958, 0.605687762634627,
NA, 0.683435292388886, NA, 0.615207645374612, NA, 0.739018717912613,
NA, 0.79341715750565, NA, 0.769473075819677, NA, NA, NA, NA,
NA, 0.6555264303958, 0.605687762634627, NA, 0.683435292388886,
NA, 0.615207645374612, NA, 0.739018717912613, NA, 0.79341715750565,
NA, 0.769473075819677, NA, 0.773437871034822, NA, NA, NA, 0.6555264303958,
0.605687762634627, NA, 0.683435292388886, NA, 0.615207645374612,
NA, 0.739018717912613, NA, 0.79341715750565, NA, 0.769473075819677,
NA, 0.773437871034822, NA, 0.760395526989779, 0.607323723612999,
0.6555264303958, 0.605687762634627, NA, 0.683435292388886, NA,
0.615207645374612, NA, 0.739018717912613, NA, 0.79341715750565,
NA, 0.769473075819677, NA, 0.773437871034822, NA, 0.760395526989779,
NA), .Dim = c(17L, 6L), .Dimnames = list(c("1", "3", "4", "5",
"6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16",
"17", "19"), c("1", "2", "3", "4", "5", "6")))
明確な形式では、テーブルは次のようになります。
1 2 3 4 5 6
1 0.6555264 0.6555264 0.6555264 0.6555264 0.6555264 0.6555264
3 0.6056878 0.6056878 0.6056878 0.6056878 0.6056878 0.6056878
4 NA 0.5768758 NA NA NA NA
5 0.6834353 NA 0.6834353 0.6834353 0.6834353 0.6834353
6 NA 0.5028490 NA NA NA NA
7 0.6152076 NA 0.6152076 0.6152076 0.6152076 0.6152076
8 NA 0.5160288 NA NA NA NA
9 0.7390187 NA 0.7390187 0.7390187 0.7390187 0.7390187
10 NA 0.4860231 NA NA NA NA
11 0.7934172 NA 0.7934172 0.7934172 0.7934172 0.7934172
12 NA 0.4892819 NA NA NA NA
13 0.7694731 NA 0.7694731 0.7694731 0.7694731 0.7694731
14 NA 0.4648519 NA NA NA NA
15 NA NA NA 0.7734379 0.7734379 0.7734379
16 NA 0.4995469 NA NA NA NA
17 NA NA NA NA 0.7603955 0.7603955
19 NA NA NA NA 0.6073237 NA
ここで、列の見出しは Tube.numbers で、行のタイトルは Measurement.locations です。
x軸に測定位置、y軸に標準偏差(表の計算値)を使用して、各チューブ番号(6行)を1つのグラフにプロットしたいと思います
私はmatplotを試しましたが、同じ数の行と列を持たなければならないという必要性を回避できないようです。私はそれをループで行うことを好みます。
誰でも私を助けることができますか?
とても有難い!