similar to: split-plot design?

Dear all, I am trying to use the Matrix package to do some calculations on rather large and sparse matrices. An example of such a matrix is given below. mig<-0.2 side<-10 np<-side^2 mig.mat<-matrix(0,np,np) diag(mig.mat[1:(np-side),(side+1):np])<-mig/4 diag(mig.mat[(side+1):np,1:(np-side)])<-mig/4 diag(mig.mat[-np,-1])<-mig/4

Hi all, I am wondering whether there is a small bug in the binom.test function of the ctest library (I'm using R 1.6.0 on windows 2000, but Splus 2000 seems to have the same behaviour). Or perhaps I've misunderstood something. the command binom.test(11,100,p=0.1) and binom.test(9,100,p=0.1) give different p-values (see below). As 9 and 11 are equidistant from 10, the mean of the

hi! I apologize in advance if this is a newbie dumm question, but I really can't figure it ou. I have lists of sumeric and character data on some URLs, which look like this: <photo id="5876248819" owner="13716719 at N04" secret="faf9bb7f52" server="5264" farm="6" title="our rose garden" ispublic="1"

Dear All, I am quite new to R and am having a problem trying to run a linear model with random effects/ a regression- with particular regard to my variable lengths being different and the models refusing to compute any further. The codes I have been using are as follows: vc<-read.table("P:\\R\\Testvcomp10.txt",header=T) >> attach(vc) > > family<-factor(family) >

Hi ! We study the effect of several variables on fruit set for 44 individuals (plants). For each individual, we have the number of fruits, the number of flowers and a value for each variable. Here is our first model in R : y <- cbind(indnbfruits,indnbflowers); model1 <-glm(y~red*yellow+I(red^2)+I(yellow^2)+densite8+I(densite8^2)+freq8_4+I (freq8_4^2), quasibinomial); - We have

Dear all, I'm puzzled by the definition of r-squared used by 'lm' when the model is without intercept. The help for summary says: >r.squared: R^2, the ``fraction of variance explained by the model'', > > R^2 = 1 - Sum(R[i]^2) / Sum((y[i]- y*)^2), > > where y* is the mean of y[i] if there is an intercept and > zero

Hi, I am new to this list as a poster, but a reader for some time. I've using R for several weeks now, and I have a lot of questions about certain procedures. Here I go: I want to test if there are differences in the time spent by pollinators visiting flowers of a given plant species, according to a number of experimental manipulations made on those flowers. All experimental

I'm sure there must be various peak-finding algorithms out there. Not knowing of any, I have written one myself*, but I thought I'd ask to see what's out there. Basically, I have a 2-dimensional data set and I want to identify local peaks in the data, while ignoring "trivial" peaks. My naive algorithm first identifies every peak and valley (point of inflection change

Sorry to all that are angry about the form of my previous mail. I didn't realise what would happen :((. Here it is in (hopefully) plain text (if my mailer doesn't spoil it again): ############## Dear developers, I have a problem with some discrepancy between R 1.2.1 for Windows and R 1.2.2 (and less) for Linux. While trying to correct the wilcox.test (see my previous bug report) I

R users, How come R doesn't have a mascot? Linux has one and so does LaTeX, so shouldn't R? I personally think that associating a "friendly face" with R would be a good thing for R (one letter names can be quite intimidating). I apologize if this is addressed in the FAQ. I searched the FAQ as well as the mailing list archives and checked ?mascot but to no avail. ;-)

Hello, After putting together interaction code that worked for a single pair of interactions, when I try to evaluate two pairs of interactions( flowers*gopher, flowers*rockiness) my computer runs out of memory, and the larger desktop I use just doesn't go anywhere after about 20 minutes. Is it really that big a calculation? to start: mle2(minuslogl = Lily_sum$seedlings ~ dnbinom(mu = a,

I buy flowers at a local market on a fairly regular basis. The flower vendors post their prices and if I want to buy only one or two flowers I will generally get the posted price. From time to time I want to buy large quantities of flowers, and sometimes a vendor will give me a better price than their posted price for the bulk order, but more often I have to offer them a higher price than the

Dear all, can anybody help me with the program below? The function predict.lda seems to be defined but cannot be used by errortest. The R version is 1.7.1 Thanks in advance, Stefan ---------------- library("MASS"); library("ipred"); data(iris3); tr <- sample(1:50, 25); train <- rbind(iris3[tr,,1], iris3[tr,,2], iris3[tr,,3]); test <- rbind(iris3[-tr,,1],

Your spf record is broken: dovecot.org. 39942 IN TXT "v=spf1 a -all" -- Jim Flowers <jflowers at ezo.net> -- This message has been scanned for viruses and dangerous content by MailScanner, and is believed to be clean.

Hello, I have difficulties to deal with multilevel model. My dataset is composed of 10910 observations, 1237 plants nested within 17 stations. The data set is not balanced. Response variable is binary and repeated. I tried to fit this model model<- glmmPQL( y ~ z1.lon*lun + z2.lat*lun + z1.lon*lar + z2.lat*lar + z1.lon*sca + z2.lat*sca +z1.lon*eta + z2.lat*eta, random = ~ lun + lar + sca

How come I allways immediately run into difficult stuff when I''m trying some new programming language? Am I blind for the simplicity of Ruby/Rail, which a see must be there? Anyway, don''t try and answer this rhetoric question. I have got another one for you, seemingly difficult. I''m struggling with a n:m relationship (in a database, that is) and its join table.

When I run ogginfo on a .ogg file I created with oggenc (RC3), it lists some things that I have questions about: 1. Ogginfo lists bitrate_average and bitrate_nominal. Now, bitrate_average (I assume) would mean the overall bitrate average for the entire file but what does bitrate_nominal mean? 2. How is the serial number created? Can I tell something specific about an .ogg file by looking at it?

On Thu, 13 Oct 2005, Misha Brukman wrote: >> How does this tentative plan sound: we have two more weeks of >> development, then start the release processing part on about Oct 31 >> (spooky!). > > In that case, it needs a cool "release name". :) Suggestions welcome! -Chris -- http://nondot.org/sabre/ http://llvm.org/

>> In that case, it needs a cool "release name". :) > > Suggestions welcome! Names of flowers are numerous but not exactly cool. I do not think any other software project has used flower names. Aaron

Hello, It's not easy to express clearly what I have in mind by stochastic growth models. I've been working with a couple of titles -Stochastic Models in Biology, for one, but both date back a decade or more. I'll try to illustrate the model a little more. A reasonably comparable situation is Bolker's analysis of how many flowers are created on a given plant among a patch of