Search Archives

See the Thomas Lumley home page There is code at http://faculty.washington.edu/tlumley/weave.html best, vito ----- Original Message ----- From: "Matthias Fischer" <Matthias.Fischer at wiso.uni-erlangen.de> To: <r-help at stat.math.ethz.ch> Sent...

Are you looking for Re() and friends? Toy examples: > abs(Re(3+4i)) [1] 3 > abs(Re(-3+4i)) [1] 3 > abs(Re(3)) [1] 3 see ?complex for further details on complex numbers in R vito ----- Original Message ----- From...

dear lucas, yes you are right, segmented does not handle 'lars' objects. Out of curisity, are you interested in selecting the number of breakpoints or in selecting additional covariates with linear parameters? vito Il 06/06/2012 0.01, Lucas...

Dear Aur??lie, I think that for *fixed* (i.e. assumed known) amount of smoothing, you can use a simple LRT by comparing the two candidate models. BTW, have a look to the mgcv or gam packages for a general...

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...

Dear all, It appears that the function is.factor() returns different results when used inside the apply() function: that is, is.factor() fails to recognize a factor.. Where is the trick? many thanks, vito > df1<-data.frame(y=1:10...

Dear Mario, I don't know whether the $\ell$ symbol is available .. However you can use the LaTeX psfrag/pdfrag packages to convert tags in latex symbols.. hope this helps, vito Mario dos Reis wrote: > Is it possible to use...

dear John, The package segmented can help you. ?plot.segmented includes a Poisson example best, vito "Sorkin, John" <jsorkin at som.umaryland.edu> ha scritto: > Is there an R package that will perform a piecewise continuous > Poisson regression? I want...

As arima.sim() simulates stationary ARMA "errors" if your underlying model is additive I think you can type, for instance, just: x<-1:100 #time variable mu<-10+.5*x #linear trend y<-arima.sim(length(x), model=list(ar...

As far as I know, the Hosmer-Lemeshow test is not a good choice, as it depends on the groups to be formed to compute the statistic. The design and/or hmisch packages by F.Harrell should include some alternative...

Hi all, I am experiencing the following quite strange (at least in my knowledge) problem in a simple function like the following: fn<-function(y,v=2){ n<-length(y) y<-y[(v+1):(n-v)] plot(y,type="l...

Hi there, I'm interested in estimating likelihood- (or simply deviance- for GLM) based CI. I use the following code, but is there a more efficient way to do it? obj<-glm(y~x+z....) #fit the full model beta...

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...

There are at least two R packages dealing with changepoint estimation, segmented and strucchange. Two possible relevant papers are available: 1)Journal of Statistical Software for strucchange (2002, Vol.7, Issue2) 2)Rnews for segmented (2008, 8/1: 20-25...

Dear *all*, (how many.. seeing that this appears to be the first message? :-)) My question concerns LMM and it is not strictly related to developments of the lme4 package. In a LMM framework is it possible to set a heterogeneity...

Dear all, Two of the arguments of my function fn(), say, are a glm-obj and a variable in the data.frame where there are the variable of the obj too. Is there a way to build a function that...

dear Christofer, Try the following d<-expand.grid(a=1:3,b=1:4) d$freq<-rpois(12,5) o<-glm(freq~factor(a)+factor(b)+I(a*b), family=poisson, data=d) vito Il 02/04/2012 9.34...

This does not answer exactly to your question, anyway..: If you are planning to use latex, the package movie15 allows to include media files in your document (to be processed via pdflatex) vito Bruno C. wrote: > With version 8 of...

hi feng, a possible solution is b1<-apply(a,2,list) and possibly lapply(b1,unlist) if you want exactly the output equal to "list(a[, 1], a[, 2])" best, vito Il 18/01/2011 13.53, Feng Li ha...

Dear all, Some days ago I posted the same message below, but unfortunately with no reply. So, apologizes for this my re-sending, but I hope there is now someone on-line that can help me.... In using the lme...