similar to: warning message for lmer model with poisson family

Suppose I wanted to plot the barley data like this: dotplot(variety ~ yield | year+ site, barley, strip = strip.custom(style = 4)) The factor levels are far too long for that to be useful. I can overcome that problem if I shorten the levels like this: dotplot(variety ~ yield | year + site, barley, strip = strip.custom(style = 4, factor.levels =

Dear all, I have a question on the use of GAM with repeated measures. My dataset is as follows: - a number of study areas where bird abundance has been determined. Counts have been performed in 3 consecutive years and there were 2 counts per year (i.e. in total 6 counts). - a number of environmental predictors that do not change over year Xi). When using a GLM, a repeated measures design would

Hi R experts, I have been using ReML as follows... model<-lmer(late.growth~mtf+year+treat+hatch.day+hatch.day:year+hatch.day:treat+ mtf:treat+ treat:year+ year:treat:mtf+(1|fybrood), data = A) then I wanted to plot the results of the three way interaction using lines() as follows... tmp<-as.vector(fixef(model)) graph1<-plot(mtf,fitted(f2), xlab=list("Brood Size"),

Hola a todos, Estoy intentando realizar un análisis de medidas repetidas en R y me gustaría pediros opinión porque estoy un poco atascada. Estoy utilizando las funciones lme y lmer. Tengo datos de censos de grillos (efi) realizados semanalmente (week) en tres parcelas (plot) diferentes, cada una de las cuales tiene seis transectos (transect) donde se realizan dichos censos. Los datos son de dos

# Dear R list, # # I'm needing help with lattice, regression and respective lines. # My data is below: bra = gl(2, 24, label = c('c', 's')) em = rep(gl(3, 8, label = c('po', 'pov', 'ce')), 2) tem = rep(c(0, 0, 30, 30, 60, 60, 90, 90), 6) tem2 = tem^2 r = rep(1:2, 24) y = c(40.58, 44.85, 32.55, 35.68, 64.86, 51.95, 42.52, 52.21,

Hola, El problema que tienes es en esta sentencia: ymd_hms(TEM$TIMESTAMP) No la estás asignando a nada. Lo suyo es que hagas algo así: TEM$myTimeStamp <- ymd_hms(TEM$TIMESTAMP) Y crees el intervalo (seq...) on esta nueva variable. El error que obtienes es que "TEM$TIMESTAMP" al leerla del fichero es un string, pero en la lectura automáticamente (salvo que le digas que no lo

Un ejemplo con un modelo más simple: He especificado este modelo: >formula(m2.pile) ger ~ tem + pot + time Si hago predict me da: >predict(m2.pile,newdata=data.frame(tem=25,pot=0,time=3),type="response") 0.08243262 Extraigo los coeficientes: > coef(m2.pile) (Intercept) tem pot time -1.89521331 -0.02303313 4.74499714 0.02043222 Ahora calculo la

As you can see: time is from 1 to 1460, 1 means they the measure Tem for the globe at the first 6 hours , 2=after 12 hours .....; and so on 1460 *6 = 8760 hours which equals 1 year first of all I want to convert all data from kelvin to degree second I want to convert from 6 hourly to daily and make it like this raw data what i want to make time

Dear All, I have conducted an experiment in order to examine predation pressure in the surroundings of potential wildlife road-crossing structures. I have documented predation occurrence (binary?) in these structures and calculated several possible explanatory variables describing the spatial heterogeneity in several scales. At the landscape scale I have calculated the percentage of different

This is no doubt trivial but after searching the help files and the web, I cannot seem to find it. 1) How do I convert 'hgt' into 'HGT' in R? 2) How should I have used the help facilities to find this? At the end of the day, all I want to do is case insensitive string matching... i.e. 'if ("HGT" == 'hgt') print('this should be true')' I tried

Estimad en s errer en s He hecho este modelo glm m1.pile<-glm(ger~tem+pot+time+I(tem^2)+I(tem^2):pot ,family="binomial" ,data=long.PILE ) Que nos da la probabilidad de germinación de una semilla en función de tem (Temperatura), pot (Humedad del suelo) y time (Tiempo que la semilla pasa en esas condiciones). Ahora quiero, para diferentes tem, pot

Oh where to begin? I work for a company that has not taken our Linux or Unix server configuration management very seriously. We are primarily a Windows shop but long story short upper-management has taken the advice of our Windows Server admins and purchased BigFix which is now IBM Tivoli Endpoint Management (TEM). I have tried for over a month to get one of my SuSE Linux Enterprise Server (SLES)

Oi F?bio. Pelo IP tamb?m n?o acessa. A minha suspeita de que era um problema ligado a interface est? errada. Aleat?riamente usu?rios nas duas interfaces perdem o acesso. Ainda assim, somente usu?rios com Windows XP tem esse problema. Vou testar suas indica??es. F?bio Rabelo wrote: > > > 2010/3/23 Leonardo Carneiro - Veltrac <lscarneiro at veltrac.com.br > <mailto:lscarneiro

The US water year extends from 01 October yyyy-1 through 30 September yyyy and is referenced by the year starting on the included 01 January yyyy. I'd like to be able to find the annual means for the water year. To do so I've taken the input date-time, which is in the usual format "1991-10-07 10:35:00" changed it by: w$d<-as.POSIXct(w$date.time) Now I can add an

Hello, I have produced some segmented regressions with the segmented package by Viggo Mutteo. I have some example data and code below. I want to annotate the individual segments with the slope parameter (actually it would be nicer to annotate with 1000*slope and add some small amount of text as well). How can I do it? Reading the docs for segmented I can access all of the slope parameters via a

Dear All, After playing with lmer for couple of days, I have to say that I am amazed! I've been using quite some multilevel/mixed modeling packages, lme4 is a strong candidate for the overall winner, especially for multilevel generzlized linear models. Now go back to my two-level poisson model with cross-classified model. I've been testing various different model specificatios for the

Hi All, I am using the glm.nb function from the MASS package (current version) to fit and compare GLMs with negative binomial error distributions. My question is: what is the appropriate method to use in the anova function to compare models? If only one fitted object like m<-glm.nb(number<-p+sal+temp,data=data) is specified in the anova function (anova(m)), a fixed theta is used to

Sauda??es. Como que se faz um conector RJ45 em uma ponta e RJ11 e outra. Pretendo conectar a linha de um ATA em uma placa Khomp KFXO IP. A ponta que tem o conector RJ45 est? crimpada com a sequencia 568B e vai ser conectada na placa Khomp, mas a ponta RJ11 eu n?o sei como deve ficar. Li alguns manuais na internet mas n?o entendi ao certo como tem que ser feito. -- Att.* *** Luis H. Forchesatto

Dear Group: I have been trying to write an R function within which the function groupedData() would be used. The following is a sample program: ####################### sid<-rep(1:6,times=2) time<-c(1:12) trt<-rep(letters[1:3],times=4) tem<-data.frame(sid,time,trt) test1.fun<- function(dat=dat) { groupedData(time ~ trt|sid, data=dat) } test1.fun(dat=tem)

Hello, I migrate my server POP3/IMAP recently the Courier-IMAP for Dovecot and am having problems when the user configures your e-mail client to leave a copy of mail on the server. All mail in the inbox, are brought all the time that you check your mailbox input via POP3, causing a large traffic of unnecessary data every time. In my opinion you should bring only the new mail, unread. Can you