Search Archives

Dear all, Below there are two, simple - I suppose, questions on using pmatch(): > pmatch("xx", c("cc","xxa")) [1] 2 > pmatch("a", c("cc","xxa")) [1] NA > pmatch("xx", c("cc","xxa","xxb")) [1] NA I would like that the...

dear all, for some reasons I would like to add some new components (vectors, matrices) to a "glmerMod" fit. I am not familiar with S4, so borrowing from S3 classes my first tentative was o <- glmer(y ~ ..) o at new...

Dear all, I'm interested in fitting time-series linear models with I(1) errors. Namely given y_t=a+b*t+u_t the random term u_t are such that u_t-u_{t-1}=e_t~iid...

dear all, Given a matrix A, say, I would like to apply a bivariate function to each combination of its colums. That is if myfun<-function(x,y)cor(x,y) #computes simple correlation of two vectors x and y...

A few days ago I uploaded to CRAN a new package called segmented. The package contains functions to fit (generalized) linear models with `segmented' (or `broken-line' or `piecewise linear') relationships between the response and one or more explanatory variables...

dear all, apologizes for this off-topic question. I am looking for a "ecological" dataset (n>100, say) including measurements of one or more growth variable and age. Could anyone to suggest the R package/URL where I can find...

Dear all, I writing my first package and everything seems to work (at least up to now) However when I try to build documentation (.dvi or .pdf), using Rcmd Rd2dvi.sh --pdf mypack.Rd I get a mypack.pdf whose...

Dear all, I would like to write a function like: myfun<-function(x,fn) {xx<-exp(x); x*fn(xx)} where fn is a symbolic description of any function with its argument to be specified. Therefore myfun(5,"2+0...

see ?array best, vito ----- Original Message ----- From: "Francisco do Nascimento Junior" <fnj at cin.ufpe.br> To: "R-help" <R-help at stat.math.ethz.ch> Sent: Wednesday, January 29, 2003 10:50 AM Subject: [R] Array of 3D Hi...

dear all, When I use all.vars(), I am interest in extracting only the variable names.. Here a simple example all.vars(as.formula(y~poly(x,k)+z)) returns [1] "y" "x" "k" "z" and I would like to...

Hi, The package nlme by Pinero & Bates provides both lme() and nlme(). The sintax for lme() is lme(y~x1+x2, random=~1+x1|z, data) where y is the response x1 and x2 two fixed effects and x1 the...

A few days ago I uploaded to CRAN a new package called segmented. The package contains functions to fit (generalized) linear models with `segmented' (or `broken-line' or `piecewise linear') relationships between the response and one or more explanatory variables...

dear all, This is a probably a silly question. If I type > grep("x",c("a.x" ,"b.x","a.xx"),value=TRUE) [1] "a.x" "b.x" "a.xx" Instead, I would like to obtain only "a.x" "b...

Hi all, I'm looking for some expert and/or volunteer which is intersted in discussing with me about Extended-Quasi-Likelihood methods and relevant implemantation in R. The aim is conjoint modelling of mean and dispersion. in GLM. It...

dear Szymon, it is a bug (in the new version), thanks. It depends on the flat underlying relationship you are trying to estimate with a small sample.. I will correct it as soon as possible. Meanwhile you can use o1...

Dear all, I've just downloaded the miniR.exe, miniR-1.bin,....,miniR-8.bin in order to update R to 1.6.2 (I cannot download the file of 19Mb !). The installation seems to have worked fine, but I...

dear all, for some reason I am intersted in updating a glm taking variables from its model argument, namely: > dati<-data.frame(y=runif(10),x=1:10) > obj<-glm(y~x,data=dati) > obj$model[,c("A","a:b...

dear all, Could anyone explain me the behaviour of median() within by()? (I am running R.2.7.0) thanks, vito > H<-cbind(rep(0:1,l=20),matrix(rnorm(20*2),20,2)) > by(H[,-1],H[,1],mean...

Zoran Loncarevic wrote: > > Is there a way to fit mixed proportional odds models in R? As far as I know, no. (anyway have a look to J Lindsey's packages, I don't know) However MIXOR and friends at http...

dear all, Is it possible to extract the response vector from a fitted arima object? For instance in glm it is allowed, by: obj.glm<-glm(y~x) obj.glm$y #the response vector In arima I can't find...