similar to: Contrastes con el paquete survey (svycontrast)

Dear forum members, What is the formula to calculate denominator degrees of freedom (den df) for nonlinear mixed-effect models with covariates? My model is similar to a CO2 uptake example from Pinheiro and Bates (2000, page 376). In this CO2 dataset, there are two treatments and two types (84 observations in total), but den df for each parameter of the model is 64. Isn’t it too high? Your

Hi, I ran the example on pp. 799-800 from Machael Crawley's "The R Book" using package survival v. 2.36-5, R 2.13.0 and RStudio 0.94.83. The model is a Cox's Proportional Hazards model. The result was quite different compared to the R Book. I have compared my code to the code in the book but can not find any differences in the function call. My results are attached as well as a

Hi Rajiv, On Mon, Jul 20, 2009 at 4:01 PM, Rajiv Rajaian <rajiv.grid@gmail.com> wrote: > > Hi  deshantm > > This is Rajiv, from Anna university India. > sorry for disturbing you.. > > I ve installed xen-3.2.3-0xs  in Red Hat Enterprise Linux Server release 5.. > & I ve successfully created guest domain with Fedora 4  image file > While booting the guest domain

Dear all, I am reading the book "statistic analysis of gene expression microarray data" by Terry Speed . On page16-17, there is an example that I can not repeat. The content can be found using below linkage:

Hi, Dear R community, These is a paper of "Wheat Grain Yield Response to N Application Evaluated through Canopy Reflectance". This paper of "Materials and Methods" said "A single regression equation relating A and B was derived for each location using PROC NLIN (SAS Inst. 1990). Slopes and intercepts for the independent regressions were compared with an *approximate

Hello to all the people in the forum, I'm trying to perform a Scheirer-Ray-Hare test in R, and despite I've found some posts talking about this test, there no one in which I can found how to do so. Can anyone explain me how to do the test? Thanks Pablo -- View this message in context: http://r.789695.n4.nabble.com/Scheirer-Ray-Hare-tp3818476p3818476.html Sent from the R help mailing

Dear R experts, I am trying to analyze data from an article, the data looks like this Patient Age Sex Aura preCSM preFreq preIntensity postFreq postIntensity postOutcome 1 47 F A 4 6 9 2 8 SD 2 40 F A/N 5 8 9 0 0 E 3 49 M N 5 8 9 2 6 SD 4 40 F A 5 3 10 0 0 E 5 42 F N 5 4 9 0 0 E 6 35 F N 5 8 9 12 7 NR 7 38 F A 5 NA 10 2 9 SD 8 44 M A 4 4 10 0 0 E 9 47 M A 4 5 8 2 7 SD 10 53 F A 5 3 10 0 0 E 11

Hola a todos. Estoy utilizando la función lmer del paquete lme4 para ajustar un modelo mixto. Tengo varias variables en mi data.frame, unas son a nivel individual y otras a nivel de comarcas. Listo algunas. ingre_6 : Ingresos (nivel individual) iscs_a : un indicador sintético resumen de otras variables, calculado mediante componentes principales. sau_com : superficie agraria útil de cada una

Dear list colleagues, I'm trying to come up with a test question for undergraduates to illustrate comparison of means from a complex survey design. The data for the example looks roughly like this: mytest<-data.frame(harper=rnorm(500, mean=60, sd=1), party=sample(c("BQ", "NDP", "Conservative", "Liberal", "None", NA), size=500,

Hello all, This is something that I am sure has a really suave solution in R, but I can't quite figure out the best (or even a basic) way to do it. I have a simple linear regression that is fit with lm for which I would like to estimate the x intercept with some measure of error around it (confidence interval). In biology, there is the concept of a developmental zero - a temperature under

Using the survey package I find it is convenient and easy to get estimated proportions using svymean, and their corresponding estimated standard errors. But is there any elegant/simple way to calculate corresponding confidence intervals for those proportions? Of course +/- 1.96 s.e. is a reasonable approximation for a 95% CI, but (incorrectly) assumes symmetrical distribution for a proportion.

Dear R users, I know, this is the second time i return on this topic. Sorry, but this analysis is of great value for me, and i hope someone can help me. I need to model a time-varying effect in a Cox model. Briefly explained here: http://books.google.com/books?id=9kY4XRuUMUsC&lpg=PP1&hl=it&pg=PA147#v=onepage&q&f=false

Hello, svyvar from the survey package computes variances (with standard errors) from survey design objects. Is there any way to compute standard deviations and their standard errors in a similar manner? Thanks a lot, Sebastian

Hi all, I am running a set of logistic regressions, where we want to use some weights, and I am not sure whether what I am doing is reasonable or not. The dependent variable is turnout in an election - i.e. survey respondents were asked whether or not they voted. The percentage of those who say they voted is much higher than the actual turnout, probably due both to non-response bias and social

I'm trying to use the survey package to calculate a risk difference with confidence interval for binge drinking between sexes. Variables are X_RFBING2 (Yes, No) and SEX. Both are factors. I can get the group prevalences easily enough with result <- svyby(~X_RFBING2, ~SEX, la04.svy, svymean, na.rm = TRUE) and then extract components from the svyby object with SE() and coef() to do the

Hola a todos, Ajusté un modelo lineal mixto en el cual tengo una interacción significativa entre dos facotres, cada uno con 3 niveles (osea un diseño 3x3). Esta es la tabla de medias             P1                P 2             P3 M1      220.66       311.85        260.80 M2      348.57        89.33         191.71 M3      103.57      381.37         511.62 Con el fin de interpretar

Estimados usuaRios: deseo hacer una consulta, debo realizar una prueba no paramétrica a un conjunto de datos que posee una variable de respuesta con tres factores que poseen interacción, o al menos debo probar ello. Es decir, realizar una prueba Scheirer-Ray-Hare o de Contrastes de rangos alineados. Estuve trabajando el paquete ARTool aunque no me queda claro

Hello, folks! Currently me with one other guy are trying to play with clang. The proposal may seem stupid, excuse me, if it was already discussed, we just want to try to implement something useful which seems absent for now. Ok, the idea. It seems interesting to try to make lexer a little bit more efficient in terms of macro expanding by applying partial expansion of macros. the idea is that

I am using R (version 2.13.0; 2011-04-13). Platform: x86_64-pc-mingw32/x64 (64-bit). The messages that I see is “This application has requested the Runtime to terminate it in an unusual way. Please contact the application's support team for more information” I am using the following libraries: > library(package="gmodels") > library(package="gtools")

Using 10G interfaces on kernel 3.2 + xen 4.1.2 we''re seeing: dom0 ~9.2Gbps domU *~2.5Gbps* dmesg on domU: XENBUS: Device with no driver: device/vif/0 Is this normal? Thanks Kristoffer -- View this message in context: http://xen.1045712.n5.nabble.com/10gig-in-domU-tp5634876p5634876.html Sent from the Xen - User mailing list archive at Nabble.com.