similar to: adding horizontal lines to a trellis plot

Greetings, \Sexpr{} has worked on MikTex with earlier versions of R, I can not seem to get it to work with R-2.4.0. I run Sweave with MikTex using the following statement in Run under Accessories. latex -include-directory="C:\Program Files\R\R-2.4.0\share\texmf" "C:\Documents and Settings\Suzette\Desktop\MyFile\MyProgram.tex" Any advice would be much appreciated, Suzette

Greetings, I implemented BRugs to run the EWOC model with a cohort size N=1. I output the simulation data using bugsdata(data), where data is the following list. > data $Dose [1] 140 $Y [1] 0 bugsdata(data) puts out the file data.txt as follows. list(Dose=1.40000E+02, Y=0.00000E+00) after I type the following line modelData("data.txt")

I am wondering how to estimate the survival curve for a particular case(s) given a coxph model using this example code: #fit a cox proportional hazards model and plot the #predicted survival curve fit <- coxph( Surv(futime,fustat)~resid.ds+strata(rx)+ecog.ps+age,data=ovarian[1:23,]) z <- survfit(fit,newdata=ovarian[24:26,],individual=F) zs <- z$surv zt <-

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

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

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)

I have a data vector x. When I try qqline(x) I get the following error: Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : plot.new has not been called yet And a blank plot appears. Can anybody help? What am I doing wrong? Thanks, Scotty _________________________________________________________________ Contest [[alternative HTML version deleted]]

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

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

I'm working over an ssh connection without X11 graphics. I'm making a plot, the first stage of drawing which takes a long time. I want to experiment with adding details. Here is what I was hoping to do, which results in error. ## Draw the master plot on png dev 2 png(file="master.png") plot(1:10) ## Save a copy on png dev 3 png(file="copy1.png") dev.set(2)

Greetings, Does anyone know how to get an id number in the little header above each individual plot within a trellis plot? The default seems to be to print the word id and add a line indicating on a linear scale where the current id sits. Thanks in advance for any help you can send, Suzette ================================= Suzette Blanchard, Ph.D. Research Scientist Frontier Science

I am not able to download packages ("xtable") from CRAN not even just a .zip file, I get that the page can not be found. Please could you help? Suzette ================================= Suzette Blanchard, Ph.D. Research Scientist Frontier Science

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

Greetings, I would like to use a cubic spline or smoother to model the fixed effects within nlme. So far the only smoother I have been able to get to run successfully in nlme is smooth(). I tried smooth.spline: fixed=list(lKa~1,lCL~smooth.spline(BSA, df=3)) the error I got was the following. Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : invalid

Greetings, I have a simulation of a nonlinear model that is failing. But it does not fail til way into the simulation. I would like to look at the run that is failing and maybe I could if I could capture the seed for the failing run. The help file on set.seed says you can do it but when I tried rs<-.Random.seed print(paste("rs",rs,sep=" ")) I got 626 of them so I

Greetings, I used latex type code in my Rd files. The pdf version using R CMD Rd2dvi --output=PKtools.pdf --pdf --title="PKtools" PKtools man came out quite nice. However, my current HTML version does not tex the latex so there is latex code in the files and looks bad. example problem code: AIC$_c$ or \begin{itemize} \item NLME: \begin{itemize} \item population level:

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)

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="")

Hello, I've some questions concerning log-transformations and plotting of the regression lines. So far as I know is it a problem to log-transform values smaller than 1 (0-1). In my statistics lecture I was told to do a log(x+1) transformation in such cases. So I provide here a small example to explain my questions: # Some example data for testing a1