Dear all,

I learn something from my colleagues about the beauty of Bayesian
network. I want to use the deal package to analyze my genetic data. I
have 3575 subjects. For each subject, I have 20 genetic marker
information (discrete, three levels) and his or her plasma cholesterol
levels (continuous). I want to construct the markov blanket of plasma
cholesterol levels. I tried the deal in my computer and computer
cluster. Due to so many genetic markers, I can not successfully pass the
calculation of the prior probability of my network. The maximum
variables I could put into the network is about 12. I can not add more.

The algorithms (gs, iamb, hc) in bnlearn package can not deal the
discete and continuous variables at the same time. Could someone give me
some suggestions how I should do it? 

Thanks in advanve,

Kevin


	[[alternative HTML version deleted]]

