similar to: Proc Mixed to R

Dear R Community, I am attempting to transpose a dataset from rows to columns but am stuck. I have tried using reshape() with little luck, possibly due to the categorical nature of the data. For example: id<-c(1,2,2,3,3,3) author<-c("j","k","k","l","l","l")

Hello, I have a couple questions regarding generalized linear mixed models specifically around fitting the random effects terms correctly to account for any pseudo-replication. I am reading through and trying to follow examples from Zuur et al. Mixed Effects Models and Extensions in Ecology with R, but am still at bit unsure if I am specifying the models correctly. Background information: Our

I have a nested ANOVA, with a fixed factor "tmt" nested within "site" (random). There are missing values in the data set. aeucs, tmt and site have been defined as objects I have tried: model1=lme(aeucs~tmt,random=~1|tmt/site) I get the following error message Error in na.fail.default(list(aeucs = c(0.833333333, 1, 1, 1, 1, 1, 0.666666667, : missing values in object I

I have a local library which depends on the expm library. The expm library is loaded into my personal space and I have the environment variable R_LIBS_USER set appropriately. The command "library(expm)" works just fine from the command line, and in fact the package works if I do the source() and dyn.load() commands by hand. The following sequence works: tmt% R CMD build

Dear all, I would like to use the values in vegaggr.BLMCMR02$colony str(vegaggr.BLMCMR02) `data.frame': 1678 obs. of 3 variables: $ vegtype : Factor w/ 27 levels "2010","2020",..: 3 4 5 19 4 5 19 5 $ colony : Factor w/ 406 levels "0","1","10","100",..: 1 1 1 1 2 2 2 $ Totvegproparea: num 0.00055 0.03956 0.95705

Martin, I am aware of --vanilla; I use it myself for some testing. In this case R_LIBS_USER was set externally (part of my login) and does not involve any of the R scripts. That means it is inherited by any subprocess. For example: tmt1495% R --vanilla --no-environ R version 3.3.1 (2016-06-21) -- "Bug in Your Hair" Copyright (C) 2016 The R Foundation for Statistical Computing

Dear All, I wonder if anyone can please offer any advice on a model including 2 fixed effects and 1 random effect, as well as a covariate? The experimental design is as follows: I have a two by two factor design, where the two factors, Age (A) and Group size (G), both have 2 levels (old or young, and 1 or 3 respectively), and I am interested in the effect of these factors upon a continuous

I have a data frame containing the Id, Mother, Father and Sex from about 10,000 animals in our colony. I am interested in graphing simple family trees for a given subject or small number of subjects. The basic idea is: start with data frame from entire colony and list of index animals. I need to identify all immediate relatives of these index animals and plot the pedigree for them. We're

Hi, I have a data set with 999 observations, for each of them I have data on four variables: site, colony, gender (quite a few NA values), and cohort. This is how the data set looks like: > str(dispersal) 'data.frame': 999 obs. of 4 variables: $ site : Factor w/ 2 levels "1","2": 1 1 1 1 1 1 1 1 2 2 ... $ gender: Factor w/ 2 levels "0","1":

I'm not quite sure whether this is a plain browsing problem, or a wins problems, or ... I had wins running on one machine (Samba(2.0.5a)/SCO Unix), but my main file server on another (Samba (2.2.1a)/Linux). I moved the wins support to the Linux server, and at the same time I made the Linux server the Domain Master Browser, but now some devices have dissapeared from the network. The devices

Hi, I have been trying to re-order several items in a trellised barchart display in lattice, but can't seem to figure it out. ###sample code, Stage and Colony have 2 and 3 levels respectively. barchart(Activity ~ Percent | Stage + Colony, data = Percent.df, horizontal = TRUE, layout = c(2,3), xlab = "Percent Time Engaged in Activity", ylab = "Activity") This

I have a local library with functions that interrogates an institution-specific web API, so is not of interest to anyone outside of Mayo. For some reason the R CMD INSTALL command fails. See below: Build the library, then install it. tmt-local2127% R CMD build dart * checking for file ?dart/DESCRIPTION? ... OK * preparing ?dart?: * checking DESCRIPTION meta-information ... OK * installing

Trying to plot multiple lines from a simple matrix with headers (eight observations per column). I will be doing a number of these, with varying numbers of columns, and do not want to enter the header names for each one (I got frustrated and just wrote them out, which did work). Data reads fine, first plot is fine, but when i use the code at the bottom for a for i loop it tells me that x and y

A user asked me about this and I can't figure it out. tmt% R R Under development (unstable) (2018-04-09 r74565) -- "Unsuffered Consequences" Copyright (C) 2018 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) > library(survival) > data(cgd0) Warning message: In data(cgd0) : data set ?cgd0? not found ---- The data set is present and can be

Greetings, I have "inherited" a cDNA macroarray dataset that is structured as follows. Three different stressors were tested. For each stressor, there are two treatments (control and stressed). For each treatment, two biological replicates exist, and these are paired (i.e., there is a stressed array for colony A and a control array from this same colony). For one of these samples,

I've downloaded Samba 2.2.6 (binaries) from SCOs Skunkware site, but when I try and run smbclient I get a message : dynamic linker : smbclient : error opening /usr/local/lib/libncurses.so.4 Killed I've checked and there is no libncurses.so.4 (or libncurses.<anything> for that matter) on the system. I've also failed to find anything relevant in either the list archive or SCOs

This is likely yet another instance of round off error, but it caught me by surprise. tmt% R --vanilla (headers skipped, version 3.0.2 on Linux) > load('qtest.rda') > length(temp) [1] 3622 > max(temp) >= quantile(temp, .98) 98% FALSE I can send the file to anyone who would like to understand this more deeply. The top 3% of the vector is a single repeated value. Terry

Hello, I am Mudit Raj Gupta, fourth year student of M.S. (Hons.) Chemistry and B.E. (Hons.) Electronics and Instrumentation at BITS-Pilani ( http://www.bits-pilani.ac.in/). I am interested in *Machine Learning and Computaional Intelligence*. I have an interest in implementing various existing algorithms and developing modified/new algorithms related to machine learning and computational

Hello, I have been attempting to run a poisson glmm using lme4 for some time now and have had a lot of trouble. I would say 9 times out of 10 I receive the following warning: CHOLMOD warning: %h Error in mer_finalize(ans) : Cholmod error `not positive definite' at file:../Cholesky/t_cholmod_rowfac.c, line 432 My data are counts of microbe colony forming units (CFUs) collected from

I am getting failure of build and check, for an Rd file that has a long argument list. Guess diagnosis: a quoted string beyond a certain point in the argument list is fatal. Example: Use the function below, create an Rd file for it with prompt(). Move the .Rd file to the man directory (no need to edit it) and try building dart.control <- function(server=c("production",