similar to: Metafiiles into Word R 1.9.0

Last week (Wed 9/10/2003, "regression questions") I posted a question regarding the use of gnls( ) and its dissimilarity to the syntax that nls( ) will accept. No one replied, so I partly answered my own question by constructing indicator variables for use in gnls( ). The code I used to construct the indicators is at the end of this email. I do have a nagging, unanswered

Dear R-helpers, I would like to change the default value of a functions that is call by a package. The function deparse() is called by a set of functions from a library and because I don't think is appropriate to change the library, I am trying to change the default value in deparse(). At the moment the only solution I have come across is to place a function in my working environment like

Dear R-helpers, Would anybody have a compiled binary version (for Windows) of the Dopt library that is in the development site?. best regards, IOsu -------------------------------------------------------------- Jes?s Mar?a Fr?as Celayeta Dublin Institute of Technology School of Food Sci. and Env. Health. Cathal Bruha St., Dublin 1. Ireland Phone: +353 1 4024459 Fax: +353 1 4024495

Dear R-helpers, I need some help to produce a set of contour plots that I am trying to make in order to compare surfaces between the levels of a factor. For example: library(lattice) g <- expand.grid(x = 60:100, y = 1:25, ti = c("a","b","c")) g$z <-with(g, (-1e-4*x-1e-3*y-1e-5*x*y)*(ti=="a") + (1e-2*x-1e-3*y-1e-4*x*y)*(ti=="b") +

I got a question from a fellow PhD student that work with spectrum analysis in Excel and now he has lots of spectrums that needs to be smoothed, which would be nice to be able to do in batch. Is there an R package that can do: Savitzky-Golay smoothing for reflectance spectral data or a function that does something similar. _______________________________________ Henrik Andersson

Dear R-helpers, I’ve got a summary of results from a by() call that I am making with a list of more than two of factors not very different from the example in the by() help page require(stats) by(warpbreaks[, 1], warpbreaks[, -1], summary) The result of the command gives a list of the form wool: A tension: L Min. 1st Qu. Median Mean 3rd Qu. Max. 25.00 26.00

Hi there, I'm trying to fit a growth curve to some data and need to use a weighted least squares estimator to account for heteroscedasticity in the data. A weights argument is available in nls that would appear to be appropriate for this purpose, but it is listed as 'not yet implemented'. Is there another package which could implement this procedure? Regards, Robert Brown

Dear R-helpers, I am trying to reshape a data frame that is in a wide format into a long format one. I have come across the reshape() function, but I have been very successful using it. As an example: ll <- data.frame(time.1=seq(0.1,1,length=10),L1=rnorm(10),time.2=c(seq(1,7),rep(NA ,3)),R2=c(rnorm(7),rep(NA,3)),time.3=c(seq(1,1.5,length=3)+5,rep(NA,7)),M3=c (rnorm(3),rep(NA,7))) which

Dear R-helpers, I have some problem with a model I am trying to run with nlmeODE. I have a covariate that I want to include in the model and I can't find in the documentation how to include it in the model so that it can calculate it. >formula(trial) PME ~ time | Ident abcnoniso <- list(DiffEq=list( debdt = ~ -ksol*pH*eb , desdt = ~

A couple of months ago I downloaded R 1.6.0's binary for Windows, which I have since been using without any particular problem, that is on one machine running Windows 98. Last week I set up a local network and have since tried to use R remotely, i.e. R is installed on the previous machine and nowhere else. On the client's side, the GUI opens correctly and allows me to do simple

Dear Colleagues, Our group is also working on implementing the use of R for pharmacokinetic compartmental analysis. Perhaps I have missed something, but > fit <- nls(noisy ~ lsoda(xstart, time, one.compartment.model, c(K1=0.5, k2=0.5)), + data=C1.lsoda, + start=list(K1=0.3, k2=0.7), + trace=T + ) Error in eval(as.name(varName), data) : Object

Dear R-experts I would like to read all files from a directory, the files have names "myname0001.txt" etc. I paste the directory plus file names and use "read.delim()". My problem is that some file names are missing, so I get an error and my program stops. Is there a way to check for a null pointer analogous to C, so that I can simply skip non-existent filenames? Please do

Does someone out there have an example of R-code for a split-split plot ANOVA using aov or another function? The design is not balanced. I never set up one in R before and it would be nice to see an example before I tackle a very complex design I have to model. Thanks, Mike Mike Saunders Research Assistant Forest Ecosystem Research Program Department of Forest Ecosystem Sciences University of

Hi all I am doing some stats work for a group of biologists who require windows metafiles (*.wmf) for their publications. To create these, I appear to have two choices: 1. Restart my machine in Windows and use savePlot 2. Keep my machine in linux, save as another format, then convert. I'd rather stay in linux; but how do I get wmf files? I looked at using ImageMagick's convert,

Hello R-Help, I've got a bit of an issue with WMF's. I am working on WindowsXP and outputting WMF format images. I then take the WMF format images and insert them in PowerPoint. I take the PowerPoint and convert it to PDF. The WMFs are nothing special. Just the typical x-y plot with a gray dotted grid added and a few matlines. The WMFs, when produced look great! They show

For the first error message: > win.metafile(file = "//.../plot1.wmf", + width = 8.5, height = 6.25) > lset( list( background = list(col = "white"))) Error in get(x, envir, mode, inherits) : variable "win.metafile://.../plot1.wmf" was not found > traceback() 4: get(device) 3: trellis.device(device = .Device, new = FALSE) 2: trellis.par.get(item) 1:

Hi, I recently encountered an error using the command "savePlot" when trying to save the third of 3 open graph windows. After successfully saving and closing the first two windows I receive an "Invalid device number in savePlot" error. The following is copied and pasted from an example session to illustrate the behaviour: > windows() > hist(rnorm(100)) > windows()

Hi, I have been having issues producing plots in JPEG format, using type = "cairo" to get better anti-aliasing. When trying to set the physical size with units = "cm" or units = "mm", the width is set correctly but the height is not - it looks like the height is simply treated as pixels regardless of the 'units' argument. Example: x = 1:10 y = 2 * x

gday R friends, I was wondering if anyone could help me with some graphics issues please. I'm processing R in batch on redhat6.2 with some graphics to be imported into msword on NT. For scaling within my doco I want some vector type graphics. win.metafile doesn't seem to be supported on redhat which gives me postscript. postscript is ideal but doesn't contain a preview so I

Dear all, Has anyone written functions to compute (and ideally project / plot onto a plane) confidence ellipsoids for q linear model parameters, q > 2 ? Greetings and best wishes, Renaud -- Dr Renaud Lancelot, v?t?rinaire CIRAD, D?partement Elevage et M?decine V?t?rinaire (CIRAD-Emvt) Programme Productions Animales http://www.cirad.fr/presentation/programmes/prod-ani.shtml ISRA-LNERV