[R] Plackett-Dale Model in R

Dear R users, Can someone inform me about a library/function in R that fits a Plackett-Dale model ? Thanks in advance Pryseley - [[alternative HTML version deleted]] __

[R] lme convergence

Dear R-Users, Is it possible to get the covariance matrix from an lme model that did not converge ? I am doing a simulation which entails fitting linear mixed models, using a for loop. Within each loop, i generate a new data set and analyze it using a mixed model. The loop

[R] Help: lme

Good day R-Users, I have a problem accessing some values in the output from the summary of an lme fit. The structure of my data is as shown below (I have attached a copy of the full data). id trials endp Z.sas ST 1 1 -1 -142.42884 1 1

Re: [R] Help: lme

suggestions from Doug Bates and myself. Is there some reason the post by Andrew Robinson is not working? Harold -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Pryseley Assam Sent: Thursday, June 01, 2006 6:25 AM To: r-help@stat.math.ethz.ch Subject: [R

[R] Query: lme output

Dear R-Users I have a problem accessing some values in the output from the summary of an lme fit. I fit the model below: ggg - lme (ST~ -1 + as.factor(endp):Z.sas + as.factor(endp), data=dat4a, random=~-1 + as.factor(endp) + as.factor(endp):Z.sas|as.factor(trials),

[R] query: lme

Good R-users, I have difficulties accessing the variance components for an lme fit when the variance covariance matrix of the random effects is not positive definite. For example, i fit the following model: ggg - lme (ST~ -1 + as.factor(endp):Z.sas + as.factor(endp),

[R] query: lme

Dear R Users I have difficulties accessing the variance components for an lme fit when the variance covariance matrix of the random effects is not positive definite. Can anyone inform me on how to get by this ? Thanks in advance Pryseley

[R] Pairewise Likelihood

Dear R-users Can anyone inform me of a library or more specifically functions that can maximise (or calculate) a Pairwsie likelihood from a data. Better still, i would like to know if there is a function (library) that fits regression models based on pairwise likelihoods.

[R] Log Cholesky parametrization in lme

Dear R-Users I used the nlme library to fit a linear mixed model (lme). The random effect standard errors and correlation reported are based on a Log-Cholesky parametrization. Can anyone tell me how to get the Covariance matrix of the random effects, given the above mentioned parameters

[R] Three level linear mixed models

Hello R-users Is it possible to fit a three level linear mixed effect model in R? If anyone has an idea or sample code, i will appreciate it very much if i can receive it. I am reading the book by Pinheiro and Bates but have not come across that yet! Kind regards

[R] lme and gls : accessing values from correlation structure and variance functions

Dear R-users I am relatively new to R, i hope my many novice questions are welcome. I have problems accessing some objects (specifically the random effects, correlation structure and variance function) from an object of class gls and lme. I used the following models:

[R] Help with lme and correlated residuals

Dear R - Users I have some problems fitting a linear mixed effects model using the lme function (nlme library). A sample data is as shown at the bottom of this mail. I fit my linear mixed model using the following R code: bmr -lme (outcome~ -1 + as.factor(endpoint)+

[R] Linear mixed model with correlated residual

Dear R- users Does any one know how to fit a linear mixed model such that the residuals ( grouped by a variable say gender) are correlated and have a covariance matrix (in this case a 2 by 2 covariance matrix). Thanks in advance Pryseley

[R] Help with lme code

Dear R - Users I have some problems fitting a linear mixed effects model using the lme function (from the library nlme). A sample data is as shown at the bottom of this mail. I fit my linear mixed model using the following R code: bmr -lme (outcome~ -1 + as.factor(endpoint)+

[R] Help with classes and generic function

Hello R-Users I am pretty new to R and forgive me if my questions are childish , nevertheless i need help. I have some problems while trying to grasp the idea of a generic function. I understand that i have to build methods with naming syntax function.class, where class is the

[R] Help with functions

Dear R-users I intend to create a function which calls some smaller other functions in return. Some of these smaller functions all call some functions. I do not know a good way to do this. I tried using the source() function to include the smaller functions within the main functions

Re: [R] Help with R: functions

Dear R-users Thanks for all the help/suggestions, they have been most helpful. Best regards Pryseley - [[alternative HTML version deleted]] __

[R] Help with R: functions

Hello R-users I am new to R and trying to write some functions. I have problems writing functions that takes a data set as an arguement and uses variables in the data. I illustrate my problem with a small example below: sample data #-- visual24-rnorm(30,3,5)

[R] Help with mixed effects models

Dear R-users I have problems using lme The model i want to fit can be viewed as a two-level bivariate model Two-level bivariate: bivariate (S coded as -1,T coded as 1) endpoint within trial OR It can equivalently be considered as a three-level model.Three-level: endpoint