similar to: Censored outcomes - repeated measures and mediators

Hi, I have a dataset in which I would like to select rows based on matching conditions and return the maximum value of a variable else return one row if duplicate counts exist. My dataset looks like this: PGID PTID Year Visit Count 6755 53121 2009 1 0 6755 53121 2009 2 0 6755 53121 2009 3 0 6755 53122 2008 1 0 6755 53122 2008 2 0 6755 53122 2008 3 1 6755 53122 2009 1 0 6755 53122 2009 2 1 6755

Hello all, A really noddy question for you all: I''m trying without success to do some subhypothesis testing. Using simple anova model, with a toy dataset from a book. I have four factors A,B,C,D, and wish to test mu_C = mu_D. This is what I have tried: > contrasts(infants$group,how.many=1) <- c(0,0,1,-1) > contrasts(infants$group) [,1] A 0 B 0 C 1

Hello! I have problem with mediation analysis. I can do it with function mediate, when I have one mediator. But how I can do it if I have one independent variable and one dependent variable but 4 mediators? I have try function mediations, but it dosen't work. If I use mediate 4 times, each for every mediator, is it same? I want to know what is total mediate effect for 4 mediators. t.Mete

I simply put, > NEVER=subset(infants$bwt,ISNO1) > UNTILPREGNANT=subset(infants$bwt, ISNO2) > ONCENOTNOW=subset(infants$bwt, ISNO3) and I wanna combine those three. I do it like ISNO=NEVER&UNTILPREGNANT&ONCENOTNOW and R tells me 1: In NEVER & UNTILPREGNANT : longer object length is not a multiple of shorter object length 2: In NEVER & UNTILPREGNANT & ONCENOTNOW

Hi I have three steps like this: Given ''a user called "$name" of type $user_type and status $user_state in organisation "$organisation_name"'' do |name, user_type, user_state, organisation_name| $runs ||= 1 puts "MOO #{$runs}" $runs += 1 # ... end Given ''a user called "$name" with username

Hi all, Two questions: 1. Is there a way to evaluate models from lmer() with a poisson distribution? I get the following error message: library(lme4) lmer(tot.fruit~infl.treat+def.treat+(1|initial.size),family=poisson)->model plot(fitted(model),resid(model)) Error: 'resid' is not implemented yet Are there any other options? 2. Why doesn't the function estimable() in gmodels

Hello, I have been trying for weeks to compile a 64-bit Rblas. I started with ATLAS where I have had success in the past, but 64 bit was not behaving, and as each compilation takes between 9 and 12 hours, "test, check, and revise" was not going to be really viable. I therefore switched to OpenBLas (OPBL). I was successful in compiling R-3.0.1 and an OPBL-based BLAS for Windows 64bit

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

Many readers on this list probably don't know that S Poetry had an old home. S Poetry is a book in pdf about programming in the S language. It was written mainly with S-PLUS 3.x in mind when R was still an infant relative to what it is now. Some parts of the book are still quite accurate for R. In other places R is rather different from what is in the book. The FAQ on the differences

Ok, after familiarizing myself with clone it appears to me this is a bug in compiler-rt. >From the clone man page: In Linux 2.4 and earlier, clone() does not take arguments ptid, tls, and ctid. The source file passes those arguments without any fencing to check the Linux version. Also, ptid, tls, and ctid are only used in conjunction with certain flags (e.g., CLONE_PARENT_SETTID), but none

>From k.j.mcconway@open.ac.uk Wed Mar 28 16:46:49 2001 Hardly crucial, but I've come upon a potential error in the documentation of the 'swiss' datafram in the R base package. The description accurately matches what is said in the Mosteller and Tukey source quoted, but according to the data archived at Princeton (links from http://opr.princeton.edu/archive/eufert/switz.html), the

Welcome to "Ghetto" America! …...Yes, you heard me; I said "Ghetto" America! What is "Ghetto" America, I'm glad you asked. Before I go on please be warned, this explanation may not sit well with many of you, however it is no less true! "Ghetto" America is a society of individuals comprised primarily of African Americans who represent the

+Sergey Matveev On Wed, Jul 10, 2013 at 2:00 AM, Andy Jost <Andrew.Jost at synopsys.com> wrote: > Ok, after familiarizing myself with clone it appears to me this is a bug > in compiler-rt.**** > > ** ** > > From the clone man page:**** > > ** ** > > In Linux 2.4 and earlier, clone() does not take arguments ptid, tls, and > ctid. **** > > ** ** >

fixed @ R-forge. New version should appear on CRAN soon. Thanks for the report! Torsten > > ---------- Forwarded message ---------- > Date: Sat, 3 Sep 2011 23:56:35 +0200 > From: Ulrich Halekoh <Ulrich.Halekoh at agrsci.dk> > To: "r-help at r-project.org" <r-help at r-project.org> > Subject: [R] glht (multcomp): NA's for confidence intervals using

Hej, Calculation of confidence intervals for means based on a model fitted with lmer using the package multcomp - yields results for calpha=adjusted_calpha - NA's for calpha=univariate_calpha Example: library(lme4) library(multcomp) ### Generate data set.seed(8) d<-expand.grid(treat=1:2,block=1:3) e<-rnorm(3) names(e)<-1:3 d$y<-rnorm(nrow(d)) + e[d$block]

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:

No, you're quite right, I was wrong. I must have been very confused when I put up the first 'bug' report, #886. I noticed my error a couple of days after and put up a comment to that effect on R-bugs, but it seems to have vanished (or alternatively I did that wrong as well). I agree with you on the criterion for the choice of districts as well. There remains a minuscule discrepancy

At times we want to convert a two-sided formula to a one-sided formula. In S we can do this by dropping the second entry in the formula. In R that object no longer has a formula class. R> ttt <- score ~ age | Infant R> class(ttt) [1] "formula" R> length(ttt) [1] 3 R> ttt[-2] [[1]] ~ [[2]] age | Infant R> class(ttt[-2]) NULL R> do.call("~",

We are trying to read a sas export file into R. (Works fine on Mac by ftping the export file from solaris box and then importing to R with read.xport, this gives rational numbers etc.) Trying to do this on an IBM Power p655 running linux, reading the same sas export file as used on MAC, in this case the results are not rational. Many of the numbers are interpreted as INF etc. Is it possible

R 2.15.1 OS X and Windows Colleagues, I have a file that looks that this: TABLE NO. 1 PTID TIME AMT FORM PERIOD IPRED CWRES EVID CP PRED RES WRES 2.0010E+03 3.9375E-01 5.0000E+03 2.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 1.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 2.0010E+03 8.9583E-01