search for: oates

Dear R-Community! The example "oats" in MASS (2nd edition, 10.3, p.309) is calculated for aov and lme without interaction term and the results are the same. But I have problems to reproduce the example aov with interaction in MASS (10.2, p.301) with lme. Here the script: library(MASS) library(nlme) options(contrasts = c("contr.treatment", "contr.poly")) # aov: Y ~

Dear List, I have been using R to create an xyplot using the panel function within lattice libraries. This plot is based on the data supplied in R named 'Oats'. The graph represents oat yield by nitro level with an overlay of each variety of oats for each nitro level. I have three questions regarding this graph: 1) I cannot seem to specify the type of symbol used by the plot, even though

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, ...) { /

Full_Name: Kevin Wright Version: 1.8.0 OS: Windows 95 Submission from: (NULL) (170.54.59.160) This bug also happens to me using R 1.9.0 on Windows 2000. Took me a long time to create a reproducible bug, but I think I have succeeded. I suspect my test function has a bug, but I don't see anything wrong. Plus, the bug only shows up when selecting certain filenames. Nearest I can speculate,

The following code produces a legend ("key") that mentions the unused levels of Block. library(MEMSS) xyplot(yield~nitro, subset=(Block=="I" | Block=="II"), data=Oats, group=Block, auto.key=T) and adding "drop.unused.levels=T" does not fix it. And in fact even the following does not solve the problem: xyplot(yield~nitro,

On Thu, Mar 29, 2018 at 12:48 PM, Asif Iqbal <vadud3 at gmail.com> wrote: > > > On Thu, Mar 29, 2018 at 7:21 AM, Steven Tardy <sjt5atra at gmail.com> wrote: > >> A STATEFUL firewall with ?ip any any? can and will still block asymmetric >> communications due to the firewall keeping track of state (hence tha name >> stateful firewall). >> >>

Hi! I am trying to implement the Morris One at a Time (OAT) sensitivity analysis technique in R using the package 'sensitivity'. The OAT and similar other techniques in the package requires data frame X (which stores the design or parameter combinations) to be populated by some function in the package such as fast99 or morris. These functions take only the range and step size of the

Hello, Are there any access functions to the various residual variables that should result from a repeated measures ANOVA ? MyAOVObject$residuals does not exist, and simply printing MyAOVObject gives a very long print of all fields in the result list, many of which I can't see what they are exactly : $error.qr$qraux, for instance. What I would like basically is to inspect those residuals

I am trying to do Tukey HSD comparisons on a repeated measures expt. I found the following example on r-help and quoted approvingly elsewhere. It is broken. Can anyone please tell me how to get it to work? I am using R 2.4.1. > require(MASS) ## for oats data set > require(nlme) ## for lme() > require(multcomp) ## for multiple comparison stuff > Aov.mod <- aov(Y ~ N + V +

A STATEFUL firewall with ?ip any any? can and will still block asymmetric communications due to the firewall keeping track of state (hence tha name stateful firewall). Tcpdump on your servers /other/ NICs and you?ll see the tftp traffic leaving your server on some other NIC (probably on with the default route). The upstream firewall will then block the tftp response if it never saw the tftp

When we call a lattice function such as xyplot, to what extent does the "data" designation cause the function to look inside the "data" for variables? In the examples below, the "subset" argument understands that "Variety" is a variable in the data. But the "scales" argument does not understand that "nitro" is a variable in the data.

Dear all, Consider a completely randomized block design (let's use data(Oats) irrespoctive of the split-plot design it was arranged in). Look: library(nlme) fit <- lme(yield ~ nitro, Oats, random = ~1|Block, method="ML") fit2 <- lm(yield ~ nitro + Block, Oats) anova(fit, fit2) gives this: Model df AIC BIC logLik Test L.Ratio p-value fit 1 4 624.3245

I am using R in my undergraduate engineering statistics courses, in which I am currently discussing one-way analysis of variance. Many of the data sets for the exercises or examples give the data in the form on one column per treatment level. To use the aov or lm functions such data should be converted to one column with all the response observations and a companion column with indicators of the

Dear R users, I have a question about the patterned variance-covariance structure for the random effects in linear mixed effect model. I am reading section 4.2.2 of "Mixed-Effects Models in S and S-Plus" by Jose Pinheiro and Douglas Bates. There is an example of defining a compound symmetry variance-covariance structure for the random effects in a split-plot experiment on varieties of

On Thu, Mar 29, 2018 at 7:21 AM, Steven Tardy <sjt5atra at gmail.com> wrote: > A STATEFUL firewall with ?ip any any? can and will still block asymmetric > communications due to the firewall keeping track of state (hence tha name > stateful firewall). > > Tcpdump on your servers /other/ NICs and you?ll see the tftp traffic > leaving your server on some other NIC (probably

have you checked that tftp is added to hosts.allow. syslog may be reporting libwrap errors, libwrap is trcpwrappers regards peter On 11 April 2018 16:57:04 "Asif Iqbal" <vadud3 at gmail.com> wrote: > On Thu, Mar 29, 2018 at 12:48 PM, Asif Iqbal <vadud3 at gmail.com> wrote: > > > > > > > On Thu, Mar 29, 2018 at 7:21 AM, Steven Tardy <sjt5atra at

Hi use Rs, I have a csv file: "1989-90","1990-91" Barley,23,34 Oats,15,16 Which I want to turn into: year, Barley, Oats 1 "1989-90", 23, 15 2 "1990-91",34,16 Transpose doesn't quite do it, is there a standard way? Cheers, Geoff Russell

Hello, can anybody tell me how to produce a plot like the one in http://cran.r-project.org/web/packages/lme4/vignettes/Implementation.pdf on page 13, Figure 6? The data is stored in: library(nlme) data(Oats) Cheers -- View this message in context: http://r.789695.n4.nabble.com/Help-needed-in-reproducing-a-plot-tp4119603p4119603.html Sent from the R help mailing list archive at

Hi all, I have a problem in that when I plot points that have a high degree of precision, some significant rounding seems to occur, resulting in uneccessary overlap of my points. Is there a way to specify the resolution or precision in plotting functions? Is there an underlying grid I have to modify somehow? Many Thanks, Jon -- Mares eat oats and does eat oats and little lambs eat ivy. A

Oats[Oats$Variety %in% c("Victory", "Golden Rain"),] or subset(Oats, Variety %in% c("Victory", "Golden Rain")) Stefano -----Messaggio originale----- Da: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch]Per conto di elyakhlifi mustapha Inviato: luned? 23 aprile 2007 9.56 A: R-help at stat.math.ethz.ch Oggetto: [R] extract from