search for: cornwallis

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))

Can you provide the summary(m2) results? > -----Original Message----- > From: Simon Pickett [mailto:S.Pickett at exeter.ac.uk] > Sent: Wednesday, August 16, 2006 7:14 AM > To: Doran, Harold > Cc: r-help at stat.math.ethz.ch > Subject: [SPAM] - RE: [R] REML with random slopes and random > intercepts giving strange results - Bayesian Filter detected spam > > Hi again,

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

Hi, This patch adds R600/SI disassembly text to compiled object files, when a code dump is requested, to assist debugging in Mesa clients. Here's an example of the output in a Mesa client with a corresponding patch and RADEON_DUMP_SHADERS set: Shader Disassembly: S_WQM_B64 EXEC, EXEC ; BEFE0A7E S_MOV_B32 M0, SGPR6 ; BEFC0306

On Wed, Oct 09, 2013 at 08:06:42PM -0500, Jay Cornwall wrote: > Hi, > > This patch adds R600/SI disassembly text to compiled object files, when > a code dump is requested, to assist debugging in Mesa clients. > > Here's an example of the output in a Mesa client with a corresponding > patch and RADEON_DUMP_SHADERS set: > > Shader Disassembly: > >

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi, I have the following system set up: Server: Ubuntu 7.04, Samba 3.0.24, http://www.esuna.co.uk/~jay/misc/smb.conf Client: Vista Ultimate 64-bit, NetBIOS over TCP/IP disabled, LMHOSTS lookup disabled. Share mapped to drive Z:. Up until a few days ago this was working fine. Read/write transfer rates of over 20MB/s onto a RAID-1 array, speedy

Hi all, I have been bumbling around with r for years now and still havent come up with a solution for plotting reliable graphs of relationships from a linear regression. Here is an example illustrating my problem 1.I do a linear regression as follows summary(lm(n.day13~n.day1+ffemale.yell+fmale.yell+fmale.chroma,data=surv)) which gives some nice sig. results Coefficients:

Hi, I've encountered a problem invoking LLVM-generated x86_64 code on Windows from MSVC-compiled code. There seems to be disagreement between the two compilers over the calling convention for passing a 12 or 16 byte structure by value. e.g. struct aggregate { int a; int b; int c; int d; }; int test(aggregate byVal) { ... } LLVM chooses to pass byVal with a and b packed into rcx

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

Hi, I would like to (if possible) set the default width and height for graphs at the start of each session and have each new graphic device overwrite the previous one. I only know how to do this using windows(width=,height=...) which opens up a new plotting device every time, so I end up with lots of graphs all over the place until I get the one I want! Thanks in advance, Simon Simon Pickett

Hi, I am an R beginner and having problems structuring my REML models. I have a model with y=weight x1=time x2=timesquared id=individual identity I need to structure the model such that in the random effects there is a constant intercept for all individuals but a separate individual slope for both x1 and x2 (a coefficient score for every individual).

Hi everyone, I have been using REML to derive intercepts and coeficients for each individual in a growth study. So the code is m2 <- lmer(change.wt ~ newwt+(newwt|id), data = grow) Calling coef(model.lmer) gives a matrix with this information which is what I want. However, as a test I looked at each individual on its own and used a simple linear regression to obtain the same information, then

Hi R experts, I have been using ReML as follows... model<-lmer(late.growth~mtf+year+treat+hatch.day+hatch.day:year+hatch.day:treat+ mtf:treat+ treat:year+ year:treat:mtf+(1|fybrood), data = A) then I wanted to plot the results of the three way interaction using lines() as follows... tmp<-as.vector(fixef(model)) graph1<-plot(mtf,fitted(f2), xlab=list("Brood Size"),

Hi All, I would like to separate the media traffic from the signalling. Can Asterisk send and receive media (rtp) traffic from a secondary network interface? Thanks, Harel ________________________________ This electronic message and any files transmitted with it are confidential and intended solely for the use of the individual or entity to whom they are addressed. If you are not the named

Hi, I have been dabbling with str() to extract values from outputs such as lmer etc and have found it very helpful sometimes. but only seem to manage to extract the values when the output is one simple table, any more complicated and I'm stumped :-( take this example of the extracted coeficients from a lmer analysis... using str(coef(lmer(resp3~b$age+b$size+b$pcfat+(1|sex), data=b)))

Hi all, I have been using princomp() recently, its very useful indeed, but I have a question about how to specify the rows of data you want it to choose. I have a set of variables relating to bird characteristics and I have been using princomp to produce PC scores from these. However since I have multiple duplicate entries per individual (each bird had a varying number of chicks), I only want

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).

...20 GB disk space and 64 mb memory - PCMCIA adaptor. ** I would like to start a conversation via e-mail with someone on your team to get me back on track. ** Can you help. Steven Steven L. Sprague Large Enterprise Direct Advisory Project Manager (919) 486-9879; T/L 526-9879 IBM Corporation 3039 Cornwallis Road Bldg 629 Office Q101 RTP, NC 27709 Internet: ssprague@us.ibm.com