similar to: ZINB multi-level model using MCMCglmm

Dear colleagues: I am working with binary mixed models using lme4 package and everything works fine until I try to spatialize the fitted models using grids predictors. I have used the "predict" function from RASTER package and the RSAGA functionalities, but always I have the same error: "Error in object$xlevels : $ operator not defined for this S4 class". I think the

Hi! I am trying to know which habitat variables most affect bird counts in a radius of 100m. I obtained bird counts in 2751 spatial points, and measured percentage of 21 habitat variables in these points. I applied a mixed model using the "lmer" function to these data, but I do not know how to select the best model using AIC here. Is there a way to do this automatically with R?

Hola: Soy nuevo en esto de R, llevo toda la mañana perdida con esto, a ver si alguien me puede ayudar, Cuando genero el gráfico con el siguiente archivo: library(gplots) data <- read.table("data.dat",h=T) postscript(file = "join10kb.ps", horizontal = FALSE, onefile = FALSE,family = "Helvetica", width = 5, height = 4) par(mar=c(5.1, 4.1, 0.1, 0.1)) d <-

I am analyzing some data obtained after measuring some parameters at different times in samples obtained from many subjects. The model is quite simple: aov(parameter ~ Time + Error(Subject/Time)) Now I want to make a pairwise comparison between the levels of Time. However, I have not find how to do such a thing. I cannot use TukeyHSD or pairwise.t.test, I supposse. Maybe using contrasts? Could

Hola buenos días, me presento, me llamo Miguel y 'soy de' y 'vivo en' Galicia. Soy profesor de secundaria (Bachillerato Adultos) y llevo 15 días estudiando R a un buen ritmo, pero todavía me faltan miles de cosas. He visto que R facilita, no solo el análisis de datos y que posee una potencia en cálculos estadísticos a cualquier nivel, sino gran caudal de recursos para Data Mining,

Hola amigos, buenos días Esto sigue avanzando. Estoy ahora tratando de establecer la conectividad de R con Excel y he seguido dos viídeos en YouTube de una profesora llamada Bebilda que se explica muy bien. Con ambos tengo el mismo problema así que planteo uno de ellos, éste es: https://youtu.be/_pSJQO_9I4k?list=PL7DA3FDA21A1A6310 Pero os lo explico brevemente: Cargo el paquete RODBC >

Hi All, I have a very simple problem. I have several files in a same directory. I would like to send for an object only the files that finish in ".sens.". I execute the command below, files <- dir(pattern="*.sens") but it includes all of the files that have "sens", independent of they be in the end or in the middle of the name of the file. How could I solve

Dear all, I generated a covariance matrix and I would like to generate a 1D plot of the data that auto-correlate. any suggestions? Thanks, Marcelo -- Marcelo Andrade de Lima UNIFESP - Universidade Federal de São Paulo Departamento de Bioquímica Disciplina de Biologia Molecular Rua Três de Maio 100, 4 andar - Vila Clementino, 04044-020 Lab +55 11 55764438 R.1188 Cell +55 11 92725274

Hello, I'm trying to write the following matlab code into R: N = zeros(n-1); for i=2:(n-1) N(1,i) = 1/(pi * (i-1)); end for i=2:(n-2) for j=i:(n-1) N(i,j) = N(i-1,j-1); end; end for i=2:(n-1) end for j=1:i N(i,j) = -N(j,i); end; any suggestions? Thanks can i just add the following line to my calculation N=1/(pi*(i-1) -- Marcelo Andrade de Lima UNIFESP - Universidade Federal

Hello, I would like to plot the results of a LDA analysis plotting the discriminant scores with the decission boundaries on it with rggobi. I have GGobi already installed on my computer. I have three classes, so the plot would be LD1xLD2 plus the decission boundaries. Here there is the code I use make the plot: library(MASS) data <- zgcppr273K.pca$x[,1:7] Tumor <-

Dear Dimitris, I have exactly the same problem than you, Do you get some solution? Thanks, Lola Lola Gadea Profesora titular de Economía Aplicada/Lecturer in Applied Economics Universidad de Zaragoza/University of Zaragoza (Spain) lgadea@unizar.es <http://estructuraehistoria.unizar.es/personal/lgadea/index.html>http://estructuraehistoria.unizar.es/personal/lgadea/index.html Grupo de

Dear Dimitris, I have exactly the same problem than you...Do you have some any solution? Thanks, Lola Lola Gadea Profesora titular de Economía Aplicada/Lecturer in Applied Economics Universidad de Zaragoza/University of Zaragoza (Spain) lgadea@unizar.es <http://estructuraehistoria.unizar.es/personal/lgadea/index.html>http://estructuraehistoria.unizar.es/personal/lgadea/index.html Grupo

Dear Sr. of R-project: In the document "R: Regulatory Compliance and Validation Issues. A guidance Document for the Use of R in Regulated Clinical Trial Environments", in the point 6.3 say that exist a set of validation test, that are localized in the test subdirectory. I search the place in the page www.r-project.org and not found this, can you say me where is the set?. I work in a

Buenas, Sigo aprobechando par aaprender en las vacaciones de alumnos y nuevamente recurro a Uds. Sé que la perspectiva puede ser erronea, o sea que les pido que entiendan lo que quiero decir, más que lo que digo. Tengo una lista de matrices todas de idénticas dimensiones, y necesito hacer operaciones entre ellas. Ocurrencia: creí que podría hacer un array, de tal manera que , así como puedo

Muchas gracias, Javier, por tu respuesta. Me sobrevalorás. Mi última clase de matemática como tal fue en 4º de secundaria, y era la matemática que nos daban a los de letras. Ni te cuento hace cuántos años. No puedo seguir la mayoría de lo que comentás sin hacer un trabajo forense. Así que agarré lo que entendí... y me parece que, o lo que yo quiero hacer son matrices de tres dimensiones (si es que

Dear Mr.Zeileis & all. (1)     Thx for your reply. Yes, I am talk about the function zeroinfl() from the package "pscl". I want to use Newton Raphson to get parameter             estimation ZINB, so I try this: ----------------------------------------------------------------------------------------------------------------------------------         > zinb <- zeroinfl(y

Dear users, Does anyone have any idea how to generate standardised beta coefficients for a ZINB model in R to compare which explanatory variables are having the greatest impact on the dependent variable? Thanks, Jeremy [[alternative HTML version deleted]]

Dear all, I am modeling the incidence of recreational anglers along a stretch of coastline, and with a vary large proportion of zeros (>80%) have chosen to use a zero inflated negative binomial (ZINB) distribution. I am using the same variables for both parts of the model, can anyone help me with R code to compute overall marginal effects of each variable? My model is specified as follows:

Dear all.. I have a respon variable y. Predictor variable are x1, x2, x3, x4, x5 (1) What is the syntax to get paramater estimation of ZINB Model by Newton Raphson (not BFGS) (2) What syntax to plot probability of observed & predicted of ZINB Thx. Regards Krist. [[alternative HTML version deleted]]

I want to generate 1,000 random samples of sample size=1,000 from ZINB. I know there is a rnegbin() to generate random samples from NB, and I know I can use the following process: do i=1 to 1000 n=0 do i=1 to 1000 if runi(1)>0.1 then x(i) = 0; else x(i)=rnegbin(); n=n+1; if n>1000 then stop; end; output; end; Anybody can help me out with the R code? Thanks very much ahead of time.