Search Archives

dear all, I do not if it is a nonsense question.. Is it possible in the R session to get the name of the current .Rdata file that I ran? I mean: suppose I double click the file myfile.Rdata...

dear all, Is the following intentional? Am I missing anything in documentation? d<-data.frame(y=rnorm(10,5,.5),exp=rnorm(10), age=rnorm(10)) formula(lm(exp(y)~exp+age, data=d)) #--> exp(y) ~ exp + age formula(lm...

Dear all, The help file for the generic function vcov states "Classes with methods for this function include: 'lm', 'glm', 'nls', 'lme', 'gls', 'coxph' and 'survreg' (the last two in package 'survival')." Since, I am not able to use vcov...

No, as far as I know. However the brlr package performs bias-reduced logistic regression that might be useful for you. (for instance yields ML estimates when there is separation in data and even CMLE do not exist) See the...

Dear all, is it possible to save warning messages in the object created by any function causing the warning message itself? For instance coxph() in survival package fits Cox regression model and returns the fitted model with possible warning messages...

Dear all, I found the following strange behavior of is.null() in the x component of a lm object. *However note that I'm running R-1.5.1, so probably the error (if someone) has been fixed in the...

Dear R-users, Probably this is rather general question and I apologize for this: I'm using R and other softwares on Win and I would like to install Linux on my PC in order to be able to install...

Dear all, In R<=1.6.2 I usually used the following code (in DOS prompt) to run R in batch mode C:\documents> Rcmd BATCH myfile.R and I could see the results (including warning messages) in the file...

> Hello: > > Is there a function that will do exact logistic regressions along the > line of LogXact? > No, as far as I know, there is no package to perform exact logistic regression (namely calculate the complete permutational distribution of the statistics...

I don't remember this test exactly; however testing homogenity of the odds ratio in K 2x2 tables is equivalent to testing interaction x:z in the logistic model: y~x+z+x:z where y is the (binary) response...

Dear all, Is there *any* reason explaining what I describe below? I have the following line myfun(x) If I type them directly in R (or copy/past), it works.. However if I type in R 2.6.1 > source...

Dear Thierry, I do not know fdrtool::monoreg(), however I think that the simple R^2 is meaningless in thise case. An unconstrained (non-monotonic) fit will have always a higher R^2 and more degrees-of-freedom. A better...

Dear all, In using the lme() function from the nlme package, I would like to specify a particular correlation structure for the random effects. For instance, for a 3 by 3 matrix, say, I am interested in assuming a `quasi...

dear all, hopefully it is of interest for someone.. I have written a set of functions to fit segmented mixed models, i.e. piecewise linear relationships with random changepoints in a likelihood/frequentist framework. The functions seem to work reasonably...

Dear all, I have the following (rather) strange problem.. For some reasons, I finally work with a variable whose name includes an R function, "a.log(z)", say. And that is a problem when I call it in a formula...

Hi all, I am interested in usign lmer/lmer2 to fit (G)LMM where the variance of the random effects is blockdiagonal or simply multiple of identity. Namely: G=blockdiag(s1*I_1,...,s_kI_k) G=s*I In...

dear s, If you are willing to use B-splines placing p knots, say, (rather than smoothing splines placing n knots), you could have a look to the functions splineDesign() or spline.des() (argument derivs) in the package "splines". vito...

dear all, I am trying to use an user-specified link function in glmer() (via the argument family), but the error is Error in famType(glmFit$family) : unknown link: ?my_link? my_link is a proper link function (of class...

See the xtable package. best, vito ----- Original Message ----- From: "Daniel E. Ho" <deho at fas.harvard.edu> To: <r-help at stat.math.ethz.ch> Sent: Tuesday, November 05, 2002 6:30 AM Subject: [R] LaTeX Output? > Hello all, > > Does...

dear all, I would like to get the lme call without fitting the relevant model. library(nlme) data(Orthodont) fm1 <- lme(distance ~ age, random=list(Subject=~age),data = Orthodont) To get fm1$call without fitting the model I use call...