Displaying 11 results from an estimated 11 matches for "fitlms".
Did you mean:
films
2006 Nov 09
1
Extracting the full coefficient matrix from a gls summary?
Hi,
I am trying to extract the coefficients matrix from a gls summary.
Contrary to the lm function, the command fit$coefficients returns
only the estimates of the model, not the whole matrix including the
std errors, the t and the p values.
example:
ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group <-
2006 Jan 11
3
how to obtain "par(ask=TRUE)" with trellis-plots
Dear alltogether,
how can a delay like possible with par(ask=TRUE) be attained while using
trellis-plots within a loop or something like that?
the following draws each plot without waiting for a signal
(mouse-klick), so par() does not work for that:
library(nlme)
for(i in 1:3)
{
fitlme <- lme(Orthodont)
par(ask=TRUE) # does not work with trellis....
print(
2012 Nov 07
5
Calling R object from R function
Hi,
Can you please help me with this please?
What I am trying to do is call a vector from R function and used in the new
function
So I create 4 functions with these arguments
M11 <- function(TrainData,TestData,mdat,nsam) {
ls <- list()
I have few statments one of them is
vectx <- c(,1,2,3,4,5,6,6)
vectz <- c(12,34,5,6,78,9,90)
and then................
ls(vectx=vtecx,vectz=vectz)
2007 Jan 12
2
Magnitude of trend in time series
Hello,
I am analyzing some climate time series data using the Mann Kendall package
and was wondering if there was a way to calculate the trend using Sen's
nonparametric estimator slope in R?
Thank you in advance,
Barry
_________________________
Barry Baker, Ph.D.
Global Climate Change Initiative
The Nature Conservancy
2424 Spruce St., Suite 100
Boulder, CO 80302
Tel: (303)-541-0322
Fax:
2010 Mar 15
3
the problem about sample size
Hi all:
I am a user of "JM" package.
Here's the problem of "sample size".
The warning is:
Error in jointModel(fitLME, fitSURV_death, timeVar = "time", method = "piecewise-PH-GH") :
sample sizes in the longitudinal and event processes differ.
According to the suggestion of "missing data",I use the same data set(data_JM) without any
2012 Oct 05
0
jointModel error messages
I contacted the package developer and that lead to me removing events at time
0 (or subjects with only 1 longitudinal measurement). I then still had the
error message "Can't fit a Cox model with 0 failures" which I have managed
to avoid by adding 1.8*10^(-15) to all my survival times, any number greater
than this also works but nothing smaller! Any explanation of this would
help!
2005 Dec 12
2
convergence error (lme) which depends on the version of nlme (?)
Dear list members,
the following hlm was constructed:
hlm <- groupedData(laut ~ design | grpzugeh, data = imp.not.I)
the grouped data object is located at and can be downloaded:
www.anicca-vijja.de/lg/hlm_example.Rdata
The following works:
library(nlme)
summary( fitlme <- lme(hlm) )
with output:
...
AIC BIC logLik
425.3768 465.6087 -197.6884
Random effects:
2012 Nov 16
2
R-Square in WLS
Hi,
I am fitting a weighted least square regression and trying to compute
SSE,SST and SSReg but I am not getting SST = SSReg + SSE and I dont know
what I am coding wrong. Can you help please?
xnam <-colnames(X) # colnames Design Matrix
fmla1 <- as.formula(paste("Y ~",paste(xnam, collapse=
2011 Nov 22
0
plotting output from LME with natural cubic spline
I have used LME to fit a mixed effects model on my data. The data has
274 subjects with 1 to 6 observations per subject. Time is not linearly
associated with the outcome, so I used ns to fit a natural cubic spline
with 3 auto knots. Subject and the natural cubic time of spline are both
treated as random effects. This model has run without any problem, but
now I would like to plot trajectories for
2007 Sep 12
0
constructing an lm() formula in a function
I'm working on some functions for generalized canonical discriminant
analysis in conjunction with the heplots package. I've written a
candisc.mlm function that takes an mlm object and computes a
candisc object containing canonical scores, coeficients, etc.
But I'm stumped on how to construct a mlm for the canonical scores,
in a function using the *same* right-hand-side of the model
2003 Jan 29
1
Scoping rule problem -- solved
Thanks to some comments from Brian D. Ripley, I found my error:
I should not have given a data argument to lm() after creating a
formula-object. This obviously confused things...
Thanks again, I've really learnt again a bit more on R-programming...
Cheers, Winfried
---------------------------------------------------------------------
E-Mail: Winfried Theis <theis at