similar to: Computing row means for sets of 2 columns

Dear everyone, I am doing some work about species richness estimation. Nonparametric estimation (such as Chao1, Jacknife1) can be done just using function "specpool()" and "estimateR()" in package "vegan". The problem is that I can not found any functions for parametric estimation (such as MMMeans, MMruns, Michaelis-Menten). Do you know any function for doing this?

The result of rowsum() in R doesn't have the dimnames I'd expect, e.g.: > rowsum(matrix(1:12, 3,4), c("Y","X","Y")) [,1] [,2] [,3] [,4] 1 2 5 8 11 2 4 10 16 22 whereas S-Plus gives the more useful result: [,1] [,2] [,3] [,4] X 2 5 8 11 Y 4 10 16 22 This is because R's rowsum() code gives the

Hi everyone, I have a question about functions. I have two functions: Sampling=function(fname,Total,nSample,nPatient){......return(list(Gmean,Gsd))} Power=function(alfa,m1,m2,s1,s2,n1,n2){....return(powe)} I want to use "Gmean" and "Gsd" which are the returned values from "Sampling" function, in "Power" function. i.e. when I call "power" I want

Dear all, I am using "coxme" function in Kinship library to fit random treatment effect nested within centre. I got 3 treatments (0,1,2) and 3 centres. I used following commands, but got an error. > ugroup=paste(rep(1:3,each=3),rep(0:2,3),sep='/') > mat1=bdsmatrix(rep(c(1,1,1,1,1,1,1,1,1),3),blocksize=rep(3,3),dimnames=list(ugroup,ugroup)) >

Hi, The result of rowsum() in R doesn't have the dimnames I'd expect, e.g.: > rowsum(matrix(1:12, 3,4), c("Y","X","Y")) [,1] [,2] [,3] [,4] 1 2 5 8 11 2 4 10 16 22 whereas S-Plus gives the more useful result: [,1] [,2] [,3] [,4] X 2 5 8 11 Y 4 10 16 22 This is because R's rowsum() code gives

Hello I'm having trouble being able to plot more than one set of data using the plotmeans function. Plotmeans is where you can plot means and standard errors easily. Here is part of my code plotmeans(Mmeans ~ Month, data=BW.hydat.05HG001, n.label=FALSE)#one set of data par(new=TRUE) plotmeans(Mmeans ~ Month, data=Mouth.hydat.05HH001, col="red")#the second set of data Basically

I'm using the following model to do an analysis faicout <- glm(cbind(events,patnums-events) ~ as.factor(treat) + as.factor(numtrial), family = binomial ) Is this example there are 4 treatments . In the glm object I can find the contrasts of the main treats vs the first i.e. 2v1, 3v1 and 4v1 ... however I would like to get the complete set including 3v2, 4v2, and 4v3 ... along with the

Dear R users: I'm new to R and am trying to fit a mixed model Cox regression model with coxme function. I have one two-level factor (treat) and one covariate (covar) and 32 different groups (centers). I'd like to fit a random coefficients model, with treat and covar as fixed factors and a random intercept, random treat effect and random covar slope per center. I haver a couple of

Hello all, I am a 2 month newbie to R and am stumped. I have a data set that I've run multivariate stats on using the manova function (I included the data set). Now it comes time for a table of effect sizes with significance. The univariate tests are easy. Where I run into trouble filling in the table of effect sizes is the Mahalanobis D as an effect size. I've included the table so

Hello, I have many subsets of x. I want to get the standard deviation for each subset with the same factor levels. For "FUN=mean" and "FUN=median" I am using "ave()". Can anybody tell me the match of "ave()" for using "FUN=SD"? At the beginning I used aggregate(), also for mean and median. But aggregegate make ""arithmetical

Hi! I am studying HTB. I used the topology showed below: 10.1.0.1 ------ 10.1.14.25 (eth1) / 192.168.10.10 (eth0) ----- 192.168.10.11 (host 1) (router) (host 2) All machines use Debian Etch. Has a Iptables masquerading rule to eth1 in the router machine . The NICs are 100 Mb/s. The host 1 has Apache 2 and a file with 670 MB (CD ISO image) to download. When I used

Already outdated, but rounded ;) I literally just opened the 5yrs logo, changed the text and then done: inkscape -z -o logo/fish-10yrs.{png,svg} cp {logo,website}/fish-10yrs.svg and then updated the rest of the files. Signed-off-by: Martin Kletzander <mkletzan@redhat.com> --- Makefile.am | 2 + logo/fish-10yrs.png | Bin 0 -> 65790 bytes logo/fish-10yrs.svg |