similar to: Using factors in R

This is only partly a question about R, as I am not quite sure about the underlying statistical theory either. I have fitted a non-linear mixed-effects model with nlme. In the fixed part of the model I have a factor with three levels as explanatory variable. I would like to use Tukey HSD or a similar test to test for differences between these three levels. I have two grouping factors:

Could someone explain the meaming of this error message from augPred: > augPred(area3.pen.nlme, primary=~day) Error in predict.nlme(object, value[1:(nrow(value)/nL), , drop = FALSE], : Levels 1,2,3 not allowed for block > predict.nlme(area3.pen.nlme) does not produce an error. area3.pen.nlme was created with: > area3.pen.nlme <- nlme(area ~ SSlogis(day, Asym, xmid, scal),

Y aún así veremos si puedes apañarte... R no es muy buen gestor de memoria. Además de lo que te ha dicho Carlos, lo mejor es "programar bien", sabiendo cuando R hace copias de los objetos (cosa que hace muy a menudo), para evitarlo al máximo. Otra cosa que puede hacerse es "adelgazar los objetos" si no necesitas mucha de la información que contienen. Por ejemplo, si haces un

Dear all, I have tried to calculate a GLMM fit with lmer (lme4) and glmmPQL (MASS), I also used glm for comparison. I am getting very different results from different functions, and I suspect that the problem is with our dataset rather than the functions, but I would appreciate help in deciding whether my suspicions are right. If indeed we are attempting the wrong type of analysis, some

I am learning R, and although I have looked in the documentation, I may be asking something obvious. Sorry, if that is the case. In a split-plot design if there is a missing subunit summary gives me a table with two rows for the same factor, one in the error within section and one in the section using error between units. With no data missing the table is "normal". How does one

On Wed, 20 Sep 2006, Alexandre Sobral wrote: > Willian, > > Tnks It´s ok with unset command. > Fantastic! > I´m trying implement the Samba for AIX with Active Directory Integration. > Do you have any idea where can i find a How to about this. > I´ve already search from internet but I haven´t found yet. > > I know that Johann Kappacher in Vienna, Austria just

Vaya, polémica. ¿R puro es R base? Porque van ya por 6618 paquetes en CRAN http://cran.r-project.org/web/packages/index.html Y nadie sabe bien cuántos hay en github. Por cierto, ¿alguno conoce un código para hacer una estimación? El 06/05/2015 a las 16:09, Freddy Omar López Quintero escribió: 2015-05-06 11:00 GMT-03:00 Pedro Concejero Cerezo <pedro.concejerocerezo en

Y con fread.. directamente sobre el ".rel"... > datIn <- fread("DRUG-AE.rel") > dim(datIn) [1] 6821 8 Saludos, Carlos Ortega www.qualityexcellence.es El 29 de septiembre de 2015, 18:45, Pedro Concejero Cerezo < pedro.concejerocerezo en telefonica.com> escribió: > No tiene cabecera. Tendrás que ponerle nombre a las variables > Prueba este código,

Hola de nuevo, el archivo .rel esta en una carpeta comprimida .zip, yo la descomprimo con win rar. Como no sea ese el motivo por el que no me da el fichero completo? El 30 de septiembre de 2015, 10:30, MªLuz Morales <mlzmrls en gmail.com> escribió: > Pues soy gafe entonces, > no lo entiendo: > > > datIn <- fread("C:\\Users\\iphealthMariluz\\Documents\\Proyecto

Hola, Las presentaciones de la sesión sobre ecuaciones estructurales (SEM) así como los vídeos -como ya anunció Carlos ayer- están subidos a la página del Grupo de Madrid: http://r-es.org/Grupo+de+Inter%C3%A9s+Local+de+Madrid+-+GIL+Madrid Saludos! -- Pedro Concejero BigData - Analytics @ Telefónica I+D<http://www.tid.es/es/Paginas/default.aspx> Part of Telefónica

Para los perezosos como yo, para el talleR de mapas hoy en Madrid: install.packages(c("maps", "sp", "maptools", "rgdal", "rgeos", "mapproj", "RgoogleMaps", "plotGoogleMaps",

Si en la línea 1856 de tus datos reemplazas |"cotton-wool" spots| por |cotton-wool spots| funciona (y fread lee todas las línas). Suena a bug en el paquete (porque las comillas que no son vecinas de separador las gestiona correctamente). Igual quieres comentarle al autor del paquete el problema para que le eche un ojo. De todos modos, te recomendaría que comprobases el problema

Hello, Thanks for the reply Rui and for pointing out that I forgot to attach my code. Please find attached in this email my code and data. Thanks in advance. Best regards, Pedro Gerhardt Gavronski. On Fri, Feb 23, 2024 at 5:50?AM Rui Barradas <ruipbarradas at sapo.pt> wrote: > > ?s 16:34 de 22/02/2024, Pedro Gavronski. escreveu: > > Hello, > > > > My name is Pedro

Hola Pedro. El INE cambió los ficheros de microdatos no hace mucho, aquí dejo como se haría ahora, (utilizando MicroDatosEs). Lo que cambia es la función para recodificar. http://rpubs.com/joscani/unemplrate El 15/10/14 a las #4, Carlos Ortega escribió: > Hola Pedro, > > Acabo de recordar que hace poco José Luis Cañadas (participa en esta lista) > publicó un enlace suyo a un

The data came through fine, the program was a miss. Can you paste the program into a ".txt" document like a notepad file and send that? You could also paste it into your email IF your email is configured to send text and NOT html. TIm -----Original Message----- From: R-help <r-help-bounces at r-project.org> On Behalf Of Pedro Gavronski. Sent: Friday, February 23, 2024 5:00 AM To:

Hi Im testing wxRuby. I have a problem when ruby try to load the load wx. Im newbie to ruby, maybe a make a simple mistake : Best regards. Install wxRuby-------------------------------------------------------------- pedro@la-vaca-azul:~$ sudo gem install wxruby2-preview Need to update 2 gems from http://gems.rubyforge.org .. complete Select which gem to install for your platform (i486-linux)

Dear R-list, I am using the rpart/mvpart-package for selecting a right-sized regression tree by 10-fold cross-validation. My question: Is there a possibility to find out for every observation in which of the ten folds it is lying? I want to use the same folds for validating another regression method (moving averages) in order to choose the better one. Thanks a lot, Pedro

Hi Pedro, If you keep that same margins for the second plot: par(mar=c(10,0,6,6)) barplot(fhist$counts,axes=FALSE, space=0,horiz=TRUE,col="lightgray") it looks reasonably well aligned to me. Because you are plotting the counts of the values in Simulation, the ordinate (vertical axis) of the bar plot is in quite different units from that of the plot on the left side. Jim On Wed, Jun

Do you think I could use the SelectionDAGBuilder in a pass to accomplish this for me? if there is already code that does this I don't really think I should have to do this again, unless it's not possible to access that code from the stage where I want to use it? On Thu, Dec 8, 2011 at 12:56 AM, Pedro Ferreira <pedro.ferreira at imgtec.com>wrote: > I was unaware of its removal.

I am fitting a logistic model to binary data. The response variable is a factor (0 or 1) and all predictors are continuous variables. The main predictor is LT (I expect a logistic relation between LT and the probability of being mature) and the other are variables I expect to modify this relation. I want to test if all predictors contribute significantly for the fit or not I fit the full