Displaying 20 results from an estimated 120 matches similar to: "incorrect handling of NAs by na.action with lmList (package nlme) (PR#13658)"
2012 May 04
0
oddsratio and some basic help on epitools
Here is a working snippet.
library(epitools)
mat <- matrix(c(10,15,60,25,98, 12,10,70,28,14, 9,11,68,10,12
,8,13,20,11,58) ,ncol=2)
colnames(mat) <- c("treatmentA","treatmentB")
row.names(mat) <- paste("Cond",rep(1:10,1))
dimnames(mat) <- list("Condition" = row.names(mat), "instrument" =
colnames(mat))
> mat
instrument
2012 May 04
0
oddsratio epitool and chi-square
Here is a working snippet.
library(epitools)
mat <- matrix(c(10,15,60,25,98, 12,10,70,28,14, 9,11,68,10,12
,8,13,20,11,58) ,ncol=2)
colnames(mat) <- c("treatmentA","treatmentB")
row.names(mat) <- paste("Cond",rep(1:10,1))
dimnames(mat) <- list("Condition" = row.names(mat), "instrument" =
colnames(mat))
> mat
instrument
2012 May 04
0
epitools question
Here is a working snippet.
library(epitools)
mat <- matrix(c(10,15,60,25,98, 12,10,70,28,14, 9,11,68,10,12
,8,13,20,11,58) ,ncol=2)
colnames(mat) <- c("treatmentA","treatmentB")
row.names(mat) <- paste("Cond",rep(1:10,1))
dimnames(mat) <- list("Condition" = row.names(mat), "instrument" =
colnames(mat))
> mat
instrument
2009 Oct 26
1
explalinig the output of my linear model analysis
Hi,
I am new in statistics and i manage to make the linear model analysis but i
have some difficulties in explaining the results. Can someone help me
explalinig the output of my linear model analysis ? My data are with 2
variables habitat (e,s) and treatment (a,c,p) with multiple trials within.
Thank you in advance
Call:
lm(formula = a$wild ~ a$habitat/a$treatment/a$trial)
Residuals:
Min
2006 Jun 09
0
interaction terms in regression analysis
G'day,
My problem is I'm not sure how to extract effect sizes from a nonlinear
regression model with a significant interaction term.
My data sets are multiple measurements of force response to an agonist
with two superimposed treatments each having two levels.
This is very similar to the Ludbrook example in Venables and Ripley.
The experiment is that a muscle is exposed to an agonist
2009 Jul 30
3
What is the best method to produce means by categorical factors?
I am attempting to replicate some of my experience from SAS in R and assume
there are best methods for using a combination of summary(), subset, and
which() to produce a subset of mean values by categorical or ordinal
factors.
within sas I would write
proc means mean data=dataset;
class factor1 factor2
var variable1 variable2;
RUN;
producing an output with means for each variable by factor
2005 Jan 24
4
lme and varFunc()
Dear R users,
I am currently analyzing a dataset using lme(). The model I use has the
following structure:
model<-lme(response~Covariate+TreatmentA+TreatmentB,random=~1|Block/Plot,method="ML")
When I plot the residuals against the fitted values, I see a clear
positive trend (meaning that the variance increases with the mean).
I tried to solve this issue using weights=varPower(),
2003 Sep 08
0
lmList with NAs
Hello R-Helpers,
I was trying to use the lmList function to get the lmList graphic similar to Pinheiro and Bates (pg 33). I did not have a problem creating the graphic when I used the Orthodont data frame or 2 other data sets when there are no missing values.
My data has missing values. Do I need to remove the missing values before the lmList function will work?
for a small example:
> a
2006 Aug 18
0
lmList and missing values
Dear all,
I have a question on handling of missing values in lmList. My data set
have continuous predictor and response, x and y, and a grouping variable
group.id. All these variables have NAs and the data set also has several
other variables that also contains NAs.
To create the lmList-object seems to work fine:
y.list <- lmList(y ~ x | group.id, data=mydata, na.action=na.omit)
However,
2008 Oct 06
1
question on lmList
Hi list,
Using the lmList function from ?nlme package? I get the following error message:
?Error in !unlist(lapply(sum.lst, is.null)) : invalid argument type?
# this is the syntax used
# fm_cream <-lmList(MULTDV~TIME|cod, data=zz, na.action=na.omit)
# summary(fm_cream)
When I generate a simple example the function works, but on the real data I get that message. I imagine that it is a
1998 Oct 21
0
anovalist.lm or anova.lmlist?
In R, we currently have the functions
anovalist.lm
and anova.glmlist
S / S-plus has
anova.lmlist
anova.glmlist
On the other hand, the [n]lme package (library) of Doug Bates and Jose
Pinheiro has an "lmList" class and an anova.lmList(.) method for that.
We are considering to use
anovalist.lm
and anovalist.glm
instead of the S/S-plus names mentioned above.
These functions
2015 Sep 20
0
summary(lmList) issues
Over on StackOverflow someone is trying to get this code to work:
library("nlme")
fm2Pixel.lis<-lmList(pixel~day+I(day^2)|Dog, Pixel)
summary(fm2Pixel.lis)
(results: Error in `[<-`(`*tmp*`, use, use, ii, value = lst[[ii]]) :
subscript out of bounds)
http://stackoverflow.com/questions/32680306/plotting-individual-confidence-intervals-for-the-coefficients-in-the-lmlist-fit
2018 Jul 28
1
possible bug in plot.intervals.lmList
Dear R-devel members,
I think I've found a minor bug in plot.intervals.lmList. ( The guide
https://www.r-project.org/bugs.html suggests to report it here, as I
do not have a bugzilla account.)
Here is a minimal reproducible example to demonstrate the problem:
fm1 <- lmList(distance ~ age | Subject, Orthodont)
plot(intervals(fm1),ylab="a")
This results in: "Error in
2003 Feb 13
1
sorting in lmList object
Hi all,
Forgive me if this is an obvious one....
I want to make a plot of confidence intervals from an lmList object with a
collection of simple linear models (lm(y~x)) using:
plot(intervals(mylmList))
and sort the plot by increasing mean values for the intercept. Is there an easy
way of doing this? I've tried the "order()" and "sort.list()" functions, but I
suspect
2004 May 01
1
changes to y-axis labels in lmList intervals plot
Dear List,
I am plotting lmList objects using plot(intervals()) in nlme
package. I want to make changes to the y-axis labels. When I
try to change cex of y-axis labels using the following:
fm1 <- lmList(distance ~ age | Subject, Orthodont)
plot(intervals(fm1), scales=list(y=list(cex = .7)))
I receive: Error in bwplot(formula = group ~ intervals | what,
data = structure(list( : formal
2004 Sep 13
1
lmList for glm
Greetings all,
Is there a way to do the equivalent of a lmList (fit separate models for each subject) in the package nlme for a glm?
Thank you very much,
Peter B.
[[alternative HTML version deleted]]
2004 Oct 11
0
lmList - strange problem
I am using lmList, which takes 'Data' partitioned according to the
levels of a grouping factor (gf) and individual 'lm' fits are obtained
for each 'data'partition, using a model defined as in "lm". So my call
to lmList looks something like:
> mg.lis <- lmList(strmbf ~ age1c + gender1 + sysbp.clinic +
+ diabp.clinic + ldl1 + mets.total + bmi1c |
2005 May 13
1
error in plot.lmList
Hello,
in R-2.1.0 I'm trying to prodice trellis plots from an lmList object as described in the help for plot.lmList. I can generate the plots from the help, but on my own data plotting fails with an error message that I cannot interpret (please see below). Any hints are greatly appreciapted.
kind regards,
Arne
> dim(d)
[1] 575 4
> d[1:3,]
Level_of_Expression SSPos1 SSPos19
2005 Oct 04
1
lmList error message
It seems a problem sending this out. I am resending
it.
Hello Rs,
After running lmList, I've gotten an error message
"Error in lm.fit(x, y, offset = offset, singular.ok =
singular.ok, ...) : NA/NaN/Inf in foreign function
call (arg 4)" I'm not able to understand what it
means. Could anyone help me with it?
Thanks all in advance.
Kevin
### commands and output ###
2007 Feb 11
2
Suppresing default text in pairs.lmList() in package = nlme
I would like to suppress the text 'Scatter Plot Matrix' that appears
under the plot. Could someone please suggest how?
_____________________________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400 Charlottesville, VA 22904-4400
Parcels: Room 102 Gilmer Hall
McCormick Road Charlottesville, VA 22903
Office: B011