similar to: an(other) anova question

On Sun, 10 Oct 2004, RenE J.V. Bertin wrote: > Could I ask some hints/help in interpreting the following ANOVA results, > please? This concerns an experiment where I study the incidence and > severity of motion sickness. I have Sickness.norm, a subjective > discomfort/sickness estimate, normalised to 0..1, the session time T > (normalised to 0..1 and binned in 0.2 wide bins) and a

I believe you are right, but can you please explain why anyone would want to fit this model? It differs only in the coding from aov(y ~ a + b + Error(c), data=test.df) and merely lumps together the top two strata. There is a much simpler fix: in the line if(intercept) nmstrata <- c("(Intercept)", nmstrata) remove the condition (and drop the empty stratum later if you

I think there are two bugs in aov() that shows up when the right hand side of `formula' contains both `-1' and an Error() term, e.g., aov(y ~ a + b - 1 + Error(c), ...). Without `-1' or `Error()' there is no problem. I've included and example, and the source of aov() with suggested fixes below. The first bug (labeled BUG 1 below) creates an extra, empty stratum inside

Hi, I have a problem with aov(). I used it many times, but now I have new data, tried to use it in the same way. I get I strange Error message that i can't understand (see below). I guess it is caused by incompatible data types. The dataframe also looks ok to me. I use: spk

Dear list members, I am analysing my microarray data using limma package. Now I encounter several problems. Looking forward to your suggestions! Question 1: During the process of background correction using method="normexp", four warning messages appeared as "NaNs produced in: log(x)" (as you can see in the program posted below). What does that mean? How will it effect the

Dear All, I am analysing a dataset on levels of herbivory in seedlings in an experimental setup in a rainforest. I have seven classes/categories of seedling damage/herbivory that I want to analyse, modelling each separately. There are twenty maternal trees, with eight groups of seedlings around each. Each tree has a TreeID, which I use as the random effect (blocking factor). There are two

Prof Brian Ripley <ripley@stats.ox.ac.uk> writes: > I believe you are right, but can you please explain why anyone would want= =20 > to fit this model? When analysing data from mixture designs, the variables add up to a constant, and it is often preferrable to fit models without intercept term. We often have experiments where a mixture design (typically raw materials) is combined

Dear R colleagues: Am trying to fit a simple NL model to determine Economical Optimum Nitrogen Rates. The segmented (quadratic + plateau) model only works with some y's, in some cases I get a "singular gradient" error. I'll appreciate any ideas in how to solve the singular gradient error. Thanks, Jose # The following code works using yield2 in the nls model but not using

Hi list! While regular backup solutions like amanda or bacula are very good at their job, ie keeping point2point copies of the files currently on disk, I find them less suited for archiving - having unused files move to tape in duplo and stay there until requested. I've even read of multi-tier solutions - move to slower disks after a week and further on to tape after a month. Does anyone

Hi, OK, I will describe it in more detail.In the meantime I installed the newest debian package of r-base etc. and the error message changed from:Error in "names<-.default"(`*tmp*`, value = nmstrata) : names attribute must be the same length as the vectorto:Error model is singular in: aov(rati ~ cs * spk + Error(sub/(cs * spk)), data = rating.df) The model is a repeated

Hello, I'm attempting to find multiple breakpoints in an association of my response variable (R.AUC) with two explanatory variables ('s.size' and 'bedekking'). The association between 's.size' and 'R.AUC' shows a plateau, but the value when this plateau is reached is differs for different values of 'bedekking'. Initially I thought these different

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 This package implements tools for manipulationg digital images and the Propagation Separation approach for smoothing digital images. Version 0.6.1 is a major update with many new features: Input/Output: - - we now keep exif information as comment of an image ~ this information is evaluated to recover the original adimpro-object ~ used to produce

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 This package implements tools for manipulationg digital images and the Propagation Separation approach for smoothing digital images. Version 0.6.1 is a major update with many new features: Input/Output: - - we now keep exif information as comment of an image ~ this information is evaluated to recover the original adimpro-object ~ used to produce

Hi there are some predefined self-start functions, like SSmicmen, SSbiexp, SSasymp, SSasympOff, SSasympOrig, SSgompertz, SSflp, SSlogis, SSweibull, Quadratic, Qubic, SSexp (nlrwr) Btw, do you know graphic examples for this functions? The SSexpDecay (exponential decay) for y = (y0 - plateau)*exp(-k*x) + plateau from

I am simulating sickness among a group of families. Part of the task is to randomly draw who in the family will be sick, randomly drawing from family ID's where Dad =1, Mom = 2, Kid1 = 3, Kid2 = 4., etc. My census of Dads is of the form shown below. Dad_ID Spouse (Y=1;N=0) #Kids #People_Becoming_Sick 1 1 0 1 2

Dear All, I would like to make a quadratic with a plateau model in R. Is there a package in R doing this? The bentcableAR package seems won't work. The link below describes what I am looking for in R exactly: http://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_nlin_sect033.htm -- Thanks so much! Orange help.ly2005@gmail.com [[alternative HTML version

Hi, I'm interested in doing a repeated measures anova using aov. The procedure is nicely described in section 6.7.1, pp. 24-27 of Baron and Li's "Notes on the use of R for psychology experiments and questionnaires," and I've reproduced their example exactly. My own problem is almost identical to theirs: rawdat<-c(1.6530074e+001, 1.2124254e+001, 1.0040371e+001,

Dear Lists, Here are my latest posts on retrenching and locating facilities better. Do you have anything in academic papers that use data mining- analytics for 1) Better Employee ROI quantitative measurement or cutting costs rationally based on data rather across the board. 2) Locating facilities based on pure quant based approaches than local supply demand conditions. Specific application for

Hello everyone, I working in a public health project and we have created a Decision Tree for categorical variables usign the package rpart. Our goal is to develop a model (Using the ROC tool) in order to predict presence/ausent of  diabetes and get a better understanding  of what are the important factors in a particular chilean population.  There are some importants variable that we have found. 

Dear All, I am doing experiments in live plant tissue using a laser confocal microscope. The method is called "fluorescence recovery after photo-bleaching" (FRAP) and here follows a short summary: 1. Record/ measure fluorescence intensity in a defined, round region of interest (ROI, in this case a small spot) to determine the initial intensity value before the bleaching. This