similar to: [Fwd: Re: Several simple but hard tasks to do with R]

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

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

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, I would like to exactly reproduce a figure like the 1.5 or 1.9 or 4.13 from the book Mixed effects models in S and S-Plus. Not for the sake of it, but because I have my own data I would like to plot in that fashion (no colors) If I write plot(ergoStool) I can get a good informative plot with colors, but I would like to have a B&W one instead. I've played a little with

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

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

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

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

Hi all, I would like to add to each panel of a bwplot a coloured central band, centered on the mean of the values, being its width +- 2% of the mean itself. I know how to add lines, i.e. something like bwplot(X ~ Y|FACTOR data=my.df, panel= function(x, y){ panel.bwplot(x, y) panel.abline(v = mean(x, na.rm = T) - mean(x, na.rm = T) * 0.02 panel.abline(v =

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,

Dear all, A hopefully simple question: how do I round a series of values (held in an object) to the nearest 5? I've checked out trunc, round, floor and ceiling, but these appear to be more tailored towards rounding decimal places. Thanks, Steve _________________________________________________________________ icons.

Dear R-users: I have some problems working with lme function, and i would be glad if anyone could help me. this kind of analysis i was used to do with PROC MIXED from SAS, but i would like to move to R, for many reasons... So, the problem is: Imagine the I have 3 factors: fact_A, fact_B and fact_C: The latter I would assume that is random, and the rest of them are fixed. Analysing the

Hi everyone! I've recently begun to learn R for my job as the IT department suffers from lack of funding for new software. I was talking to the guy in charge of Requisitions and have found out the budget for books is in great shape. So, I'm curious what books people know of that have R examples and are good for: 1.) Uni and Multivariate Time Series Analysis/Forecasting 2.) GLMs (at