Displaying 3 results from an estimated 3 matches for "meandata".
2011 Feb 03
0
Need advises on mixed-effect model ( a concrete example)
...Sylvain Cl?ment
University of Lille Nord de France
######################################
1st APPROACH Using Mean RT of participants as my DV.
######################################
# first fit a model with single fixed effects of condition and groups
result.lme1<-lme(RT~CONDITION+GROUP,data= meandata,random=~1|SUBJECT,
method="ML")
plot(effet.lme) # residuals are goodlooking
#then I try to add an interaction term in the model
result.lme2<-lme(RT~CONDITION*GROUP,data= meandata,random=~1|SUBJECT,
method="ML")
#I then compare both models
anova(result.lme1,result.lme2)
#g...
2009 Jul 01
1
Average of data files in a directory
...average when you have many data vectors in a directory. I have done this:
------------------------------------
setwd("/.../data/")
listfiles <- list.files(pattern=".pre") # list all datafiles in the
directory
loadfiles <- lapply(listfiles,read.csv) # read the csv files
meanData <- lapply(loadfiles,mean) # intended to compute the element-wise
average.
------------------------------------
Which does not produce the right result. Let me explain myself better... The
data files have the form (500 rows):
Time,Pressure
0,300
1,320
2,350
And I want to compute the average time...
2012 May 02
0
MCMCglmm priors including phylogeny
...= list(V = diag(2)*x, n=2),
G = list(G1=list(V = diag(2)*x, n=2)))
#model1
m1<-MCMCglmm(cbind(y1.logit,y2.logit) ~ trait:Ecotype - 1,
random=~us(trait):animal, rcov=~us(trait):units, family=c("gaussian","gaussian"),
prior=prior1, data=meandata, pedigree=meantree, nodes="TIPS",
thin=100, nitt=150000, burnin=30000, verbose=F)
#prior model2
prior2<-list(R = list(V = diag(2)*x, n=2))
#model2
m2<-MCMCglmm(cbind(reac.prob.logit,inert.prob.logit) ~ trait:Ecotype - 1,
rcov=~us(trait):units, family=c(&...