Search Archives

Dear all, I want to predict a new observation with nlme using fixed and random effects estimates from a previous fit. How can I do this? (My problem is different from the Loblolly example) Thanks a lot! Peter

Dear list, is there a function available which provides an estimate of the hazard function based on a cox proportional hazard model? I only found the cumulative hazard and the survival function as survfit options. Thanks for your help Peter

Dear list, I am looking for a way to fix the scale parameter when fitting a generalized linear model with gamma errors and log link. Is there something like "SCALE" such as in GLIM? As always thanks a lot. Peter

Dear all, I would like to simulate data from a nlme model using fixed effects and the variances of the random effects distribution. I scanned the help files and the mailing lists with no succes. Thanks for your help! Peter

Dear all, I am trying to use xyplot inside a function plotme<-function(dataframe) { xyplot(x~y|z,data=dataframe) } x,y,z are members of the data frame. When calling function library(lattice) plotme(dataframe) nothing happens. (R-1...

Dear all, within the help file for nonlinear regression an upcoming option for optional weights is mentioned. Is this already available? Thanks a lot! peter -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ...

Dear all, does any of you know a source for R code for the generation of correlated binary random numbers? Thanks a lot Peter Schlattmann -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ...

Dear all, is there a way to extract the individual variance-covariance matrix of the parameter estimates of the i-th individual using nlsList ? Thanks a lot Peter -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik...

Dear all, I have a problem with the readline editor of R since it doesn't work. Setting ./configure --with-readline does not the fix the problem. Does anyone have an idea how to fix this? Many thanks in advance...

Dear all, one perhaps stuoid question: If I want to update from 1.4.1 to 1.5.0 under Linux do I not to rebuild the whole software from scratch or is there an update mechanism as in "update...

Dear all, I tried to create a groupedData object, where the grouping factor is not ordered. Here ist the code: library(nlme) test<-groupedData(conc~Time|Subject,order.groups=F,data=as.data.frame(Theoph)) > getGroups(test) Levels: 6 < 7...

Dear all, I want to include covariates into the self-starting linear regression functions SSbiexp and SSfol. Is there a way to do this or do I need to define my own functions depending on the covariates? My second question...

Dear all, I need to stack the columns of a matrix with a variable number of columns. I am aware of new<-c(a[,1],a[,2]) if I have a matrix 'a' with e.g. two columns. For a...

Dear all, I try to perform weighted linear least squares wtih using glns. The weights are something like ww<-rep(c((0.1,0.9),n)) my dataframe is named test and I add ww to test by test$ww...

Dear all, I am stuck with the following problem: Using the Indometh data set I want to plot the predicted values of a simple bixeponential model together with the observed concentrations. I managed to plot the predicted values: datafi<-Indometh...

Dear all, I am trying to gain some experience with the function gnls from the nlme package. I tried to model the Theophyline data by trying to model the presumed dependency of the clearance on the body weight. This is...

Dear all, when performing a meta analysis I have two results obtained with logLik and lme, which I do not quite understand. The results are based on these data: study or var 1 0.10436 0.299111 2 -0.03046...