similar to: SAS PROC NLMIXED into R

Hello friends of R list, Im a physician and Im not that good in statistics. I have posted similar email in the epi-sig list but no one aswered so far. Im cunducting a systematic review on two diagnostic test for a particular tropical disease. I found a free software that make most of the analysis called MetaDiSc. However there is a paticular analysis that I wuould like to try that it is not

To account for likely differences between families in naturalization rates, we fitted a generalized linear mixed model, using PROC NLMIXED in SAS10, with the naturalization rate per genus (that is, the number of naturalized species in a genus as a proportion of the total number of introduced species in a genus) as the response variable, a variable coding genera as containing at least one native

I'm trying to run a zero-inflated regression for a continuous response variable in R. I'm aware of a gamlss implementation, but I'd really like to try out this algorithm by Dale McLerran that is conceptually a bit more straightforward. Unfortunately, the code is in SAS and I'm not sure how to re-write it for something like nlme (if at all possible - with conditions etc). Does

Hi, I am trying to convert the following NLMIXED code to NLME, but am running into problems concerning 'Singularity in backsolve'. As I am new to R/S-Plus, I thought I may be missing something in the NLME code. NLMIXED *********** proc nlmixed data=kidney.kidney; parms delta=0.03 gamma=1.1 b1=-0.003 b2=-1.2 b3=0.09 b4=0.35 b5=-1.43 varu=0.5; eta=b1*age+b2*sex+b3*gn+b4*an+b5*pkn+u;

Hello! Can anyone help me. I am using the posterior.mode from the result of glmmML. It apears to be different from the BLUe estimate of the RANDOM statement in PROC NLMIXED in SAS. Why is that? Thank you Ronen [[alternative HTML version deleted]]

Dear All, reading an article by Rodolphe Thiebaut and Helene Jacqmin-Gadda ("Mixed models for longitudinal left-censored repeated measures") I have found this program in SAS proc nlmixed data=TEST QTOL=1E-6; parms sigsq1=0.44 ro=0.09 sigsq2=0.07 sigsqe=0.18 alpha=3.08 beta=0.43; bounds $B!](B1< ro < 1, sigsq1 sigsq2 sigsqe >= 0; pi=2*arsin(1); mu=alpha+beta*TIME+a i+b i*TIME;

Hello, I am writing a meta-analysis on the complication and reoperation rates after 5 treatment modalities of a distal radius fracture. I have a code to compare the complication and reoperation rates. Currently it is using a mixed effects model. Is it possible to change the code so a random effects model is used? Thank you very much, Jay R code library(meta) library(readxl) All <-

The code in your mail in a mangled mess, since you posted in HTML. Please configure your email client to send emails in plain text. Could you explain what exactly you mean by "Currently it is using a mixed effects model. Is it possible to change the code so a random effects model is used?" Best, Wolfgang >-----Original Message----- >From: R-help [mailto:r-help-bounces at

Hello R masters, I have sent this same message to other lists and none so far could give some light. I was trying to use odfWeave to generate a report from R and Im getting an error that I think is related to latin characters. I looked around and did find some stuff related to this problem about Sweave http://labmoluscos.wordpress.com/2010/02/18/sweave-latex-character-encoding/ but did not find a

?s 23:47 de 31/01/2024, Roy Mendelssohn - NOAA Federal via R-help escreveu: > HI All: > > We are trying to figure out a problem that is occurring with a package, and we need a non-NOAA person with a Windows computer with the latest R to test for us what is failing (but works on Macs and Linux from different sites). Part of the problem is there is an Akamai service in between that we

Hi, I have two tables with identical fields(assignments, and matches). However, only one of them accept NULL values for its fields. I want to be able to find out all the matches for any assignment. The search should retrun the match if the corresponding filed is NULL as well. I wrote down the following code but it doesn''t return anything. I wonder what I am doing wrong. Any help is

Hi everyone, I connected a data in FileMaker server with odbcConnect. When I call the data "CANCERS" using sqlFetch, it looks okay. However, the number of obs was different with the actual number. If I read the same data from Microsoft Access, there are 656 obs. in it but srt() shows me 600 obs. Does anyone know why this happened and how to read all obs. from R? Thank you in advance.

HI All: We are trying to figure out a problem that is occurring with a package, and we need a non-NOAA person with a Windows computer with the latest R to test for us what is failing (but works on Macs and Linux from different sites). Part of the problem is there is an Akamai service in between that we feel may be messing up things. To do this you would need to be willing to: 1. install the

Hi everybody: I?m trying to rewrite some routines originally written for SAS?s PROC NLMIXED into LME4's glmer. These examples came from a paper by Nelson et al. (Use of the Probability Integral Transformation to Fit Nonlinear Mixed-Models with Nonnormal Random Effects - 2006). Firstly the authors fit a Poisson model with canonical link and a single normal random effect bi ~ N(0;Sigma^2).The

Do Users of Nonlinear Mixed Effects Models Know Whether Their Software Really Works? Lesaffre et. al. (Appl. Statist. (2001) 50, Part3, pp 325-335) analyzed some simple clinical trials data using a logistic random effects model. Several packages and methods MIXOR, SAS NLMIXED were employed. They reported obtaining very different parameter estimates and P

My thanks to Drs. Armstrong, Bates, Harrell, Liaw, Lumley, Prager, Schwartz, and Mr. Wang for their replies. I have pasted my original message and their replies below. After viewing http://www.itl.nist.gov/div898/strd/ as suggested by Dr. Schwartz, it occurred to me that it might be educational to search for some data repositories on google. I was able to find some,though I'm sure many of

Hi there, I am using SAS Proc NLMIXED to maximize a likelihood with multivariate normal random effects. An example is the two part random effects model for repeated measures semi-continous data with a cluster at 0. I use the "model y ~ general(loglike)" statement in Proc NLMIXED, so I can specify a general log likelihood function constructed by SAS programming statements. Then the

hi Jay, Consult a local statistician. Statistics is not you think is (namely simple computations, R and probably plotting..). regards, vito Jay Zola <jayjay.1988 at hotmail.nl> ha scritto: > Hello, > > > I am medical student, writing a meta-analysis on complication and > reoperation rates after the five most common treatments of distal > radius fractures. I have

I'd like to estimate an ordinal logistic regression with a random effect for a grouping variable. I do not find a pre-packaged algorithm for this. I've found methods glmmML (package: glmmML) and lmer (package: lme4) both work fine with dichotomous dependent variables. I'd like a model similar to polr (package: MASS) or lrm (package: Design) that allows random effects. I was

Hello, I am medical student, writing a meta-analysis on complication and reoperation rates after the five most common treatments of distal radius fractures. I have been busy with the statistics for months by my self, but find it quite hard since our classes were very basic. Now I want to compare the treatment modalities to see if there are significant differences. Using R I was able to