similar to: on trellis.par.set/get (reproducing figures from Pinheiro & Bates)

Hi R users, I'm here trying to understand correlated residuals in nonlinear estimation. I'm reading/studying the book Bates, D. M. and D. G. Watts, (1988), /Nonlinear regression analysis and its applications/, Wiley, NY. pages 92-94, trying to reproduce the figures and to find out the code in R to perform the necessary calculations. I also consulted Pinheiro and Bates, but without

Dear R users, today I've got the following problem. Here you are a dataframe as example. There are some SAMPLES for which a CONCentration was recorded through TIME. The time during which the concentration was recorded is not always the same, 10 points for Sample A, 7 points for Sample B and 11 for sample C Also the initial concentration was not the same for the three samples. I would like

Dear R-users, consider the two following outputs, ## 1 and ## 2 \begin{Scode}{Setup, echo = FALSE, print = FALSE, eval = TRUE} with(expand.grid(Fact1 = 1:3, Fact2 = 1:40), table(Fact1, Fact2)) ## 1 xtable(with(expand.grid(Fact1 = 1:3, Fact2 = 1:40), table(Fact1, Fact2))) ## 2 \end{Scode} The first line with(expand.grid(Fact1 = 1:3, Fact2 = 1:40), table(Fact1, Fact2))

Dear R-users, I have the following problem In a lab experiment I have to mix three solutions to get different concentrations of various molecules in a cuvette I've used R to calculate the necessary ?liters for each of the level of the experiment and I must confess that it is more useful and easier to achieve the results than using spreadsheets. But there's a problem. Imagine that

Dear R users, I have the following two problems, related to the function sub, grep, regexpr and similia. The header of the file(s) I have to import is like this. c("y (m)", "BD (g/cm3)", "PR (Mpa)", "Ks (m/s)", "SP g./g.", "P (m3/m3)", "theta1 (g/g)", "theta2 (g/g)", "AWC (g/g)") To get rid of spaces and

Hi there, this is my first post to the list. I'm a newcomer to the linux world. I started using it on a regular basis on May and I must say I'm quite comfortable with it, even if I have to re-learn a lot of things. But this is not a problem, I will improve my knowledge with time. My main problem now is that I bought myself a new computer (described at the end of this e-mail) and I was

Dear R-users. I eventually bought myself a new computer with the following characteristics: Processor AMD ATHLON 64 DUAL CORE 4000+ (socket AM2) Mother board ASR SK-AM2 2 Ram Corsair Value 1 GB DDR2 800 Mhz Hard Disk WESTERN DIGITAL 160 GB SATA2 8MB I'm a newcomer to the Linux world. I started using it (Ubuntu 7.10 at work and FC4 on laptop) on a regular basis on May. I must say I'm

Dear R-users. I eventually bought myself a new computer with the following characteristics: Processor AMD ATHLON 64 DUAL CORE 4000+ (socket AM2) Mother board ASR SK-AM2 2 Ram Corsair Value 1 GB DDR2 800 Mhz Hard Disk WESTERN DIGITAL 160 GB SATA2 8MB I'm a newcomer to the Linux world. I started using it (Ubuntu 7.10 at work and FC4 on laptop) on a regular basis on May. I must say I'm

Dear R-Users, I'm currently studying the book Statistical model in S by Chambers J.M and Hastie T.J.. At page n 3 there's a plot showing the means of a variable at each of the levels of the factors of an experiment. I hope to be able to reproduce it here by ASCII art. B6 S | | Thin | | | | | | B3 __|_____|_____|__......etc | |

Dear R users, this is a follow up of this message http://tolstoy.newcastle.edu.au/R/e6/help/09/05/13897.html I'm reproducing the core of it for convenience. > // > / data(Oats, package = "MEMSS") / > / tp1.oats <- xyplot(yield ~ nitro | Variety + Block, / > / data = Oats, / > / panel = function(x, y, subscripts, ...) { /

Dear R-users, my present problem is related to lattice. I would like to put some text in each panel, namely a p-value. I therefore wrote a simple panel function as reported here below. I'm able to write the value in each panel at the maximum value of y for each panel, but this obviously overlay the text to the points. What I'm looking for is to write the text always in the same

Dear R-users, the problem I'm facing today is to write a multicolor main title in a plot. I would like to have 3/4 colors to help people in the audience to visually realize what is what in a 4x2 plot array On top of each plot I would like to have two lines. The 1st line will indicate what is plotted The 2nd line will indicate a a reference Graphically speaking I would like a result similar

Dear R-users, after seeking for help in R-search I did not find any hint on my particular problem. Countless help on "true" overlay, but nothing on this. Please consider the following: par(mfrow = c(2, 1)) T <- seq(0, 20, by = 0.01) ## PLOT 1 plot(T, 30*exp(-0.65*T), type = "l" ) points(T, 30*exp(-0.26 * T), type = "l", lty = "F8") points(T,

Kevin Wright wrote: > library(nlme) > m2 <- gnls(conc ~ t1*(1-t2*exp(-k*time)), > data = df.Chloride, > start = list( > t1 = 35, > t2 = 0.91, > k = 0.22)) So my error was to use nls instead that gnls. Thanks a lot, Kevin. > summary(m2) > plot(m2) > lag.plot(resid(m2), do.lines=FALSE) >

Dear R users, my problem today deals with my ignorance on regular expressions. a matter I recently discovered. Consider the following foo <- c("V_7_101110_V", "V_7_101110_V", "V_9_101110_V", "V_9_101110_V", "V_9_s101110_V", "V_9_101110_V", "V_9_101110_V", "V_11_101110_V", "V_11_101110_V",

Dear R-helpers, the problem I'm facing today is to convince lattice to paint some areas in gray. The areas I would like to have in gray, are confidence bands I've googled around in the mailing list archives and eventually find some clues. This link is my starting point http://tolstoy.newcastle.edu.au/R/e2/help/07/04/15595.html I'm reproducing here the code for your convenience est

Deae R helpers, the problem I'm facing today is related to the manipulation of a string. The string is coming from a a porosimeter, whose control is under a complicate set-up of two computers One (running on DOS) is controlling directly the hardware, while the other (running on win XP) which process the data is connected to the first one via serial. I get the data from this last computer.

Dear R-users, please consider the following minimal example: \documentclass[a4paper,titlepage,onecolumn,12pt]{article} \usepackage[italian]{babel} \usepackage{amssymb} \usepackage[utf8x]{inputenc} \usepackage[pdftex]{graphicx} \begin{document} <<label=test, echo=FALSE, results=tex>>= df.data1 <- cbind.data.frame(A = rnorm(18), B =factor(rep(LETTERS[1:6],

Hi all, I'm trying to work out the following problem, but I can't imagine how. I have the following (much reduced & oversimplified) dataset My.df <- cbind.data.frame(PPM=c(15.78, 15.81, 15.87, 15.83, 15.81, 15.84, 15.91, 15.90, 15.83, 15.81, 15.93, 15.83, 15.70, 15.92, 15.76, 15.81, 15.91, 15.75, 15.84, 15.86, 15.82, 15.79,

Hi all, first of all thanks for the answer to my previous question on lattice some time ago. In particular to Patrick Connolly for advices on netiquette (I hope this time I'm doing right....) and to Deepayan Sarkar fro the help on lattice. Now, my nowaday problem. Please consider the following mydf<-cbind.data.frame( RESPONSE = c(rnorm(9,rep(2:4,each=3),10),