similar to: removing part of a string

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(exp(y)~., data=d)) #--> exp(y) ~ age variable 'exp' (maybe indicating "experience") is not included in the model. The same happens with 'log' (and

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 obtain "y" "x" "z" Where is the trick? many thanks vito -- ==================================== Vito M.R. Muggeo Dip.to Sc

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(): my.cc<-call("lme.formula", fixed= distance ~ age, random = list(Subject = ~age)) However the two calls are not the same (apart from the data

dear all, apologizes for this OT =========================== 28th International Workshop on Statistical Modelling (IWSM), Palermo (Italy) 8-12 July 2013. http://iwsm2013.unipa.it Dear friend, For your information, I would like to bring to your attention that deadline for submission of abstracts is FEBRUARY 18 If you are still interested to visit Palermo (and taste its specialities :-))

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. ls() returns the names of the objects in the current workspace (that is saved in myfile.Rdata). In the current R session, I would like to obtain "myfile.Rdata". Is it possible? Thanks in

dear all, apologizes for this off topic. I would like to inform you that registration and paper submission for the 28th International Workshop on Statistical Modelling (IWSM) to be held in Palermo (Italy) 8-12 July 2013 is now open at http://iwsm2013.unipa.it Register at http://iwsm2013.unipa.it/?cmd=registration and then submit your abstract. Deadlines for Abstract submission is February 4,

dear all, It appears that glmnet(), when "selecting" the covariates entering the model, skips from K covariates, say, to K+2 or K+3. Thus 2 or 3 variables are "added" at the same time and it is not possible to obtain a ranking of the covariates according to their importance in the model. On the other hand lars() "adds" the covariates one at a time. My question

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.x" How is it possible to get this result with grep()? many thanks for your attention, best, vito --

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) INDICES: 0 V1 V2 -0.2101069 0.2954377 --------------------------------------------------------------------------------------------------------------------- INDICES: 1 V1

Hello, Try this. ss1 <- "z:f(5, a=3, b=4, c='1:4', d=2)" ss2 <- "f(5, a=3, b=4, c=\"1:4\", d=2)*z" fun <- function(s) sub("(\\().*(\\))", "\\1\\2", s) fun(ss1) #[1] "z:f()" fun(ss2) #[1] "f()*z" Hope this helps, Rui Barradas On 5/21/2018 2:33 PM, Vito M. R. Muggeo wrote: > dear all, > I am stuck

I would use sub("\\(.*\\)", "()", s) It is essentially the same as Rui's suggestion, but I find the purpose to be more clear. It might also be a little more efficient. HTH Ulrik On Mon, 21 May 2018, 15:38 Rui Barradas, <ruipbarradas at sapo.pt> wrote: > Hello, > > Try this. > > > ss1 <- "z:f(5, a=3, b=4, c='1:4', d=2)" >

Dear all, It appears that MASS::polr() and Design::lrm() return the same point estimates but different st.errs when fitting proportional odds models, grade<-c(4,4,2,4,3,2,3,1,3,3,2,2,3,3,2,4,2,4,5,2,1,4,1,2,5,3,4,2,2,1) score<-c(525,533,545,582,581,576,572,609,559,543,576,525,574,582,574,471,595, 557,557,584,599,517,649,584,463,591,488,563,553,549) library(MASS) library(Design)

Dear all, This is off-topic, however I hope someone can give me useful suggestion.. Given the regression model y = b0 + b1*x + e I am interested in testing for positive coeffs, namely H0: b0>0 AND b1>0 H1: b0,b1 unconstrained It is simple to estimate the model under H0 and H1 (there are several suggestions on the Rlist about estimation but nothing about testing..) perform a likelihood

Dear list, I would like to know: 1. After I have used the R code (http://pscl.stanford.edu/zeroinfl.r) to fit a zero-inflated negative binomial model, what criteria I should follow to compare and select the best model (models with different predictors)? 2. How can I compare the model I get from question 1 (zero-inflated negative binomial) to other models like glm family models or a logistic

Dear all, I have tried to consider a correlation structure in lme (package lme4), but without success. I have used something like: > risul<-lme(y~x+ z , data=mydata, random=~ x | g, correlation = corAR1()) but the result is the same as: > risul<-lme(y~x+ z , data=mydata, random=~ x | g). Can anybody help me? Antonella ************************************************** Prof.

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("code.R") ##code.R includes the above line Error in inherits(x, "data.frame") : object "d" not found namely myfun() does not work correctly. In particular the

This drop-in patch increases the performance of the get_checksum1() function on x86-64. On the target slow CPU performance of the function increased by nearly 50% in the x86-64 default SSE2 mode, and by nearly 100% if the compiler was told to enable SSSE3 support. The increase was over 200% on the fastest CPU tested in SSSE3 mode. Transfer time improvement with large files existing on both ends

I've read up some more on the subject, and it seems the proper way to do this with GCC is g++ and target attributes. I've refactored the patch that way, and it indeed uses SSSE3 automatically on supporting CPUs, regardless of the build host, so this should be ideal both for home builders and distros. Getting the code to build right in c++ mode (checksum_sse2.cpp only) was a bit of an

Dear all, I am dealing with the following (apparently simple problem): For some reasons I am interested in passing variables from a dataframe to a specific environment, and in fitting a standard glm: dati<-data.frame(y=rnorm(10),x1=runif(10),x2=runif(10)) KK<-new.env() for(i in 1:ncol(dati)) assign(names(dati[i]),dati[[i]],envir=KK) #Now the following two lines work correctly:

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, for instance: > myname<-"a.log(z)" > dd<-data.frame("a.log(z)"=1:10,y=rnorm(10)) > o<-lm(y~1,data=dd) >