similar to: offset en bucle

Amigos de la erre. He creado mi primer bucle con for para entrenar unos modelos con GAM. La respuesta es quasipoisson porque estoy trabajando con densidades de peces. Sin embargo, tengo un problema, no se muy bien como añadir el offset a la formula siguiente cuando creo el bucle. GAM.A1 <-gam ((DYO)~s(DMA,k=4)+ s(WOD,k=4)+s(CIN,k=4)+s(DRA,k=4)+s(DBR,k=4)

Hola, buen día. Una consulta, no sé conozco en extenso todas las bondades de R, pero quería saber exactamente si hay algún paquete que permita realizar análisis ecológicos específicamente para comunidades marinas. Existen muchos índices, indicadores, estimadores en general en muchos programas que están a nivel de usuario, pero quería saber si todos estos pueden estar ya disponibles en R.

Estimados por solicitar su ayuda, estoy intentando cambiar el color de una etiqueta de un gráfico, he realizado múltiples intentos sin un resultado satisfactorio. Debido a que necesito realizar múltiples gráficos he realizado una función que dando parámetros de entrada me grafique las variables o los valores que deseo. Mi problema es en la representación de los colores. Los parámetros de la

Hola a todos, y gracias por leer mi mensaje: Llevo utilizando unos meses R, vengo de SPSS, y estoy muy sorprendido con su capacidad. Mi problema, es que estoy intentando hacer una comparación entre tres sistemas de predicción: una regresión logística, un árbol de clasificación y una red neuronal. De momento, he podido hacer las comparaciones con los dos primeros sistemas, pero mi problema es la

Buenos días. Tengo dos mapas que hago con ggplot (abajo) y necesito superponer el 2º al 1º. En el 1º se representa la variable PA, binomial (0, 1), en gris, y en el 2º la variable ExtTG7085, que es numérica y va de 0 a 200. La 2ª es un gradiente en rojo. Lo ideal sería que donde coinciden el rojo de la 2ª se viera más oscuro, pero supongo que eso no es fácil. Podría también ser en otro color, pero

Buenos días. Hago mis mapas con el código que os copio abajo, pero me gustaría que me dibujase las líneas del mapa después de poner los puntos para que se vean. Puedo hacer los puntos más pequeños, pero entonces quedan separados y yo quiero que rellenen el mapa. Gracias, como siempre, Manuel print(ggplot(legend=FALSE)+geom_path(data=map_data('world'), aes(x=long,

Hi, I just install 2.2.2 in a Red Hat 7.2 with kernel 2.4.16 on ext 3 + lvm. The problem are the locks!!!!!! I've 5 users working in .dbf + cdx (foxpro) files with an 16 bits application. It accesses to the files via DAO 2.0. A lot of times in the days the users gets 'currently locked' messages when the application open the dbf's files with read only access! I desactivated the

I have a logistic model fitted with the following R function: glmfit<-glm(formula, data, family=binomial) A reasonable cutoff value in order to get a good data classification (or confusion matrix) with the fitted model is 0.2 instead of the mostly used 0.5. And I want to use the `cv.glm` function with the fitted model: cv.glm(data, glmfit, cost, K) Since the response in the fitted

I'm working with samba 1.9.18p10-7, on a i586 with debian linux v2.1. I'm trying to personalize windows95/98 on a multiuser/shared PC instalation and I want to include a Historial folder somewhere in the home directory of every user on a centralized server. This way any user could log on from any PC (in a classroom for example) an may have their own configuration for the Internet

Hi, I am trying to add multiple variance structures such as the first example below: vf1 <- varComb(varIdent(form = ~1|Sex), varPower()) However my code below will not work can anybody please advise me? VFcomb<-varComb(varExp(form=~depcptwithextybf),varFixed(form=~FebNAO)) also if you have two variables with the same weights function would you write that as:

Dear R-help members, I have a question regarding how to use varComb function to specify a variance function for the "weights" in the gls. I need to fit a linear model with heteroscedasticity. The variance function is exp(c0+nu0*W +nu1*W^2) where W is a covariate. Initially I want to use varFunc to define my own variance function following the instruction in the Pinheiro and

Dear R-help members, I have a question regarding how to use varComb function to specify a variance function for the "weights" in the gls. I need to fit a linear model with heteroscedasticity. The variance function is exp(c0+nu0*W +nu1*W^2) where W is a covariate. Initially I want to use varFunc to define my own variance function following the instruction in the Pinheiro and Bates

Dear list, I was trying to create a VarClass for nlme to work with Fay-Herriot (FH) models. The idea was to create a modification of VarComb that instead of multiplying the variance functions made their sum (I called it varSum). After some fails etc... I found that the I was not getting the expected results because I needed to make sigma fixed. Trying to find how to make sigma fixed I run into

Hello, I am trying to design new variance structures by using fixed effects variables in combination with the VarPower function. That is, I would like to create and evaluate my own variance function in the data frame and then incorporate it into the model using varPower, with value=.5. As a start, I am trying to recreate the function of VarConstPower by introducing two new variables in the

Dear All, I am using the function gls (in the nlme package) and I would like to fit a heteroscedastic model, with different variances for each of the levels of two stratification variables. In p. 210 of Pinheiro & Bates ("Mixed effects models in S and S-Plus", 2000, Springer), the authors show the use of the "*" operator. However, that is not what I want, because it

 Hello everyone, Apologies in advance, as this is partially a stats question and partially an R question.  I have been using a GAM to model the activity level of bats going into and coming out from a forested edge.  I had eight microphones set up in a line transect at each of eight sites, and I am hoping to construct a model for each of 7 species.  My count data has a reverse J-shaped skew and

https://bugs.freedesktop.org/show_bug.cgi?id=33999 Summary: 2.6.37 - NV11 crashes X if glxgears maximised Product: xorg Version: unspecified Platform: x86 (IA32) OS/Version: Linux (All) Status: NEW Severity: normal Priority: medium Component: Driver/nouveau AssignedTo: nouveau at

Hello everybody, I have a dataset where each row has number of subjects and that gives me natural weights for the variance function. Additionally I see that variance increases with Age, which is a regressor. So using gls I have weights=varFixed(~1/n) but don't know how to include the extra effect of the regressor. Fitted values show a quadratic curve vs. age, not sure if that helps.

Hello all. I was wondering if the logLike.nls() and logLike.nlme() functions are still being used. Neither function seems to be available in the most recent release of R (1.9.1). The following is contained in the help file for logLik(): "classes which already have methods for this function include: 'glm', 'lm', 'nls' and 'gls', 'lme' and others in

  Hola,   Por favor, ¿podríais indicarme qué recursos (librerías o ideas) pueden resultar de utilidad para crear un gráfico del estilo del de la figura 3.8 del siguiente link?   http://www.tsc.uvigo.es/BIO/Bioing/ChrLDoc3.html#3.5   Actualmente estoy utilizando funciones muy básicas y la verdad es que no me encuentro muy satisfecha con el resultado.   Muchas gracias.   Eva [[alternative HTML