search for: oat

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(&quot...

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 it is included in the code below, it will change when the cod...

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, ...) { / > / # How to normalize my heatmap metric > with /...

...when selecting certain filenames. Nearest I can speculate, choose.files becomes confused after a filename that has capital letters and numbers has been chosen. Sounds crazy and I'd like to know more, but choose.files is mostly .Internal code. First, create two CSV files. library(nlme) data(Oats) write.table(Oats,"c:/oats.csv",sep=",") write.table(Oats,"c:/C20oats.csv",sep=",") Quit R and restart. Follow the steps below. Select oats.csv the first 4 times, then select C20oats.csv and then select oats.csv . Notice how the filename is corrupted the...

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, data=Oats[Oats$Block=="I" | Oats$Block=="II",], group=Block, auto.key=T) The following workaround solves it,...

...pdump: verbose output suppressed, use -v or -vv for full protocol decode > listening on any, link-type LINUX_SLL (Linux cooked), capture size 262144 > bytes > > 16:40:08.390939 IP 192.168.1.10.35553 > 192.168.1.20.69: 16 RRQ "file" > netascii > E..,J1 at .>..n./...oAt...E..#...file.netascii................... > 16:40:13.391133 IP 192.168.1.10.35553 > 192.168.1.20.69: 16 RRQ "file" > netascii > E..,N. at .>..../...oAt...E..#...file.netascii................... > 16:40:18.391220 IP 192.168.1.10.35553 > 192.168.1.20.69: 16 RRQ "...

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

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 + Error(B/V), data = oats) > Lme.mod <- lme(Y ~ N + V, random = ~1 | B/V, data = oats) > summary(Aov.mod) > anova(Lme.mod) > summary(Lme.mod) > summar...

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

...;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. What principle is at work? library(MEMSS) # The following works fine: xyplot( yield ~ nitro , data=Oats , scales=list( x=list( at=unique(Oats$nitro) ) ) , subset=Variety=="Victory" ) # But the following returns an error: xyplot( yield ~ nitro , data=Oats , scales=list( x=list( at=unique(nitro) ) ) ) Thanks for any insight Jacob A. Wegelin Assistant Professor De...

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

...this. Here is an example > library(Devore5) > data(ex10.09) > boxplot(ex10.09) > str(ex10.09) # show the structure `data.frame': 6 obs. of 4 variables: $ Wheat : num 5.2 4.5 6 6.1 6.7 5.8 $ Barley: num 6.5 8 6.1 7.5 5.9 5.6 $ Maize : num 5.8 4.7 6.4 4.9 6 5.2 $ Oats : num 8.3 6.1 7.8 7 5.5 7.2 > stack <- function(data) + data.frame(values = unlist(data), + lev = factor(rep(names(data), lapply(data, length)))) > stack(ex10.09) # try it on the example values lev Wheat1 5.2 Wheat Wheat2 4.5 Wheat Wheat3 6.0 Wh...

...ffects 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 oats. I ran the codes from the book and extracted the variance and correlation components: >library(nlme) >data(Oats) >fm4OatsB <- lme(yield~nitro, data=Oats, random=list(Block=pdCompSymm(~Variety-1))) >VarCorr(fm4OatsB) Block = pdCompSymm(Variety - 1) Variance StdDev...

...-nniany host 192.168.1.10 tcpdump: verbose output suppressed, use -v or -vv for full protocol decode listening on any, link-type LINUX_SLL (Linux cooked), capture size 262144 bytes 16:40:08.390939 IP 192.168.1.10.35553 > 192.168.1.20.69: 16 RRQ "file" netascii E..,J1 at .>..n./...oAt...E..#...file.netascii................... 16:40:13.391133 IP 192.168.1.10.35553 > 192.168.1.20.69: 16 RRQ "file" netascii E..,N. at .>..../...oAt...E..#...file.netascii................... 16:40:18.391220 IP 192.168.1.10.35553 > 192.168.1.20.69: 16 RRQ "file" netascii...

...ed, use -v or -vv for full protocol decode > > listening on any, link-type LINUX_SLL (Linux cooked), capture size 262144 > > bytes > > > > 16:40:08.390939 IP 192.168.1.10.35553 > 192.168.1.20.69: 16 RRQ "file" > > netascii > > E..,J1 at .>..n./...oAt...E..#...file.netascii................... > > 16:40:13.391133 IP 192.168.1.10.35553 > 192.168.1.20.69: 16 RRQ "file" > > netascii > > E..,N. at .>..../...oAt...E..#...file.netascii................... > > 16:40:18.391220 IP 192.168.1.10.35553 > 192.168.1.2...

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 Nabble.com.

...ts 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 kid'll eat ivy too, wouldn't you? -- Leland Palmer, from Twin Peaks ########### Jonathan Manning PhD Student Barton Group School of Life Sciences University of Dundee Scotland, UK Tel +44 1382 388707

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