similar to: Plot confidence intervals

Hi all, I'm hoping I might be able to get some help with some issues specifying priors for MCMCglmm. I'm trying to fit a gaussian glmm using MCMCglmm to a data set with two (correlated) response variables. The response variables are both logit-transformed proportions (there are a few reasons why I've chosen these with gaussian error over binomal glmm, which I won't go into).

Hello, I noticed this while attempting to run a DOS software suite. The CALL command in wcmd doesn't seem to honor the set PATH. If I put this in my environment (via regedit): c:\windows;c:\windows\system;e:\;e:\test;f:\;f:\downloads\daphne\tools Then I put the following in f:\downloads\daphne\tools\bob.bat: path c:\windows;c:\windows\system;e:\;e:\test;f:\;f:\downloads\daphne\tools

I have two large files (2+MB), one of which is a corrupted version of the other. I tried to rsync the good file to the corrupted one, hoping it would repair the few bad bytes, but it did nothing. (I tried to be careful not to specify the --archive or --update flags, either explicitly or implicitly.) The original problem occurrend on a transfer from one computer to another. I then did a number

Hola He estado tratando varias formas de colocar el subíndice para los grados de libertad del F test y la expresión adj; y sencillamente no funciona. legend("topright",legend=c(expression('r[adj]=0.40'),expression('F[3,23]=15.43'),expression('P < 0.001')),col=black, bty='n',cex=0.9) Luego al tratar de colocar las cusriva utilizando par(font=3)

Hola Estoy tratando de correr un survival analysis usando un Cox regression model. Tengo una duda respecto a la organizacion del script. Tengo una variable que es -tamano del individuo- y quiero ver si hay diferencia en sobrevivencia respecto a tamano. Como diseno de campo los tamanos fueron ubicados de forma aleatoria en bloques al azar. Cuado planteo el script tengo algo como:

Hi, I am an R newbie and one thing I am having trouble with binding variables that I have created within one data frame into a new data frame when using na.omit(). To illustrate this problem I will give the example I am working on and the approah I have been using:- data.frame1<-filepath.... attach(data.frame1) #create a new variable using a function new.variable<-rep(1,length(weight3))

Dear R users, Recent changes to the MuMIn package now means that the model averaging command (model.avg) no longer returns confidence intervals, but instead returns zvalues and corresponding pvalues for fixed effects included in models. Previously I have used this package for model selection/averaging following Greuber et al (2011) where it suggests that one should use confidence intervals from

Hello there This is the first time I have used r-help message board so I hope I have got the right address. I am trying to check the residuals of a GLMM model(run using the package lme4). I have been able to check the residiuals of REMLs in lme4 using the following: m1<-lmer(vTotal~Week+fCollar+ (1|fCat), collars) res<-resid(m1) plot(res) qqnorm(res) library(MASS) par(mfrow=c(2,3))

I have been using lmer in lme4() to analyse the effect of dropping a term from the model as below anova.name<-anova(m1,m2) Where m1 is an original model and m2 has one term removed. I can then create my own type III tables for each variable in the model. However with many variables this becomes exceedingly time consuming! anova.lme in nlme() and Anova() in car dont seem to work for lme

I've been a NUT user for a while, on Linux, and I've been very pleased with the performance of it. I have a new machine here, that runs NetBSD, and I've been having fits with that. First, NetBSD has NUT 2.0.5 in the pkgsrc collection, and I don't really want to mess with that, since it's so far out of date. So, I built NUT 2.4.1, from source. All is well there, but I

I am running a p166 hardware with suse linux 7.0. I have downloaded the ssl version from your website. I did a ./config -d no-asm I did a ./make during the compilation of /apps the gcc -o openssl ... -lcrypto -lefence -ldl the following error message: /usr/i486-suse-linux/bin/ld: cannot find -lefence collect2: ld returned 1 exit status make[2]: *** [openssl] Error 1 make[2]: Leaving

Hi, how to estimate a the following model in R: y(t)=beta0+beta1*x1(t)+beta2*x2(t)+...+beta5*x5(t)+beta6*y(t-1)+beta7*y(t-2)+beta8*y(t-3) 1) using "lm" : dates &lt;- as.Date(data.df[,1]) selection&lt;-which(dates&gt;=as.Date("1986-1-1") &amp; dates&lt;=as.Date("2007-12-31")) dep &lt;- ts(data.df[selection,c("dep")]) indep.ret1

Hi I'm trying to run a R script in a computer grid using the Rscript interpreter, but the Rscript is not returning zero (even when the scripts processes succesfully) on its exit which causes the scheduler to detect an error and not records the output. Is there any way to get the Rscript returning zero ? -- "A critical section of code is like a bathroom. Only one person is allowed

Two questions about principal components analysis in R: Q.1) Hogenraad and McKenzie (1999) used Bruce Thompson's FACSTRAP program to bootstrap the factor loadings and scores in a principal components analysis. The input to the analysis was a word-word correlation matrix derived from a frequency count of x words across n texts. This is how they described their procedure: "Finally, we

Hello everyone: It's the first time I write to this mailing list. Sorry in advance if my doubt has already been posted before, but I have been checking the archives and I haven't been able to find anything satisfactory. I am running a mixed effects model with nested effects (site and pair, referred to barn swallow nests located in different places in different farms). My dependent

The xmh program provides an Xaw-based graphical user interface to the MH Message Handling System. This version continues the tradition of only releasing once every 5 years, with the latest set of build configuration & documentation improvements. Alan Coopersmith (4): config: Add missing AC_CONFIG_SRCDIR configure: Drop AM_MAINTAINER_MODE autogen.sh: Honor NOCONFIGURE=1

--4.DBC_6C_._.B.38 Content-Type: text/html; Content-Transfer-Encoding: quoted-printable <kettering aighfpzajx alpc tnscw vzdeoh q mp fdskwzwoil wwhkkj iotl ie laaam ci oavpmqkrgzmmcva rbu></annuli> <announceq></constructorw> L<sympathyq>owe</convenew>st Co<dene>st He</gliddenr>al<slosht>th Pla</bicepsy>ns wi<bloomingtonu>th

Dear all, I was wondering if anyone could help solve a problem of a missing interaction effect!! I carried out a 2 x 2 factorial experiment to see if eggs from 2 different locations (Origin = 1 or 2) had different hatching success under 2 different incubation schedules (Treat = 1 or 2). Six eggs were taken from 10 females (random = Female) at each location and split between the treatments,

Dear R users, I am using lmer combined with AIC model selection and averaging (in the MuMIn package) to try and assess how isotope values (which indicate diet) vary within a population of animals. I have multiple measures from individuals (variable 'Tattoo') and multiple individuals within social groups within 4 locations (A, B, C ,D) crucially I am interested if there are

I am currently having problems accesing our samba msdfs shares from Windows XP SP2, Windows 2000 and WinXP SP1 seem unaffected. I have attached the smb.conf and the exports files that define the shares. Running ethereal on the different clients seems to point to the fact that XP SP2 no longer seems to be using the netbios-ssn (139) port. for GET_DFS_REFERRAL instead its using microsoft-ds