私は R の初心者であり、現在 R のベイジアン ネットワークを使用して結果を分析するコーディングを行っています。このコードには、確率とその 95% の資格間隔として結果が得られます。
これは、喫煙状況と社会階級を考慮して現在の喘息を予測するために使用しているコードの例です。私が欲しいのは、2 つのクエリから得た結果をテーブルに整理する方法です。約 48 のクエリがあり、それらすべての結果をコピーして貼り付けるのは面倒で時間がかかります。何か助けはありますか？
## probability of current asthma among non smoker in social class 1##
> cpquery(fitted.simpler_mice_hc, event = (c_asthma == 1) , evidence = (smoking_status == 0 & sei_class==1), n=100000)
[1] 0.09212951
> run_1<- replicate(100, cpquery(fitted.simpler_mice_hc, event = (c_asthma == 1) , evidence = (smoking_status == 0 & sei_class==1), n=100000))
>
> # to get the sampling quantile interval
> quantile(run_1, c(0.05,0.975))
5% 97.5%
0.08218354 0.09254333
> ### non smokers in sei2, age , gender,diabtets, dust and hypertension###
> cpquery(fitted.simpler_mice_hc, event = (c_asthma == 1) , evidence = (smoking_status == 0 & sei_class==2), n=100000)
[1] 0.08623836
>
> run_2 <- replicate(100, cpquery(fitted.simpler_mice_hc, event = (c_asthma == 1) , evidence = (smoking_status == 0 & sei_class==2), n=100000))
>
> # to get the sampling quantile interval
> quantile(run_2, c(0.05,0.975))
5% 97.5%
0.08279396 0.09144179
最初のクエリでは、0.09212951 が確率で、0.08218354 0.09254333 が間隔の下限と上限です。これらをクエリ用にテーブルまたは Excel で処理しやすいものに整理するには、それぞれをコピーして貼り付けるのではなく、どのようにすればよいでしょうか?