R のパッケージ lsmeans を使用して、治療因子のすべての観測値の平均推定値を取得しました (3 つのレベルしかないため、体系的な効果に含まれている実験計画のブロック因子のレベル全体で)。応答変数に sqrt 変換を使用しました。
したがって、Rで次のコマンドを使用しました。
最初の定義モデル
model<-sqrt(response)~treatment+block
次に、lsmeans を適用します。
model_lsmeans<-lsmeans(model,~treatment)
次に、これをプロットします
plot(model_lsmeans,ylab="treatment", xlab="response(with 95% CI)")
これにより、さまざまな処理の推定値と 95% 信頼区間を含む非常に優れたグラフが得られます。
問題は、このグラフが変換された応答用であるということだけです。
逆変換された応答 (二乗応答) でこの同じプロットを取得するにはどうすればよいですか?
新しいデータ フレームを作成して、lsmean、lower.CL、upper.CL を抽出しようとしました。
a<-summary(model_lsmeans)
New_dataframe<-as.data.frame(a[c("treatment","lsmean","lower.CL","upper.CL")])
そして、これらを二乗します
New_dataframe$lsmean<-New_dataframe$lsmean^2
New_dataframe$lower.CL<-New_dataframe$lower.CL^2
New_dataframe$upper.CL<-New_dataframe$upper.CL^2
New_dataframe
これにより、必要な見積もりと CI 境界の二乗が得られます。
問題は、上記の LS で行ったのと同じ推定値と CI のグラフを作成できないことです。
これどうやってするの？私が尋ねる理由は、私の記事に似たようなスタイルのグラフが欲しいからです。私はこの LSmeans プロットがとても気に入っていて、変換されていない応答変数を使用するのに非常に便利なので、すべてのグラフをこのスタイルにしたいと考えています。
ご助力ありがとうございます！すべてが明確であることを願っています！
敬具
ディトレフ