similar to: parameterization of factor in R

Dear all, I have a problem when I was working on Age-Period-Cohort study in R. I tried to put constraint to two coefficients on age (so that to solve the identification problem due to linear dependency). But I don't know how to do this in R(put constraint). If you could give me some suggestion, it will be very helpful! Oana

Hi, I am trying to parameterize the following mixed model (following Piepho and Ogutu 2002), to test for a trend over time, using multiple sites: y[ij]=mu+b[j]+a[i]+w[j]*(beta +t[i])+c[ij] where: y[ij]= a response variable at site i and year j mu = fixed intercept Beta=fixed slope w[j]=constant representing the jth year (covariate) b[j]=random effect of jth year, iid N(0,sigma2[b]) a[i]=random

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

Hello Everyone Both the MCMCpack and the bayesm libraries allow us to make draws from the Inverse Wishart distribution. But I wanted to find out how exactly is the Inverse Wishart distribution parameterized in these libraries. The reason I ask is the following: Now its generally standard to express Inverse Wishart as IW(0.5 * DOF,0.5* Scale). (DOF-> Degree of freedom, Scale -> Scale

Hello, I would like to perform a generalized singular value decomposition with R. The only possibility I found is "GSVD" that is based on LAPACK/BLAS. Are there other possibilities too? If not, has anybody used LAPACK/BLAS under Windows XP? How can I install them? Following [1] did not help. I hope this is the right place for my question. Thank you very much! Oana Tomescu [1]

Hello, I want to solve some two dimensional equation system with R. Some systems are not solvable analytically. Here is an example: (I) 1/n*sum{from_i=1_to_n}(Xi) = ln lambda + digamma(c) (II) mean(X) = x / lambda I want to find lambda and c, which R-function could do that task? Carsten [[alternative HTML version deleted]]

I use the lme function from the nlme library (or alternatively from the Matrix library) to estimate a random effects model. Both functions return the covariance matrix of the estimated parameters. I have the following question: Is it possible to retrieve the information matrix of such a model (ie from the fitted object)? In particular, the information matrix can be computed as a sum of individual

Hello all (especially MCLUS users). I'm trying to make use of the MCLUST package by C. Fraley and A. Raftery. My problem is trying to figure out how the (model) identifier (e.g, EII, VII, VVI, etc.) relates to the covariance matrix. The parameterization of the covariance matrix makes use of the method of decomposition in Banfield and Rraftery (1993) and Fraley and Raftery (2002) where

Dear R community, I am new to R, a reforming SAS user :) I am running R 2.10.1 on a Windows XP machine. I would like to write linear functions of my coefficient parameter estimates from a glm, but am having a difficult time understanding the parameterization R uses. In the toy example below I am running a glm on binomial data, with clones and lines within clones as fixed effects, each with 6

Dear R users, When fitting a lme() object (from the nlme library), is it possible to test interactions *before* main effects? As I understand, R conventionally re-orders all terms such that highest-order interactions come last - but I??d like to know if it??s possible (and sensible) to change this ordering of terms. I??ve tried the terms() command (from aov) but I don??t know if something

Greetings everyone, We would like to improve LNT. The following RFC describes two LNT enhancements: * Custom (extensible) metrics * Test parameterization The main idea is in document https://docs.google.com/document/d/1zWWfu_iBQhFaHo73mhqqcL6Z82thHNAoCxaY7BveSf4/edit?usp=sharing. Thanks, Elena. -------------- next part -------------- An HTML attachment was scrubbed... URL:

Dear R-helpers, I try to perform glm's with negative binomial distributed data. So I use the MASS library and the commands: model_1 = glm.nb(response ~ y1 + y2 + ...+ yi, data = data.frame) and predict(model_1, newdata = data.frame) So far, I think everything should be ok. But when I want to perform a glm with a subset of the data, I run into an error message as soon as I want to predict

Hi Elena, > First of all, now LNT uses make-style of running tests and parse results from result csv file. Are there any plans to go to cmake? There are two test drivers in LNT. "lnt runtests nt" uses the old Makefile-based system and "lnt runtests test-suite" uses the new cmake-based system. It sounds like you are using the former and should switch to the latter.

> On May 12, 2016, at 11:21 PM, Elena Lepilkina via llvm-dev <llvm-dev at lists.llvm.org> wrote: > > Hi all, > > As we understood great changes will be done in LNT, so we are waiting to new LNT version and stopped our work in LNT. > > One more question about using test-suite separately with cmake. Cmake can only build all tests and generate lit tests. After that we

I got the following warning messages when I did a binomial logit regression using glm(): Warning messages: 1: Algorithm did not converge in: glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, 2: fitted probabilities numerically 0 or 1 occurred in: glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, Can some one share your thoughts on how to

> On May 25, 2016, at 1:54 AM, Elena Lepilkina via llvm-dev <llvm-dev at lists.llvm.org> wrote: > > Hi Matthias, > > Thank you for your answer. > But can you answer for some more questions? > First of all, now LNT uses make-style of running tests and parse results from result csv file. Are there any plans to go to cmake? As James already said "lnt runtest

Shalini Raghavan 3M Pharmaceuticals Research Building 270-03-A-10, 3M Center St. Paul, MN 55144 E-mail: sraghavan at mmm.com Tel: 651-736-2575 Fax: 651-733-5096 ----- Forwarded by Shalini Raghavan/US-Corporate/3M/US on 08/16/2004 11:25 AM ----- Shalini

Dear R users; I'm looking for some hint about how to deal with the following situation: Response = Y Factor A = levels: 0, 1 Factor B = levels: 0, 1 Factor C = levels: 1,2,3,4 Model: Logistic 3-parms. where th1~1+A+C, th2~1+C; th3~1 For 'simplicity' (for me) I'm using the SAS contrast parameterization. The output looks like Beta p-value th1.(Intercept) 550

Hi Matthias, Thank you for your answer. But can you answer for some more questions? First of all, now LNT uses make-style of running tests and parse results from result csv file. Are there any plans to go to cmake? As I understood lit will run and collect all metrics, but there is no opportunity to make any settings for choosing what metrics I would like to collect. Test reports files allow to

Hi, anyone knows about any functions in R can get multidimensional integration not over a multidimensional rectangle (not adapt). For example, I tried the following function f(x,n)=x^n/n! phi.fun<-function(x,n) { if (n==1) { x }else{ integrate(phi.fun, lower=0, upper=x, n=n-1)$value } } I could get f(4,2)=4^2/2!=8, but failed in f(4,3)=4^3/3! Thanks Best, Lynette