similar to: abline() or predict.lm() when log="x"

Hello all, Is AIC calculated incorrectly in lmer? It appears as though it uses AIC = -2*logLik - 2*#parms, instead of -2*LogLik + 2*#parms? Below is output from one of many models I have tried: Generalized linear mixed model fit using PQL Formula: cswa ~ pcov.ess1k + (1 | year) Data: ptct50.5 Family: poisson(log link) AIC BIC logLik deviance 224.8466 219.19 -114.4233 228.8466

Hello, I have a basic question. Sorry if it is so evident.... I have the following data file : http://ekumen.homelinux.net/mydata.txt I need to model Y~X-1 (simple linear regression through the origin) with these data : load(file="mydata.txt") X=k[,1] Y=k[,2] aa=lm(Y~X-1) dev.new() plot(X,Y,log="xy") abline(aa,untf=T) abline(b=0.0235, a=0,col="red",untf=T)

if I want to specify y-coordinates for the heights of horizontal lines to go across a plot. x <- c(1,2,3,6,4,8,4,7) y <- c(3,2,7,4,5,4,5,6) h <- c(3,5,7) plot(x,y) abline(y=h) However I got error message: Warning message: parameter "y" couldn't be set in high-level plot() function (I tried abline(h=y) , it;s not what I want also) Do u know why?

Dear whom it may concerns, I am Jia Xu, a summer intern analyst at Citi Research. I am also a graduate student studying Financial Engineering at Cornell. I am running regression analysis with Rstudio now. I have been experiencing difficulty with the abline function. No matter I call it directly in linear regression or call it through residual plots, or plot results from lasso regression, it

Dear R-users, I am using R(a package igraph) to calculate certain topological features of networks. When I try to draw a plot between these features I get an error. Following is the code I am using : *> plot(met_eco_deg,met_eco_bet) > lmout<-lm(met_eco_bet ~ met_eco_deg) > abline(lmout) Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : plot.new

Full_Name: Stephen Eglen Version: 1.4 OS: Redhat Linux 7.1 Submission from: (NULL) (128.252.204.36) abline() produces a spurious line in addition to the correct line, at least with the X11 driver. The postscript file generated also has nan values in it, which causes an error under ghostscript. xs <- c(0, 150, 300) ys <- c(75, 40, 23) par(mfrow=c(1,2)) plot(xs,ys, log="")

I am trying to plot the results of a random effects model where each subject is allowed to have his, or her, own intercept. I use xyplot to plot the data, lme to compute the regression and then try to put a summary regression line on the xyplot. As can be seen by the output pasted below, I am getting an error message, Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : (list)

Hello, This is an (hopefully) improved question of one I posted several weeks ago. Does anyone know of a function for fitting "two-part" models? These models are designed to handle count data with so many zeroes that they can't be fit well with zero-inflated Poisson models or other 'typical' GLMs. My understanding is that they work by first fitting a binomial model to

Hello, I am not having success in a simple task. Using the effects package, I would like to add reference lines at probability values of 0.1 – 0.6 on a plot of the effects. The plot command works, but following up with an abline command produces the message “plot .new has not been called yet”, and of course the reference lines were not added. Looking through past R help lists, there was a

I really want to do this: abline( a=tan(-kT*pi/180), b=kY-tan(-kT*pi/180)*kX ) where kX,kY and kT are vectors of equal length. But I can't do that with abline unless I use a loop, and I haven't figured out the least unelegant way of writing the loop yet. So is there a way to do this without a loop? Or if I am to resort to the loop, what's the best way of doing it considering that I

Hi folks, Win7 64bit R 1.12.0 I run following command on R:- > ToothGrowth > attach(ToothGrowth) > plot(dose,len) > matrics=lm(len~dose) > abline(metrics) Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : plot.new has not been called yet Only a grey diagram is displayed without content > plot(abline(metrics)) Error in int_abline(a = a, b = b, h = h, v

Hi All, I am running a scatter plot and trying to add a best fit line. I use an abline function, but get no line drawn over the points. I also get no error. I arm using V 2.10.0 on Windows 7. Here is my code, including the SAS transport file import: require (foreign) require (chron) require (Hmisc) require (lattice) clin <- sasxport.get("y:\\temp\\subset.xpt") attach(clin)

Greetings, I'd be grateful if a good Samaritan helps me to approach this problem.... with my data, I've created the following model lm(formula = OUTCOME ~ VAR1 + VAR2) summary(model) Call: lm(formula = OUTCOME ~ VAR1 + VAR2) Residuals: Min 1Q Median 3Q Max -1.4341 -0.3621 0.1879 0.4994 0.7696 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.89020

Hello, if I do: x <- c(0.5,1,3,6,10,20,40) y <- 10-log(x)+rnorm(7,0,0.05) r1 <- lm(y ~ log(x)) plot(log(x),y) abline(r1) # I get a nice plot with the regression line almost over the points. but: plot(x,y,log="x") abline(r1) gives me exactly the same plot for the points but the regression line is completely off ! I would like the plot with the real values of X on

Hello and Happy New Year to all R-Users !!! I would like to plot a lattice graph with a logarthmic y axis and add two reference lines that is : ref<-c(0.0070, 0.0096) graph1<-xyplot(data$y1 ~ as.numeric(strptime(data$x1, format="%Y-%m-%d %H:%M:%S")) | as.character(data$Code), list(y = list(log = T)), xlab="X'", ylab="Y", panel = function(...) {

I am trying to make R-1.8.1 on (SunOS shell1 5.8 Generic_108528-15 sun4u sparc SUNW,UltraAX-i2). I did ./configure make Configure output seems ok. The make proceeds until the following line appears, repeated indefinitely (until I break): ./config.status: ./confstat28489-19881/subs.frag: cannot overwrite existing file I suspect that this may involve write permissions (and maybe the umask set in

Hi, I want to have a horizontal line on a mosaic plot with "vcd" package. This would give me an idea where is 0.5 proportion in a cell. Using "mosaicplot" function of "graphics" package, I can draw a line using "abline." But, with "mosaic" function of "vcd" package, I have tried to use "abline" function, which complains

Thanks in advance, I have searched and read the archives extensively and have been unable to find a question (solution) exactly mimicking my problem. I have created a scatterplot simulation of data through which I have drawn a linear regression in this manner: points(y~x, pch="*", col="black") lm(y~x) fm=lm(y~x) abline(fm, col="red") This works. The problem

I would like to add two horizontal lines representing acceptible drug levels to a trellis plot. I tried using abline and I get an error that "plot.new has not been called." See below. xyplot(FK~WEEK|Event1/MRN, data=FKdat.o1) abline(h=5) abline(h=10) Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : plot.new has not been called yet Any help that can be provided

Dear R-help, I have a four-part question about regression, matrices, and sandwich package. 1) In the sandwich package, I would like to better understand the meat() function. >From the bread() documentation, for a simple OLS regression, bread() returns (1/n * X'X)^(-1) That is, for a simple regression (per the documentation on bread()): MyLM <- lm(y ~ x) bread(MyLM)