similar to: Are there ANOVA for compositional data?

The compositional data xi=(x_i1,x_i2,...,x_in), for each fixed i , xij>0, and sum(xij)=1; I want to compare the mean( u_i) of several groups i.e. H0: u_1=u_2=...=u_N or H0: u_11=u_21=...=u_N1 Are there any ANOVA tpye tools to do this work in R? Thanks, WEN S Q [[alternative HTML version deleted]]

Dear People, I have a vexing problem related to pseudo-random number generation, and would appreciate any help and advice. This problem is not directly related to R, and the only reason I am posting it to this list is that my implementation is using R. Let me describe my problem by giving an example, that is close to what I am trying to do. Suppose we are given a stream of pseudo-random numbers,

Hey all - I am a newbie on mixed-effects models. I want to estimate the following model: Y_it = alpha_0t + alpha_1t*X_it + e_it alpha_0t = gamma_00 + u_0t alpha_1t = gamma_10 + gamma_11*W_it + u_1j Where Y is my outcome, X is my level-1 predictor, and W is my level 2 predictor. I am not sure if I am doing it right. Is this the correct specification of the formula? model = lmer(Y ~ X + X:Y + (

I discovered recently that the phrase "Oja median" produces no hits in Jonathan Baron's very valuable R search engine. I found this surprising since I've long regarded this idea as one of the more interesting notions in the multivariate robustness literature. To begin to remedy this oversight I wrote a bivariate version and then decided that writing a general p-variate version

Dear all, I have come across a very confusing matter regarding dist() supplied by the amap package: --- m is just a test matrix > library(amap) Loading required package: mva Warning message: package 'mva' has been merged into 'stats' > m a b c aa 0.1 0.2 0.3 bb 2.0 3.0 4.0 cc 2.0 4.0 6.0 dd 0.3 0.2 0.1 > ds<-dist(m,method="pearson") >

Dear R users, I am a mathematics postdoc at UC Berkeley. I have written a package in a Computational Algebra System named Singular http://www.singular.uni-kl.de to compute the Maximum Likelihood of a given probability distribution over several discrete random variables. This package gives exact answers to the problem. But more importantly, it gives All MLE solutions. My understanding is that

Hi, I want to get random number which is uniformly distributed on the unit disc. How can I do that with R? Best wishes, WAN WAN [[alternative HTML version deleted]]

Sorry for my fault that I just sent this mail to xen-devel first. Hi,everyone I am trying to use VHD disk format on guestVM. My environment is fedora12+ xen-4.0.1-rc6 + xen-stable/2.6.32.15+libvirt-0.8.3 I success on booting VM in VHD disk correctly using "xm create VM.cfg" and everything runs alright Where disk = [

Dear R users, i am wondering if anyone has some hints for this problem (i have not found a clear answer after searching the R-mailing list archive, 'help.search' in R, and R-Wiki, and the like...): let's assume that i have 4 periods compositional time series data: t=1, A=0.1; B=0.5; C=0.4 t=2, A=0.2; B=0.4; C=0.4 t=3, A=0.5; B=0.3; C=0.2 t=4, A=0.4; B=0.3;

again, under another subject: sorry, maybe an all too trivial question. But we have power data from J frequency spectra and to have the same range for the data of all our subjects, we just transformed them into % values, pseudo-code: power[i,j]=power[i,j]/sum(power[i,1:J]) of course, now we have a perfect linear relationship in our x design-matrix, since all power-values for each subject sum up

Suppose I've got p vectors x1, ..., xp each of length n merged in a (n x p)-matrix X. Now I want to calculate the number of distinct values of row-vectors (x_i1, ..., x_ip). Is there a way to do this in R? Thanks in advance, Jan. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Suppose I've got p vectors x1, ..., xp each of length n merged in a (n x p)-matrix X. Now I want to calculate the number of distinct values of row-vectors (x_i1, ..., x_ip). Is there a way to do this in R? Thanks in advance, Jan. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Hello, I'm wondering if anyone has benn able to successfully get g726-16 passthrouhg to work? I am wanting to use this codec instead of g729 as I'm running out of DSPs using a high complexity codec on the Ciscos. I would think it would work just as g729 does, which has been working fine for me, but it does not. G726-32 does work great however, but it's like Asterisk doesn't

Hi everyone! In my new CentOS 4.0 (final ;-) I installed apache and activates SELinux, When I tried to enter a page of some user (localhost/~user) said to me that she did not have permissions to see it, and this was incorrect, because the permissions was 755 (/home/user/html). After awhile it was happened to me to deactivate selinux and I reinitiated the machine and I work! Somebody knows why