search for: bennfine

I am trying to fit a linear mixed model of the form y_ij = X_ij \beta + delta_i + e_ij where e_ij ~N(0,s^2_ij) with s_ij known and delta_i~N(0,tau^2) I looked at the ecme routine in package:pan, but this routine does not allow for different Vi (variance covariance matrix of the e_i vector) matrices for each cluster. Is there an easy way to fit this model in R or should I bite the bullet and

I am using the Design library and robcov to compute variance-covariance matrices for clustered data regression. Is there an easy way to compute the F-test (i.e. linear hypothesis) for clustered data regression ? Thanks in advance! Benn

Hello: First, RExcel completely rocks! Hat tip to the development team. However, I don't understand something very basic about it. Suppose I put the following code in two cells and choose the Run Code option plot(runif(20),runif(20)) text(.5,.5,"Hello") I get the expected plot. HOWEVER, if I put the following in a cell and hit Run Code, I get nothing. runif(20) HOWEVER, If I

I want to fit a logistic regression model under a specified null hypothesis, say Ho:Beta_k=1 Is there a way to constrain a parameter in glm.fit ?? Thanks in advance! --------------------------------- [[alternative HTML version deleted]]

How does one specify the weights for WLS in the systemfit command ? That is, there is a weight option in lm(), but there doesn't seem to be weight option for systemfit("WLS") Thanks!

http://www.fischeredv.de/zxripc/3bg8jc&mc9en1iqo46=lgyc8dv37ec694oxxqf2enjlgwsr70 [[alternative HTML version deleted]]

I have ported some R code to C to make it faster. I can perform .Call("foobar",....) once and it works fine. Absolutely correct answer. If I put a loop inside foobar and run the main code routine more than 100 times, it crashes R. Or if I call .Call("foobar"....) seperately more than two tims it crashes R. For the most part I am doing matirx multiplies using EXP