search for: idh

Dear all, Just wondering if someone could help me out converting my code from a for() loop into a foreach() loop or using one of the apply() function. I have a very large dataset and so I'm hoping to make use of a parallel backend to speed up the processing time. I'm having trouble getting selecting three variables in the dataset to use in the foreach() loops. My for() loop code is:

Hi, I've been a bit confused by different wyas we specify random effects in lmer and MCMCglmm i just want to clear something. When I want to look for intersexual genetic correlations in the trait, is it equivalent to treat this trait for opposite sexes as separate traits and include the term idh(trait):animal - to treating this as a single trait and fitting idh(sex):animal? Do these two ways of specyfying random effects differ in this case? Cheers, sz.

Can anyone give me some advice, we have just set up a linux file server and we want to offer up a samba share of a directory, however we want it to be accessible by everyone even those who do not have an entry in the passwd file, can someone tell me what I need to do to acheive this ? many thanks in advance. Ian internal support engineer

...Specifically, the MCMCglmm course notes, available on the package's CRAN page, uses the following example on page 102: m5d.1 <- MCMCglmm(art ~ trait - 1 + at.level(trait, 1):fem + at.level(trait, 1):mar + at.level(trait, 1):kid5 + at.level(trait, 1):phd + at.level(trait, 1):ment, rcov = ~idh(trait):units, data = bioChemists, prior = prior.m5d.1, family = "zipoisson", verbose = FALSE) I have been unable to find an answer to the following questions and would appreciate any guidance: 1) Does at.level(trait, 1) index the poisson latent variable, or is that at.level(trait,2)?...

Thanks for your answers Stephen and Ben, I hope I am posting on the correct list now. I managed so far to run the multinomial model with random effect with the following command: MCMCglmm(fixed=cbind(Apsy,Mygl,Crle,Crru,Miag,empty) ~ habitat:trait,random=~idh(trait):mesh,family="multinomial12", data=dataA,rcov=~trait:units) (where multiple responses are different species, Habitat the explanatory variable and Mesh the random effect) The main question I am facing now is: Why the multinomial model fit does not provide by default a parameter es...

...bility+ELL*pctELL+Disability*pctDisability+ pctMale+factor(yr>=2008)*factor(Grade), random=~Schoolid+Teacherid+Studentid, family="categorical", data=data[data$Grade>=4,]) Error in MCMCglmm(FCATprofLevel ~ transferTypeCat + factor(yr >= 2008) + : please use idh(trait):units, us(trait):units or trait:units for error structures involving catgeorical data with more than 2 categories Does anyone know how I can fix it? Thanks -- View this message in context: http://r.789695.n4.nabble.com/How-to-write-random-effect-in-MCMCglmm-tp3567056p3567056.html Sent fr...

...632 6.467 3.09e-10 *** Mutant.IDC 98.7443 23.1632 4.263 2.55e-05 *** Mutant.IDD 97.2203 23.1632 4.197 3.37e-05 *** Mutant.IDE 118.9820 23.1632 5.137 4.49e-07 *** Mutant.IDF 241.8537 23.1632 10.441 < 2e-16 *** Mutant.IDG 107.4883 23.1632 4.640 4.80e-06 *** Mutant.IDH 105.7664 23.1632 4.566 6.74e-06 *** Mutant.IDI 517.4650 23.1632 22.340 < 2e-16 *** Mutant.IDJ 19.7777 23.1632 0.854 0.393735 Mutant.IDK 47.4240 23.1632 2.047 0.041313 * Mutant.IDL 3.2542 23.1632 0.140 0.888347 Mutant.IDM 180.9638 23.1632 7.813 5.63e-14 ***...

...try to fit a model without interactions I can't do it because I obtain the following error: > fit<-MCMCglmm(Abundancia~BA+CI, random=NULL, > family="zapoisson",data=Trucha) Error in MCMCglmm(Abundancia ~ BA + CI, random = NULL, family = "zapoisson", : please use idh(trait):units or us(trait):units or trait:units for error structures involving multinomial data with more than 2 categories or zero-infalted/altered/hurdle models I don't know where is the problem, maybe because my original data is organised as (obviously with much more data): Abundance BA...

Hi Kirill, just to notify you that your al175 driver is being removed from the NUT tree, as of 2.4.0-pre1. if you wish to see it entering the tree again, please contact the Development mailing list to talk about it. Merry Christmas and happy New Year. Arnaud -- Linux / Unix Expert R&D - Eaton - http://www.eaton.com/mgeops Network UPS Tools (NUT) Project Leader -

Hi list we have a client who complains that mails from a certain sender are always put into his spamfolder. Our strong guess is that this is an issue with the webmail (Roundcube), but I need some help in understanding the output from dovecot's mailbox-debug. In order to narrow down the source of the problem we activated the mailbox-debug from Dovecot. Now we can see that the mail is moved

Hi! I'm having a small problem with an assignment I have to do. We want to do an anova on some values, given for four different types of medicine. The different sample numbers are not identical (eg for molecule A we have 8 values, for B we have 14, etc.). What would the most elegant way of getting this info into R to do an lm be? I usually try to use data.frames but R categorically

Dear R people Could you please help. Basically, there are two variables in my data set. Each patient ('id') may have one or more diseases ('diagnosis'). It looks like id diagnosis 1 ah 2 ah 2 ihd 2 im 3 ah 3 stroke 4 ah 4 ihd 4 angina 5 ihd .............. Q: How to make three data sets: 1. Patients with ah and ihd 2. Patients with ah but no ihd 3. Patients with ihd but no

Hello list. I'm trying to rebuild the 2.6.9.67.0.20.EL kernel, but it fails even without modifications. How did I try it? Created a (non-root) build environment (not a mock ) Installed the kernel.scr.rpm and did a rpmbuild -ba --target=`uname -m` kernel-2.6.spec 2> prep-err.log | tee prep-out.log The build failed at the end: Processing files: kernel-xenU-devel-2.6.9-67.0.20.EL Checking

I wonder is someone can offer an explaination and solution to the following :- using two PC's one running SCO Unix with Samba and one using Win95 we get the following timings when copying files from the Samba machine 4Mb file 2 min 53 sec using Samba via router 46 sec using FTP via router 31 sec using Samba direct 35 sec using FTP direct 26Mb file 10 min 4 sec Samba via router 4

We have a perculiar problem and wonder if anyone can shed some light or advice ? We have a variety of machines running SCO Unix 3 .2v4.2 with samba 1.9.16, we also have an NT 4.0 server, these are all sharing up various sections of their drives. There is a Red Hat Linux running using smbmount to mount all these drives onto a directory tree, the root of this tree is shared up using samba so raher

Can anyone offer any help with improving the following timings :- Using 3 identical PC's one a win 95 client, one an NT 4.0 Server and one a Linux (Redhat 5.0, Kernel 2.0.32) with Samba (1.9.17p4) I get the following results when copying a directory tree up/down, the tree consists of 183 files in 38 directories totalling 8.5Mb. Win95 -> NT 33secs NT -> Win95

I have recently upgraded to Redhat 5.2 and Samba 2.0.3 and now our backup is only going so far (uses smbclient to tar to tape) and reporting no more room on device, I estimate it is hitting the problem at around 5-6Gb - as its a 25Gb tape and has worked fine on older versions of redhat/samba I think the tape is ok. Has any one any ideas or suggestions on what to try ? Many Thanks Ian

Can someone tell me where to find what error messages mean or enlighten me as to what smb_dont_catch_keepalive: server->data_ready == NULL means ? Many thanks Ian

Hi I started studying Bioinformatics at my engineering school (INSA de lyon!) in september. We started learning statistics and did two classes of R hands-on, but I want to learn to use R better. (As in for my every day calculations - I'm in front of a computer anyway so why should I have to take out my TI-30?) I'm having trouble finding some of the simpler stat functions though.