similar to: \Sexpr in MikTex

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

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'd like to create a S4 class contains only one type of data. However, the number of slots varies. For example, I want to create a class "a", each slots in "a" contains numeric value only. > setClass("a", contains = "numeric") If I want to create an object "a" with only one slot containing c(3,4), I will write >

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

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 everyone, Can anyone give me some hint/suggestion for unzipping a file in R? The help file for "zip.file.extract" function seems to be difficult to understand. For example, the name of my zip file is "refGene.txt.gz", which is located in "C:\Documents and Settings\Desktop\test1". I would like to unzip this file to the same directory. Here's

Hi, I would like to know about which threads will be preempted by which on my OpenSolaris machine. Therefore, I ran a multithreaded program "myprogram" with 32 threads on my 24-core Solaris machine. I make sure that each thread of my program has same priority (priority zero), so that we can reduce priority inversions (saving preemptions -- system overhead). However, I ran the following

Hi everyone. I have a problem running a Windows appilcation with wine. The application is a Acucobol program. When I run "wine wrun32.exe myprogram" it works ok. When I run "wine wrun32.exe -C config.cfg myprogram" wine says myprogram doesn't exist. Any help? -- Marcelo Arroyo

I have made many different attempts to run a R program in Batch mode. I run a DOS command like: C:\Progra~1\R\rw1061\bin\rterm.exe -q --no-restore --no-save myprogram.R out.txt What happens: 1. DOS Shell starts up well 2. R starts 3. I get the messages: "ARGUMENT 'myprogram.R' __ignored__" "ARGUMENT 'out.txt' __ignored__" 4. R Shell is ready to

hi all, I have a small script that seems to be exiting and hitting the sleep 60... The smwebsocket just connects to the web socket provided and outputs the data. This works manually. the myprogram just opens a database and writes the line... My desire is to run the smwebsocket, connect to the websocket and output the data (line by line) when we get a line that matches Location take that line

Hello everybody I'm trying to write a simple version of matlab's "quiver". The idea is that I have fluid with velocity defined on a grid. I have a matrix of x-components of velocity and a matrix of y-components and I want to see the overall flow pattern. (I work with 2D fluid mechanics problems). My first-stab function is below: quiver <- function(u,v,scale=1) # first

dear all, I created a dll file from a fortran program. How I can recompile it in R? I use "dyn.load(F:/myprogram.dll)" without any error but when I run "source(F:/myprogram.r)" this error occurs: "Error in file(file, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(file, "r", encoding = encoding) :

Dear Wu Gong and Peter Ehlers, thank you very much for your help debugging my script. Now I have a general following up question: Is there a straightforward way to rearrange the following dataset so that all first letters of each column will be combined in one column, all the second letters in a second column, all the third ones in a third column and so on, resulting in 7 columns, i.e. for

Dear List, First off, my deepest gratitude to the Sweave developers: this tool has improved my quality greatly. A question in my work I use \Sexpr{} statements scalar values and the xtable package for all manner of tables. What I'd like to do is to use a vector inline, rather than a whole separate table. Something like: %%%%%%%%%%%%%%%% begin code % Latex junk % Sweave block:

I would like to create \examples{} in the manual dynamically, and while it is possible to do this with a \Sexpr at the top level, R CMD check issues a warning for it. (See below.) Is it intentional that \Sexpr is not allowed at the top level? The Rd grammar allows this, but R CMD check does not. Is there any other way to generate/modify the \examples{} section dynamically? Thanks, Gabor In

Hi all, Currently Sweave does not format the (especially numeric) value in \Sexpr{}, and I often have to round() the numbers explicitly, otherwise I will get more than 10 digits, which is not necessary for me in most cases. Is there a way to specify the number of digits to be kept for values in \Sexpr{}? e.g. can we make \Sexpr{} respect getOption('digits')? Thanks! Regards, Yihui --