similar to: Help with lmer formula

Hi R users! I am trying to learn some multilevel analysis, but unfortunately i am now very confused. The reason: http://www.ats.ucla.edu/stat/hlm/seminars/hlm_mlm/mlm_hlm_seminar.htm http://www.ats.ucla.edu/stat/sas/seminars/sas_mlm/mlm_sas_seminar.htm and MlmSoftRev. pdf from mlmRev package. >From what i see, the first two links seem to declare the level one variable as a random part (i

I would like to fit a hierarchical regression model from Witte et al. (1994; see reference below). It's a logistic regression of a health outcome on quntities of food intake; the linear predictor has the form, X*beta + W*gamma, where X is a matrix of consumption of 82 foods (i.e., the rows of X represent people in the study, the columns represent different foods, and X_ij is the amount of

Hi, I was using the 2.4.1 R version and I had no problem saving my plots as postScript. Now that I have installed the latest version 2.6.0 I can not save any plot as postScript. When I try the following message appears: Erro: Invalid font type Al?m disso: Warning messages: 1: font family not found in PostScript font database 2: font family not found in PostScript font database What should I

Buenas tardes a tod@s, Estoy incursionando en el analisis de tweets utilizando el paquete twitteR y siguiendo http://www.webmining.cl/2012/07/text-mining-de-twitter-usando-r/ Desafortunadamente cuando ejecuto # cargar librerias library(twitteR) library(tm) library(wordcloud) # recolecta tweets de @camila_vallejo tweets = userTimeline("camila_vallejo", 2000) obtengo Error in

Dear List, I am new to TS-Modeling in R. I would like to fit an ARMA(1,1) model for a balanced panel, running Y on a full set of unit and year dummies using an arma(1,1) for the disturbance: y_it=unit.dummies+yeardummies+e_it where: e_it=d*e_it-1+u_it+q*u_it-1 How can I fit this model in R? arma() does not seem to take covariates (or I don't understand how to specify the function so that

Hi everyone, I am interested in estimating this type of random effects panel: y_it = x'_it * beta + u_it + e_it u_it = rho * u_it-1 + d_it rho belongs to (-1, 1) where: u and e are independent and normally zero-mean distributed. d is also independently normally zero-mean distributed. So, I want random effects for group i to be correlated in t, following an AR(1) process. I am

Hello, I am interested in estimating this type of random effects panel: y_it = x'_it * beta + u_it + e_it u_it = rho * u_it-1 + d_it rho belongs to (-1, 1) where: u and e are independently normally zero-mean distributed. d is also independently normally zero-mean distributed. So, I want random effects for group i to be correlated in t, following an AR(1) process. Any idea of how

Dear r-helpers, I am working with plm package. I am trying to fit a fixed effects (or a 'within') model of the form y_it = a_i + b_i*t + e_it, i.e. a model with an individual-specific intercept and an individual- specific slope. Does plm support this directly? Thanks in advance! Otto Kassi

Hello, I am using {plm} to estimate panel models. I want to estimate a model that includes fixed effects for time and individual, but has a random individual effect for the coefficient on the independent variable. That is, I would like to estimate the model: Y_it = a_i + a_t + B_i * X_it + e_it Where i denotes individuals, t denotes time, X is my independent variable, and B (beta) is the

The compositional data xi=(x_i1, x_i2,..., x_in), for each fixed i , xij>0, and sum(xij)=1; I want to compare the mean( u_i) of several groups i.e. H0: u_1=u_2=...=u_N or Hj0: u_1j=u_2j=...=u_Nj Are there any ANOVA tpye tools to do this work in R? Thanks, WEN S Q [[alternative HTML version deleted]]

Estimado Jorge Velez Hace un tiempo yo intente sobre algo parecido, pero utilizando nuestros propios correos, lástima que estoy con otros problemas para resolver y quede desactualizado, sin embargo por ahí leí que hay una modificación de política por parte de twitter, muchos lectores quedaron si acceso, esto puede ser el error en tweets = userTimeline("camila_vallejo", 2000) Aunque

Hola Jorge. Para poder correr las funciones del library twitteR tienes que tener varias cosas. Una cuenta en twitter y con esa cuenta inscribirte como developer en dev.twitter.com. Posterior a esto tienes que inscribir un app dentro de twitter y con las claves que se otorgan tienes que correr la función de autorizacion que ya te indicaron. Varias aclaraciones, la librería puede generar problemas

Hi all, I am using quantile regression with svy design. I want to retrieve summary regression statistics (std error, p-value), since I don't have any in my output: Commands: clus1_d<- svydesign(id=~cd002_co, weights=~wtper, strata=~str, data=data) bclus1<-as.svrepdesign(clus1_d,type="bootstrap",replicates=100) fit1<-