similar to: A problem with jointprior (Deal package)

Greetings all, I am trying to use Deal to learn a Bayesian Network for discrete data. I came across the following problem with jointprior function. >library(deal) > a <- read.csv("prepared.cluster1.csv") > fit <- network(a) > fit.prior <- jointprior(fit) Error in array(1, Dim) : 'dim' specifies too large an array I have tried debugging array and still

Greetings R-users, I have been using the fpc package in R to cluster my data. Speficically I am using kmeansruns clustering. I would like to know how I use R to partition data into clusters. What I am doing is as follows. # Use csv file as input ##################### wholeset = read.csv("Spellman800genesImputed.csv") # exclude first col (gene names) ##########################

Dear R-users, I have the following 2 files; A V1 V2 A 1 A 2 A 3 A 4 B 1 B 4 C 1 C 3 D 4 B V1 V2 process1 1 process2 2 process3 3 process4 4 I want to get this output C V1 V2 V3 A 1 process1 A 2 process2 A 3 process3 A 4 process4 B 1 process1 B

Hello all, I have a list of genes as follows (A) Number? GeneName 0???????????? YAR003W 1???????????? YAR007C 2???????????? YAR008W 3???????????? YBL035C 4???????????? YBL111C 5???????????? YBL112C And I have a list of gene interactions as follows (B) 0 -> 1 1 -> 5 3 -> 4 2 -> 3 I want to replace the values in (B) with the corresponding values in (A) ie, I want to achieve the

Dear all, I recently started using the deal package for learning Bayesian networks. When using the jointprior function on a particular dataset, I get the following message: >tor.prior<-jointprior(tor.nw) Error in array(1, Dim) : 'dim' specifies too large an array What is the problem? How can I resolve it? Thanks, Tim Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm

"CALL FOR PAPERS" THE ISLAMIC SOCIETY IN SOUTH EAST ASIA CONFERENCE ON STATISTICAL AND MATHEMATICAL SCIENCES BANDUNG, 25 - 26 APRIL 2003 Fakultas Matematika dan Ilmu Pengetahuan Alam Universitas Islam Bandung Jalan Tamansari No. 1, Bandung 40116 Telepon +62 22 4203368 Ext. 136 or 135 Facsimile +62 22 4263895 E-mail: - fmipa_unisba at yahoo.com - wanda at

Greetings R-users, I have the following data called mydata in a data.frame Col1 Col2 Col3 Col4 Col5 1 2 4 6 7 8 8 7 3 5 4 4 5 6 7 I want to replace the data according to the following conditions Condition 1 if data <= 3, replace with -1 Condition 2 if data >=6, replace with 1 Condition 3 if data = 4 or data =5, replace with 0

-----Original Message----- From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On Behalf Of r-help-request at r-project.org Sent: Monday, February 11, 2008 12:00 PM To: r-help at r-project.org Subject: R-help Digest, Vol 60, Issue 11 Send R-help mailing list submissions to r-help at r-project.org To subscribe or unsubscribe via the World Wide Web, visit

Dear R-users I am looking for a help in using the "deal" package. I followed the manual and paper from the author's web site to learn it, as shown below, but I could not figure out how to access the local and posterior probability of the nodes in the constructed network. library(deal) data(ksl) ksl.nw <- network(ksl) ksl.prior <- jointprior(ksl.nw) banlist <-

Hi Si Jie, Again, please send questions to the list, not me. Okay, I may have worked out what you are doing. The program runs and produces what I would expect in the rightmost columns of the result "g". You are storing the number of each test for which the p value is less than 0.05. It looks to me as though the objects storing the results should be vectors as you are only storing 100 p

Hi Jeem, First, please send questions like this to the help list, not me. I assume that you are in a similar position to sjtan who has been sending almost exactly the same questions. The problem is not in the loops (which look rather familiar to me) but in your initial assignments at the top. For instance: scale parameter=(1,1.5,2,2.5,3) produces an error which has nothing to do with the

Hi all, I'm having problems with the installation of CentOS 5 for i386 to my Dell Percision T4500 Workstation. The installation run just fine. The installation DVD had been tested before installation. After installation, the installer ask me to reboot. after reboot, when the msg for LVM saying detection of my LVMs, 1 or 2 lines after that my monitor turned off, but i believed the cpu is still

hey friends, I am trying to secure my mail server.I have enabled TLS support in postfix(version postfix-2.1.5), now I am trying to configure squirrelmail(version 1.4.4-1 rpm) for tls/ssl support.In config.php i have choosen use_imap_tls=true and use_smpt_tls=true. Moreover If I send any mail from squirrelmail there are no entries for ssl or tls in maillogs whereas If I send the mail through

Dear all, > data <- data.frame(id=seq(1:10),x=runif(10)) > data id x 1 1 0.3604464 2 2 0.4813987 3 3 0.0160058 4 4 0.7165909 5 5 0.6092248 6 6 0.2413049 7 7 0.7981568 8 8 0.6093960 9 9 0.2887064 10 10 0.3485780 > selected.id <- sample(data$id,3,replace=F) > selected.id [1] 9 7 1 I want to select data with corresponding selected.id, namely >

Dear gluster experts, We have encountered a severe performance issue related to quota feature of gluster. My underlying fs is lvm with xfs format. The problem is if quota is enabled the io performance is about 26MB/s but with quota disabled the io performance is 216MB/s. Any one known what's the problem? BTW I have reproduce it several times and it is related to quota indeed. Here's the

Hi all, I am running two node RHEL3U8 cluster of below cluster version on HP servers connected via scsi channel to HP Storage (SAN) for oracle database server. Kernel & Cluster Version Kernel-2.4.21-47.EL #1 SMP redhat-config-cluster-1.0.7-1-noarch clumanager-1.2.26.1-1-x86_64 Suddenly my active node got rebooted after analysed the logs it is throwing below errors on syslog.I want

Dear All, Can you help me, with the code below how do I obtain the fisher information from it. Is my q<-replicate(1000,x) the right way to do simulation. thank you. x<-rweibull(100,0.8,1.5) q<-replicate(1000,x) z<-function(p){ beta<-p[1] eta<-p[2] log1<-(n*log(beta)-n*beta*log(eta)+(beta-1)*sum(log(x))-sum((x/eta)^beta)) return(-log1) } zz<-optim(c(0.5,0.5),z) zz Chris

*Scholarship and Job* are posted at* **http://cambodiajobs.blogspot.com*<http://cambodiajobs.blogspot.com/2007/12/aunseed-net-scholarship-for-master.html> ** *AUN/SEED-Net Scholarship for Master* (ASEAN University Network / Southeast Asia Engineering Education Development Network) The AUN/SEED-Net Project with support from Japan International Cooperation Agency (JICA) provides scholarships

To the R users community: (We apologize for multiples copies) (Please distribute) --------------------------------------------------------------------------- 1st Call for Papers - 2nd International Symposium on Distributed Computing and Artificial Intelligence (DCAI'09) June 10th - 12th , 2009 - Salamanca - Spain http://dcai.usal.es/ The International Symposium on Distributed Computing and

Hi. Thank you very much in advance for your help. I have generated data from two simple linear models and used k-means clustering (stats4) to identify two clusters in the generated data. Next, I would like to do simple linear regression for each separate cluster. I can do this if I first use the cluster labels to define two separate data frames with the subset function. However, I would