Я пытался соответствовать одному из стандартных примеров из gRim, используя выбор назад или вперед или пошаговый выбор. Моя проблема заключается в попытке построить подходящую графическую модель. Окончательная графическая модель, кажется, не изменилась по отношению к исходной модели, с которой я начинаю. (прямые и основные эффекты или обратная / пошаговая и насыщенная модель). У кого-нибудь есть идеи о том, что я делаю не так? Это мой код:
if (!requireNamespace("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install(c("graph","Rgraphviz","gRbase","gRain","gRim"))
library(gRbase)
library(gRain)
library(gRim)
data(reinis)
msat <- dmod(~.^.,data=reinis)
mmain <- dmod(~.^1,data=reinis)
plot(as(msat, "graphNEL"))
#plot(msat) # doesn't work nowadays
mbaic <- backward(msat,k=2)
plot(as(mbaic,"graphNEL"))
mfaic <- forward(mmain,k=2)
plot(as(mfaic,"graphNEL"))
msaic <- stepwise(msat,k=2)
plot(as(msaic,"graphNEL"))