Displaying 4 results from an estimated 4 matches for "meancol".
2010 Aug 05
3
eval-parse and lme in a loop
...having trouble getting an eval-parse combination to work with lme in a for loop.
Consider the code below. The call to lme outside the loop works. The call to aov inside the loop works, but the call to to lme inside the loop does not.
The latter fails with
Error in model.frame.default(formula = ~meanCol + Group + c1 + c3 + Subject, :
variable lengths differ (found for 'Group')
and I don't know why the eval-parse portion of the code is node being replaced as it is in the call to aov. Anybody got any idea what's wrong?
rm(list=ls())
library('nlme')
df=data.frame(
Gr...
2006 Mar 28
1
variables as arguments in formulae in aov?
...--------------
roi.errs=read.table("../Group.results/allGroupsROI.errs", header=T, sep="\t");
roi.stops=read.table("../Group.results/allGroupsROI.stops", header=T, sep="\t");
# extract the names of the columns containing the means
roi.names=names(roi.errs);
meanColumns=roi.names[4:length(roi.names)];
## open the output file, split=T also spits it to the screen
sink("../Group.results/ROI.errs.txt",split=T)
cat("********************************************************************************\n");
cat("*** Summary for ROI.errs\n");...
2006 Apr 12
0
how to interpret the results of a simint call
...to interpret the results of the call to
simint. Could someone help me with this?
I've attached the code and the results it produced. There are 3 groups
(ctrl,long,short) and I'm trying to see if the data in the Mean_12
column (in the example below) differs based on group membership.
for (meanCol in meanColumns) {
roi.err.aov<-aov(eval(parse(text=meanCol)) ~ Group, data=roi.errs);
HSD.err<-simint(eval(parse(text=meanCol)) ~ Group, data=roi.errs,
type="Tukey");
cat("################################################################################\n");
cat(&quo...
2004 Sep 02
3
confidence intervals
Dear R users;
Im working with lme and Id like to have an idea of how
can I get CI for the predictions made with the model.
Im not a stats guy but, if Im not wrong, the CIs
should be different if Im predicting a new data point
or a new group. Ive been searching through the web and
in help-lists with no luck. I know this topic had been
asked before but without replies. Can anyone give an
idea of