similar to: R-help archives --- are they up-to-date?

Hola buenas. al final corri el siguiente código en mi máquina de casa. El problema es que ha habido algún cambio en la librería lme4, que hace incompatible los nuevos objetos lmer con la funcioón pamer.fnc. En este tipo de situaciones imagino que lo propio sería ponerme en contacto con el autor o intentar corregir yo mismo el código o incluso ambas. ¿Es decortes escribir al autor reportandole el

Dear all, how can I use vcovHC() to get robust/corrected standard errors from an arima() output? I ran an arima model with AR(1) and got the estimate, se, zvalue and p-value using coeftest(arima.output). However, I cannot use vcovHC(arima.output) to get corrected standard errors. It seems vcovHC works only with lm and plm objects? Is there another way I can get robust/corrected

Estimado Javier Villacampa Gonzáles Si una persona se toma el trabajo de desarrollar y compartir, como es el caso de R, seguramente conoce de problemas que tuvo que decidir mientras escibía el código, y estará agradecido porque hay gente que utiliza su aporte, su trabajo no fué a la nada, todo lo contrario. En mi caso un solo autor tuvo una respuesta negativa para con migo. Javier Marcuzzi

Javier, Creo que aquí aplica la ley de Linus que dice: "Dado un número suficientemente elevado de ojos, todos los errores se convierten en obvios". La persona que revisa y encuentra un error no necesariamente tiene que ser la misma que la que lo escribe. Una motivación muy importante al compartir un código es la de recibir los beneficios del control de calidad por parte de tus pares.

Gracias a todos por las recomendaciones. Ya me he puesto en contacto con el autor y le explicado donde esta el error. Es tan fácil como que han cambiado la estructura de los objetos tipos lmer desde la versión 2.15.3 con lo que su paquete dejo de funcionar. Por suerte la nueva estructura conserva por lo menos la parte necesaria para hacer sus funciones. Se localizar el fallo y como arreglarlo, lo

Hace unos meses os escribir para comunicaros que había un fallo en esta función. Como os prometí os comento la respuesta por si alguno está interesado en utilizar el paquete LMERconvenientsfucntions Dear Javier, The package has been updated and should work for you fine now. Note that function mcp.fnc does not return the fourth plot (dffits) anymore. We still have to figure out a way to compute

It seems that there is a problem in displaying subtitle in general, independently from multi-plot display. when I do plot (c(1,2,3), c(9,8,7), type = "l") title(main = "Main title", sub ="Sub title",cex.main=2, cex.sub = 2) subtitle doesn't get displayed > --- On Sat, 11/21/09, David Winsemius <dwinsemius at comcast.net> > wrote: >

I'm trying to use samba to serve a number of home directories to other samba Linux clients on my network. ie. I have a set of client machines that want to use the /home directory that my samba server will provide I thought I could set up a share that goes something like this [homedirs] path = /home writable = true On a client machine, I smbmount //myserver/homedirs

Hi, In this thread http://lists.digium.com/pipermail/asterisk-users/2008-October/219592.html , I wondered whether SIPPEER curcalls was working. I could test this anew today. Here are my findings : Alice, Bob and Carol ar all using SIP Phones. Whenever Alice is calling Bob, - if Carol is calling Alice, SIPPEER(Alice:curcalls) equals 0 - if Carol is calling Bob, SIPPEER(Bob:curcalls) equals 1

Hi, What is the function for inverse normal transformation? Thanks, Carol [[alternative HTML version deleted]]

Hi. Assume you have an array of person names. I want to generate results in my view that look like this: Abby is your friend or Abby, Bob, and Carol are your friends. So I''d like to say: <%= friends.to_sentence %> <%= pluralize(friends.count, "is") %> your <%= pluralize(friends.count, "friend") %> But because pluralize puts in the

Hi Every one, I have a dataframe "class" with "name", "sex", "age", "height", "Weight". if i caluclate summary statistics with the below code numSummary(class[,c("Height", "Weight")], groups=class$Name, statistics=c("mean", "sd", "quantiles"), quantiles=c(0, .25,.5,.75,1)) iam getting

Hi, Does anyone know how to take a time difference when the format of the time is as 13:22:23.586? I am trying to take the difference of time between stock transactions and need to keep the three decimal places for seconds. I have tried *diff(strptime(x, "%H:%M:%S.000))*, but apperantly that doesn't work. Carol Gao [[alternative HTML version deleted]]

Hi, How to plot the same types of graphics on the same R graphic device? Suppose that we want to plot a vector y against x (using plot for instance). How is it possible to plot y against x for different values of these two vectors on the same device so that the plots could be compared? Cheers, Carol [[alternative HTML version deleted]]

Hi All, PLEASE READ if you depend on Asterisk CDR's and support transfers. Apologies for the shout but I'm desperate to get others to agree Asterisk has a big problem with the CDR's that are generated for transfers. I can understand why not too many people are interested as transfers are complicated and messy. However for those of us having to support transfers and depending on

Hi, I got problem in using title function to create a title for multiple plots presented together by par. As can be seen in the attached file, the title is displayed truncated and the subtitle doesn't get displayed. Here is the code: par(mfrow = c(1,2)) plot(c(1,2,3), c(9,8,7)) plot(c(1,2,3), c(9,8,7)) title(main = "Main title", sub ="Sub title",outer = TRUE,

Hello, How is it possible to sort dates in R? Cheers, Carol

Hello, I am trying to install R/Bioconductor on a G4 Mac running OS X. I have successfully installed R so that a command window opens, but installation of the downloaded Bioconductor package is giving me trouble. After copying/pasting the Bioconductor installation script in to the window and typing getBioC(), I get the following error message. "Error in getBioC(): R not currently

I attempted to contact Drew Tyre, but the email I have for him is no longer in service. Hopefully someone can help. I'm using obs.error in R to model turtle occupancy in wetlands. I have 4 species and 20 possible patch and landscape variables, which I've been testing in smaller groups. > zib.out<-obs.error(y=painted,m=numvis,bp=zvars,pcovar=7) I get the following error

I previously posted about this problem on the r-help list, but was told that I should switch over to the r-sig-fedora list. So here goes: When trying to do a plot in a certain context I got an error: Error in text.default(2, 6, main, cex = cex) : X11 font -adobe-helvetica-%s-%s-*-*-%d-*-*-*-*-*-*-*, face 1 at size 16 could not be loaded This is reproducible (in my current environment)