数式または分析からのデータをシミュレートする一般的な方法はありますか?
実験計画データフレームからのデータのde novoシミュレーション。 Rに重点を置いて(ただし、他の言語のソリューションは素晴らしいでしょう)。 実験または調査の設計において、データをシミュレートし、このシミュレートされたデータを分析することにより、設計の長所と短所に関する素晴らしい洞察を得ることができます。 このようなアプローチは、統計的テストの理解と適切な使用にも不可欠です。 ただし、このプロセスはやや面倒な傾向があり、多くの場合、実験や調査でこの重要なステップをスキップするようになります。 統計モデルとテストには、データのシミュレーションに必要な情報の大部分が含まれています(分布の仮定または明示的な記述を含む)。 分析モデル(および関連する仮定(例:正常性とバランス)、因子のレベル、および有意性の尺度(p値など)が与えられた場合、シミュレートされたデータ(理想的には、 print()、predict()、simulate())。 このような一般化されたシミュレーションフレームワークは可能ですか? もしそうなら、そのようなフレームワークは現在利用可能ですか? 例、次のような関数が欲しい sim(aov(response~factor1+factor2*factor3), p.values=list(factor1=0.05, factor2=0.05, factor3=0.50, factor2:factor3=0.05), levels=list(factor1=1:10, factor2=c("A", "B", "C"), factor3=c("A", "B", "C"))) すなわち、次の一般化バージョン: sim.lm<-function(){ library(DoE.base) design<-fac.design(nlevels=c(10,3,3), factor.names=c("factor1", "factor2", "factor3"), replications=3, randomize=F) response<-with(design, as.numeric(factor1)+ as.numeric(factor2)+ as.numeric(factor3)+ as.numeric(factor2)*as.numeric(factor3)+ rnorm(length(factor1))) simulation<-data.frame(design, response)} または sim(glm(response~factor1+factor2*factor3, family=poisson), p.values=list(factor1=0.05, factor2=0.05, factor3=0.50, factor2:factor3=0.05), levels=list(factor1=1:10, factor2=c("A", …