3

私は本当の初心者で、修士論文のためにいくつかの金属チューブの材料損失に関するデータを分析しようとしています. 異なるチューブの間隔での材料損失の標準偏差を比較したいと考えています。いくつかのサブマトリックスを作成しtapply、標準偏差を計算しました。

次のスクリプトがあります。

myfunctionSD <- function(mydata) { return(sd(mydata,na.rm=TRUE))}

Alltubes <- tapply(datIn$Material.loss.interval,
                   list(as.factor(datIn$Measurement.location),
                        as.factor(datIn$Tube.number)),
                   myfunctionSD)

その出力は、列ヘッダーTube.numberと行タイトルのテーブルMeasurement.locationです。

dput(head(Alltubes))

structure(c(0.871073958553372, NA, 0.697795091282526, NA, 0.838624866472886, 
NA, 0.726992791242471, 0.807567484588899, NA, 0.598675787394729, 
NA, 0.510990323891863, 0.81510216193526, NA, 1.09844645540173, 
NA, 0.839816448199645, NA, 0.63972340253115, NA, 1.11485875917537, 
NA, 0.883318358663128, NA, 0.86706340125676, NA, 1.21565055332783, 
NA, 1.24692213662875, NA, 0.704210691776757, NA, 0.962002980998362, 
NA, 1.00703215272093, NA), .Dim = c(6L, 6L), .Dimnames = list(
    c("1", "10", "11", "12", "13", "14"), c("1", "2", "3", "4", 
    "5", "6")))

明確にするために：1、2、3、4、5、6はチューブ番号と列のヘッダーであり、1、10、11 ...は測定場所と行のヘッダーです。

データのさまざまなプロットを作成したいと思います。そのうちの 1 つは、x 軸にチューブ番号、y 軸に導出された標準偏差の広がりを示す Boxplot です。

どこでも検索して、さまざまなことを試しましたが、グラフが表示されません。

誰かが私を助けてくれることを願っています。

4

1 に答える 1

0

na.rm=TRUE必要なことを行うために特別な関数を作成する必要はありません。次のように渡すことができますtapply。

Alltubes <- tapply(datIn$Material.loss.interval,
                   list(as.factor(datIn$Measurement.location),
                        as.factor(datIn$Tube.number)),
                   sd, na.rm=TRUE)

次に、次を使用できます。

par(mfrow=c(2,2), mar=c(4,4,1,1), oma=c(1,1,1,1))

boxplot(Alltubes, names=colnames(Alltubes), xlab="Tube numbers",
        ylab="standard deviation")

boxplot(t(Alltubes), names=rownames(Alltubes), 
        xlab="Measurement locations", ylab="standard deviation")

matplot(colnames(Alltubes), t(Alltubes), xlab="Tube numbers", 
        ylab="standard deviation", type="b", lty=1, pch=19)

matplot(rownames(Alltubes), Alltubes, xlab="Measurement locations", 
        ylab="standard deviation", type="b", lty=1, pch=19)

プロット

于 2013-12-28T05:18:46.207 に答える