similar to: Plot.lme error

Dear all, In the SASmixed package there is an example of an analysis of a split-plot experiment. The model is fm1Semi <- lme( resistance ~ ET * position, data = Semiconductor, random = ~ 1 | Grp) where Grp in the Semiconductor dataset is defined as ET*Wafer. Is it possible to specify the grouping directly some way, e.g. like fm1Semi <- lme( resistance ~ ET * position, data =

I am trying to understand the proper way to encode the nesting structure for data in the context of nlme, in the specific case of individuals nested within species for which each individual is unique. I have searched through Pinheiro & Bates and also past postings, but without success. Take the Wafer data set which has 2 levels: Wafer (8 values) and Site nested within Wafer (10 values for

Hi helpeRs, I've been working through some examples in Pinhiero & Bates( 2000) trying to understand how to translate to the new Lme4 syntax but without much luck. Below is what I think I should do, but either the answers don't come out the same or I get errors. In the Oxide problems I'm particularly interested in obtaining the levels coeficients but this options no longer seems

I have one fixed effect, sor, with two levels. I have eight lots and three wafers from each lot. I have included the data below. I would like to fit a mixed model that estimates a covariance parameter for wafer, which is nested in lot, and two covariance parameters for lot, one for each level of sor. The following command fits the model that I want, except for it estimates the correlation

All, I am having trouble with a "read.table()" function that is inside of another function. But if I call the function by itself, it works fine. Moreover, if I run the script on a Mac OS X (with the default Mac OS X version of R installed, rev 2.8), it works fine. But it does not work if I run it on windows vista (also default Windows version of R, rev. 2.8). Again, both

My question relates to problems that I'm having matching lme and lmer examples in P&B. using Matix 0.995 In the Oxide example in p167-170 I can't get the level 2 coefficient estimates to match the fm1Oxide model in lme is data(Oxide,package="nlme") lme(Thickness~1,Oxide) which I translate in Lmer syntax to fm3Oxide<-lmer(Thickness~

Dear all, I am planning to use nlme library for analysis of experiments in semiconductor industry. Currently I am using "lm" but plan to move to "lme" to handle within wafer / wafer-to-wafer and lot-to-lot variation correctly. So far everything is working well, but I have a fundamentel question: NLME offers "maximum likelihood" and "restricted maximum

Dear all, Can anyone help me, my R software can not run a nested linear regression by using the lme funcion. The message that appears is  Error: could not find function "lme" I already downloaded and loaded the package, please see below. Thank you in advance for any help! Nadia. > data<-read.csv("/Users/nadiasan1/Desktop/MOE and MOR.csv")> attach(data)>

Hi all, I just wonder what singular convergence means. Thanks. Yen Lin Error in lme.formula(Data ~ 1, random = ~1 | Wafer/fie/loc, subset = Wafer == : singular convergence (7) [[alternative HTML version deleted]]

I have created a tree and want to save some of the data so that I can create a html table from it. I would like to save the output from data.ltr (see example below) to a file, but haven't found a way to do that, keeping the nice format that typing data.ltr gives me (see output below). Is there a way to do this? Example: library (maptree) library (tree)

Hello all, here's a real-world example: I'm measuring a quantity (d) at five sites (site1 thru site5) on a silicon wafer. There is a clear site-dependence of the measured value. To find out if this is a measurement artifact I measured the wafer four times: twice in the normal position (posN), and twice rotated by 180 degrees (posR). My data looks like this (full, self-contained code at

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

Paul: I may have found the issue (which is similar to your conclusion). I checked using egsingle in the mlmRev package as these individuals are strictly nested in this case: library(mlmRev) library(nlme) fm1 <- lme(math ~ year, random=~1|schoolid/childid, egsingle) fm2 <- lmer(math ~ year +(1|schoolid:childid) + (1|schoolid), egsingle) Checking the summary of both models, the output is

Hello R-Users! Does anyone know of a package to draw/analyze silicon wafer maps? Here are some examples http://www.java2s.com/Code/Java/Chart/JFreeChartWaferMapChartDemo.htm http://dp.pdf.com/site/products/wafermap/binmap.html I've google'd and searched CRAN with no luck. It seems possible for R, given the hexbin and hist2d graphs I saw at the Graph gallery.

Hi Timo, >> * LIST (\HasNoChildren) "/" "folder1" >> * LIST (\HasNoChildren) "/" "folder2" >> * LIST (\HasNoChildren) "/" {9} >> six"wafer >> . OK List completed. >> >> Note the {9} length of the following real folder name. Is this normal >> handling of special folder names ? >Yes. Well,

Hello everyone, I am ploting a groupeddata object with formula: formula(mydatausegroup) BF ~ HO | ID/Infar/Day Using this command: plot(na.omit(mydatausegroup), displayLevel=2,layout=c(10,2),aspect=2) This trellis plot does almost what I want and produces a 10x2 trellis plot, each panel is labeled as ID/Infar where infarct is either 1 or 0. And in each panel, it plots BF vs HO for each Day.

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.

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,

Warnings about /metno/internweb/fouweb/htdocs/teknisk/verktoy/R/faq.html BAD LINK: http://stat.ethz.ch/R/manual/ BAD LINK: http://nlme.stat.wisc.edu/MEMSS/ BAD LINK: http://stats.mth.uea.ac.uk/Rcgi/ BAD LINK: http://software.biostat.washington.edu/statsoft/ess/ ------------------------------------------------------------ Warnings about

Hi everyone, Using R 2.10.1 on Windows Vista. DOWNLOADED PACKAGES DO NOT INSTALL. I expect to see them in C:\Program Files\R\R-2.10.1\library These files download (see below). And they are all in zip format. What am I doing wrong? Please help. All suggestions appreciated. trying URL 'http://cran.cnr.Berkeley.edu/bin/windows/contrib/2.10/RColorBrewer_1.0-2.zip' Content type