similar to: Exact Maximum Likelihood Package

Hello, Great to see the new metafor package for meta-analysis. I would like to perform a meta-analysis in which I initially calculate the intercept of the model with a nested random-effects structure. In lme, this would be model<- lme(v3~1, random=~1|species/study, weights = varFixed(~Wt), method = "REML") where multiple effects sizes are measured for some studies and more than

Hi, I noted a discrepancy between R and openepi when I ran a fisher test with the same matrix. In R: > a=matrix(c(1,2,6,17), nrow=2) > a [,1] [,2] [1,] 1 6 [2,] 2 17 > fisher.test(a, conf.int=T) Fisher's Exact Test for Count Data data: a p-value = 1 alternative hypothesis: true odds ratio is not equal to 1 95 percent confidence interval: 0.02061498

I referred in my reply to a paper by Diaconis and Sturmfels. The exact reference is: Diaconis and Sturmfels, Algebraic algorithms for sampling from conditional distributions, Ann. Stat 26 (1998) 363-397. They cite the following: Besag and Clifford, Generalized Monte Carlo significance tests, Biometrika 76 (1989) 633-42. which actually contains your problem (section 3, Testing the Rasch model)

Dear users, I am trying to show the equation (including coefficients from the model estimates) for a gam model but do not understand how to. Slide 7 from one of the authors presentations (gam-theory.pdf URL: http://people.bath.ac.uk/sw283/mgcv/) shows a general equation log{E(yi )} = ?+ ?xi + f (zi ) . What I would like to do is put my model coefficients and present the equation used. I am an

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,

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]]

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 Hj0: u_1j=u_2j=...=u_Nj Are there any ANOVA tpye tools to do this work in R? Thanks, WEN S Q [[alternative HTML version deleted]]

My name is Michela Marignani and I'm an ecologist trying to solve a problem linked to knight' s tour algorithm. I need a program to create random matrices with presence/absence (i.e. 1,0 values), with defined colums and rows sums, to create null models for statistical comparison of species distribution phenomena. I've seen on the web many solutions of the problem, but none provides

Hi, all, Does anybody got the most general expression of a unitary matrix? I found one in the book, four entries of the matrix are: (cos\theta) exp(j\alpha); -(sin\theta)exp(j(\alpha-\Omega)); (sin\theta)exp(j(\beta+\Omega)); (cos\theta) exp(j\beta); where "j" is for complex. However, since for any two unitary matrices, their product should also be a unitary matrix. When I

Dear R-experts, I have a question on the formulas used in the gam function of the mgcv package. I am trying to understand the relationships between: y~s(x1)+s(x2)+s(x3)+s(x4) and y~s(x1,x2,x3,x4) Does the latter contain the former? what about the smoothers of all interaction terms? I have (tried to) read the manual pages of gam, formula.gam, smooth.terms, linear.functional.terms but

Dear R-experts, I have a question on the formulas used in the gam function of the mgcv package. I am trying to understand the relationships between: y~s(x1)+s(x2)+s(x3)+s(x4) and y~s(x1,x2,x3,x4) Does the latter contain the former? what about the smoothers of all interaction terms? I have (tried to) read the manual pages of gam, formula.gam, smooth.terms, linear.functional.terms but

Hello list. I'm trying to rebuild the 2.6.9.67.0.20.EL kernel, but it fails even without modifications. How did I try it? Created a (non-root) build environment (not a mock ) Installed the kernel.scr.rpm and did a rpmbuild -ba --target=`uname -m` kernel-2.6.spec 2> prep-err.log | tee prep-out.log The build failed at the end: Processing files: kernel-xenU-devel-2.6.9-67.0.20.EL Checking