search for: dum1

...ntile','.png',sep='') bitmap(file=filename) qqnorm(x) qqline(x) grid() dev.off() if (par2 > 0) { filename <- paste(server,'_',status,'-','lagplot1','.png',sep='') bitmap(file=filename) dum <- cbind(lag(x,k=1),x) dum dum1 <- dum[2:length(x),] dum1 z <- as.data.frame(dum1) z title <- paste('Lag plot (k=1), lowess, and regression line for ',server,'[',status,']',sep=' ') plot(z,main=title) lines(lowess(z)) abline(lm(z)) dev.off() if (par2 > 1) {...

I am not sure if this question has been asked before - but is there a procedure in R (in lm or glm?) that is equivalent to ENTER and REMOVE regression commands in SPSS? Thanks a lot! -- Dimitri Liakhovitski Ninah.com Dimitri.Liakhovitski at ninah.com

...8 > paneldata27<-plm.data(paneldata27,index=c("Account","Day.of.Sample")) > nlevels(paneldata27$Account) [1] 1573 > nlevels(paneldata27$Day.of.Sample) [1] 831 That is there will be 1573 individual effects and 831 time effects. > paneldata27.fe<-plm(LnMidpeak~DUM1,data=paneldata27,model="within",effect="twoways") Error in model.matrix.default(~id - 1) : cannot allocate vector of length 1854301163 Is my panel data set just too large for R? Thanks very much in advance, Peter Steele-Mosey [[alternative HTML version deleted]]

Hi, I'm trying to create a linear model for a dataset that has a breakpoint e.g. # dummy dataset x <- 1:20 y <- c(1:10,seq(10.5,15,0.5)) plot(x,y) I've modelled this using the following formula: temp <- lm(y ~ x*(x<=10)+x*(x>10)) I want to be able to omit the intercept (i.e. force the line through zero) from the first of these segments (x<=10) so that I'm only