similar to: Pairwise comparison after repeated measures ANOVA

Hi, I am Olga Viedma. I am running a Zero-inflated negative binomial (ZINB) multi-level model using MCMCglmm package. I have a doubt. Can I use the "Liab" outputs as fitted data, instead of the predicted values from "predict"? The liab outputs fit very well with the observed data, whereas the predicted values are so bad. Thanks in advance, Olga Viedma D . Olga

Dear r-help, I would like to select a subset of levels from a factor variable in a data frame and return a data frame. The data set consists of 3 variables, 2 of which are factors (Site, Fish) and one numeric (Datavalue) as follows: Site Fish Datavalue AB 2-1 2.3 AB 2-1 2.4 AB 2-1 2.2 AB 2-2 2.6 AB 2-2 2.5 AB 2-2 2.7 AB 2-3

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?

There''s been some discussion on the IRC channel about getting a site for Merb going. Since the culture around the project is pretty casual, I figured I''d just throw a few ideas together, see what you guys think. So, firstly here are three example pages, just to show off the general look and feel. http://mr-eel.com/tmp/merb/01_frontpage.gif

Hello, I am trying to submit a new package to CRAN. I had checked the packages with R-devel and R-release in windows 7 local and ubuntu 12.04 local. The package has a C file in src folder which successfully compiled in all the cases (there were few trivial warnings). I had submitted the package last week and was asked to resubmit with few changes in the description file. However upon

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 don't want to start the "How many calls can Asterisk handle?" discussion or "How many angels can stand on the point of a pin?" discussion either. But can anyone contribute some practical knowledge of systems that take in channel bank T1s or DS3s from "far away", and process the calls? I am looking for real world, been there, done that, or "check the

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

Dear Help-list, I have a dataframe containing 6 variables, 4 of which are factors, 2 numeric. I want to create another factor variable (SitePos) by combining 2 existing factors (Site and Position). I have tried a number of approaches based on trolling the R FAQs, various R webpages, etc., none of which work. One approach e.g. Data1$SitePos <- paste(Data1$Site, Data1$Position) creates the

Dear Help-list, How can I do a multiple comparison test (mct) on selected contrasts from a linear model while using packages lme4 and multcomp? I am running R 2.13.0 under Windows 7. The following linear model and mct produces a global mct of 15 paired contrasts of the combined (Site, Position) factor SitePos of which only 9 are of interest. Model.G = lmer(log10(SrCa) ~ SitePos + (1 | Eel),

Hi, While investigating the situation of visibility annotations and linkage in libc++ with the goal of removing uses of `__always_inline__`, Eric Fiselier and I stumbled upon the attached test case, which I don't think Clang compiles properly. Here's the gist of the test case, reduced to the important parts (see the attachment if you want to repro): // RUN: %cxx -shared -o

hola, tengo dos vectores de 1134 y 385 elementos que se corresponden con números de accesión de genes. Necesito saber que números son comunes o intersección de vectores. He utilizado intersect(x,y) y me da todo el rato un único valor: V1 V1.1 V1.2 V1.3 V1.4 V1.5 V1.6 V1.7 V1.8 V1.9 V1.10 V1.11 V1.12 1 AJ558305 AJ558305 AJ558305 AJ558305 AJ558305

Hi, Since I got no negative feedback (or any feedback) about the autoconf changes I did a couple days ago, I've released ices 0.1.0 at http://www.icecast.org/releases/ices-0.1.0.tar.gz This is mostly a collection of minor bugfixes, plus an update to the reencoder to use the current LAME API. Note you'll need at least LAME 3.88beta to use reencoding. From the changelog: - 0.1.0

Hello all, Would anybody tell me how to present spatial-temoral point processes in R, for example, I'd like to plot the spatial points in the sequence of their time domain? Cheers -- Pingping Zheng Department of Mathematics and Statistics Fylde College Lancaster University Lancaster LA1 4YF UK

+ Erik, who implemented DR1579 Originally, I had the warning similar to GCC's warning, but took it out due to not having DR1579 implemented in clang (warning changed in r243594) Erik in r274291 implemented DR1579, although PR27785 didn't mention anything about std::move It looks like what's happening is that Clang and GCC handles the return differently. Clang needs the std::move

I'm seeing some discrepancy between our switch traffic logs and our Icecast logs in terms of total amount of data transferred. Our Icecast logs seem to be on the low side. I was wondering how reliable the logging mechanism is in other people's experiences... We're running Icecast 1.3.0 (we can't upgrade because 1.3.10 stutters and dies like crazy on RH6) and Ices. Any info that

Hi, Could someone confirm that the following is NOT possible with Samba: -Run Samba on a Solaris box WITH NO SMBPASSWD -Have the users log in on Windows NT, using the Novell NDS client, and using ENCRYPTION -Have the users authenticated by Samba through NDS (preferred) or NIS. My understanding is that if Samba receives the password encrypted, it has no choice but to use a local smbpasswd (and

Hola tengo una tabla de tipo data.frame "datos" con la siguiente estructura, simplificada en número de niveles por variable, para crear un ejemplo más sencillo: id anio t_8a t_10a t_12a rankf8 rankf10 rankf12 1 1 100 220 220 NA NA NA 2 1 140 350 350 NA NA NA 3 2 55 165 165 NA NA NA 4 2 60 200 200 NA NA NA 5 2 100 NA NA NA NA NA 6 3 NA 350 350 NA NA NA También tengo una matriz,

Hi, since I upgraded to the new Thunderbird version 3.1 I can't establish a TLS/SSL connection anymore. But before the update Thunerbird was able to establish an encrypted session ... Maillog shows me the following now: server dovecot: imap-login: Disconnected (no auth attempts): rip=84.157.133.248, lip=192.168.1.100, TLS handshaking: SSL_accept() failed: error:14094418:SSL

Hi All, We’ve prepared a new function attribute `maynotprogress` and loop metadata `llvm.loop.mustprogress` in order to better formalize the way LLVM deals with infinite loops without observable side-effects. This is deeply related to the implicit forward progress requirements in the IR. Background: There has been a demonstrated need for clarity within the forward progress requirements in LLVM