Re: [R] References verifying accuracy of R for basic statisticalcalculations and tests

2006-07-14 Thread Jesus Frias
Hi all, There is a NISTnls package with the nonlinear regression examples in CRAN. This work might have been done already. regards, Jesus On Fri, July 14, 2006 10:57 am, Sean O'Riordain said: Please don't shoot! q: would it be a good idea to use these datasets as a basis for some

[R] Equal number of cuts in a contourplot with trellis

2006-01-12 Thread Jesus Frias
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) +

RE: [R] Split-split plot ANOVA

2005-02-01 Thread Jesus Frias
Hi Mike, *An example of the use of aov() for a split-plot is in MASS library(MASS) example(Oats) The book also gives a detailed explanation *pp 45-52 of the Pinheiro and Bates book gives you an example of the use of lme() on a split-plot. If you have a non balanced

RE: [R] Nonlinear weighted least squares estimation

2004-11-08 Thread Jesus Frias
Hi Robert You can try gnls() in the nlme package or if that is not suitable, try to follow the example in the nls() help page for weighted regression that has an example of weighted regression from MASS. from the nls help page: ## weighted nonlinear regression Treated -

[R] L term in anova.lme

2004-06-01 Thread Jesus Frias
Dear R-helpers, Could somebody provide an example of the use of the L term in the anova.lme() function to test linear combinations of different terms in a lme model?. Is it possible to make more than one test in a single call to anova.lme()?. best regards, Jesus

RE: [R] Metafiiles into Word R 1.9.0

2004-05-26 Thread Jesus Frias
-- -Original Message- From: Jason Turner [mailto:[EMAIL PROTECTED] Sent: 26 May 2004 22:11 To: Jesus Frias Subject: RE: [R] Metafiiles into Word R 1.9.0 There is one more thing that might shed some light (or not) on it: I have just found that I can actually copy-paste a plot

[R] Metafiiles into Word R 1.9.0

2004-05-24 Thread Jesus Frias
Dear R-helpers, I recently upgraded to R 1.9.0 in my computer at work and at home: 1.-The computer at home has Windows XP and Office XP and it seems to work perfectly and I copy-paste graphics perfectly. 2.-The computer at work has Microsoft Windows 2000 (5.00.2195

[R] Change default values of a function

2004-05-12 Thread Jesus Frias
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().

[R] Reshape of repeated measures data

2004-05-05 Thread Jesus Frias
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

RE: [R] Savitzky-Golay smoothing for reflectance data

2004-02-05 Thread Jesus Frias
Hi, Quote: The Savitzky-Golay smoothing turns out to be exactly equivalent to fitting data to a polynomial as described (Skoog, Holler and Nieman, Principles of Instrumental Analysis, pp 111) While R is not a tool for instrumental analysis, I am sure you can find

RE: [R] Fitting compartmental model with nls and lsoda?

2004-01-23 Thread Jesus Frias
Now, I am including Richar Upton's 2 cm model implementation and Christoffer Tornoe's nls solution (I recommend Christoffer's nlmeODE package for these problems also if multi-response data is available) The code follows: Multi-response data?, Is there a way of dealing with those

RE: [R] extracting the levels of a subset of data

2003-09-19 Thread Jesus Frias
Hiya I'd like to get only A as output... two solutions: 1.-use unique() unique(tempdata$TYPE) to get your answer 2.-or call factor again to reduce the levels: tmpd - subset(myd,clas==A) levels(tmpd$clas) [1] A B C tmpd$clas - factor(tmpd$clas) levels(tmpd$clas) [1] A regards,

RE: [R] reading non-existent files

2003-01-28 Thread Jesus Frias
Hi Christian, I suppose that you are doing this in a loop and is stopping. Have a look at the function try() ?try Try an Expression Allowing Error Recovery. Description: `try' is a wrapper to run an expression that might fail and allow the user's code to handle