similar to: Differing Variable Length Inconsistencies in Random Effects/Regression Models

Dear R users, I am trying to fit an lmer model with only random effects which is giving me the following error: Error : length(f1) == length(f2) is not TRUE In addition: Warning messages: 1: In P1L55:family : numerical expression has 390 elements: only the first used 2: In P1L55:family : numerical expression has 390 elements: only the first used I am trying to extract variance components

Dear Mark and the R list, I apologise greatly for not having thanked you earlier for your help with my last mail, which involved trying to sort out variance components using two methods: a direct looped lmer, and a two-step analysis involving regressions of AFLP markers on lmer residuals for the same result (or so I hoped). I was in India for ten days and no CRAN mirror in the vicinity made

---------- Forwarded Message ---------- Date: 18 September 2009 19:24 +0100 From: A Singh <bzwas at bristol.ac.uk> To: William Dunlap <wdunlap at tibco.com> Subject: RE: [R] Error: length(f1) == length(f2) is not TRUE Yup, they are all factors- and its still doesn't work. Getting to the stage where I can use 'summary()' is the problem- the error stalls the process before a

Dear List, I thought it would be much easier to put a second query into a second mail. I need to print 426*10000 blocks of variance components data, where 426 is the number of columns of data that have 10000 permutations of variance generated for each of them. I have tried printing out a smaller number of permutations for a smaller number of markers and that has worked. However, since a

Dear List, I need to print out each of 'k' levels of a factor 'n' times each, where 'n' is the number of elements belonging to each factor. I know that this can normally be done using the gl() command, but in my case, each level 'k' has an unequal number of elements. Example with code is as below: vc<-read.table("P:\\Transit\\CORRECT

Dear all, Does anyone have any suggestions on how to block multiple columns of data one at a time in the midst of an analysis, having specified the blocking variable? I am running a random effects model using lmer, and my data set has multiple columns. Individuals in the study are grouped into 60 families- which is the blocking factor. The random effects are markers (labeled Pxlyy below)

Dear List, I am trying to run a mixed model which, on the R console, prints output as follows: [1] "Marker" [1] "perm no." [1] NA Linear mixed model fit by REML Formula: peg.no.prm ~ 1 + (1 | family/f) Data: modeldf AIC BIC logLik deviance REMLdev 3119 3134 -1555 3112 3111 Random effects: Groups Name Variance Std.Dev. f:family (Intercept) 0.0

Hello, As a result of running linear discriminant analysis, I need to be able to plot the resulting file. I am not sure what the best way to do this is. So far I have tried regular plot("insert_file_name_here") command but the error it gives me is Error in plot.new() : figure margins too large here is sample LDA code I am working with .... library(MASS) example <-

Dear List, I'm having a curious problem with lapply(). I've used it before to convert a subset of columns in my dataframe, to factors, and its worked. But now, on re-running the identical code as before it just doesn't convert the columns into factors at all. As far as I can see I've done nothing different, and its strange that it shouldn't do the action. Has anybody

I HAVE A VERY BIG DATA OF 67 COLMS AND 25000 ROWS AND WOULD LIKE TO TRANSPOSE IT THE R HELP WAS NOT ENOUGH INFORMATION BECOUSE I AM NOT A PROGRAMMER AND FIRST TIME R USER. SO CAN YOU GIVE SOME HINTS OF CODING, AA TT GG GG CC AA TT GG GG CC AA TT GG GG CC AA TT GG GG CC AA TT GG GG CC TO AA AA AA AA AA TT TT TT TT TT GG GG GG GG GG GG GG GG GG GG CC CC CC CC CC [[alternative HTML

I have a couple of hundred American Community Survey Summary Files files containing rectangular arrays of data, mainly though not exclusively numeric. Each file is referred to as a sequence (henceforth "seq"). From these files I am trying to extract particular subsets (tables) consisting of a sets of columns. These tables are defined by three numbers (now in columns in a data frame):

Hello, I am pretty new to R. I am working on neural network classifiers and I am feeding the nnet input from different regions of interest (fMRI data). The script that I am using is this: library (MASS) heap_lda <- data.frame(as.matrix(t(read.table(file="R_10_5runs_matrix9.txt")))*100000,syll = c(rep("heap",3),rep("hoop",3),rep("hop",3))) library(nnet)

Hi, I am using nsd as my authority-only dns servers and find it fast, stable, reliable - all of the design goals you set out for, great. I have a zone that is intended only to be visible to people on my network. In bind I would enforce this with an 'allow-query' configuration directive. I checked the latest 'nsd-for-bind-users' document, and 'allow-query' is not

Dear R users, I am trying to read in a file with 105 columns, and when trying to attach it, get an error as follows: > vc1<-read.table("P:\\R\\Everything-I.txt", header=T, sep=" ", dec=".", na.strings=NA, strip.white=T) > attach(vc1) Error in attach(vc1) : variable names are limited to 256 bytes Is there a way to get around this, and make R accept the

Hello, I use the simplest of examples that somebody can think of in order to generate a matrix of random numbers from within C, calling appropriate R functions. The concrete example is below: --- file pico.c #include <stdio.h> #include <stdlib.h> #include <math.h> #include <R.h> #include <Rmath.h> #define COLM( i, j, m ) ( m*j + i) void pico ( double *y, int n,

Hi everybody, I'm having trouble getting an eval-parse combination to work with lme in a for loop. Consider the code below. The call to lme outside the loop works. The call to aov inside the loop works, but the call to to lme inside the loop does not. The latter fails with Error in model.frame.default(formula = ~meanCol + Group + c1 + c3 + Subject, : variable lengths differ (found for

Has anyone seen any attempts to give Markdown a PEG (parsing expressions grammar) recognizable by some ruby PEG parser generator grammar like Treetop? http://treetop.rubyforge.org/

Dear fellow Markdown enthusiasts, As you all might know, Markdown is a highly context sensitive language that the (often regex-based) syntax highlighting mechanisms in existing editors struggle to keep up with. I have taken John MacFarlane's excellent peg-markdown compiler and modified its parser to function as an interpreter for syntax highlighting purposes. The end result is a

Hi everybody, I'm trying to construct contrasts for an ANOVA to determine if there is a significant trend in the means of my groups. In the following example, based on the type of 2x3 ANOVA I'm trying to perform, does the linear polynomial contrast generated by contr.poly allow me to test for a linear trend across groups? doi=data.frame( Group=c( rep(1, 5), rep(2, 5), rep(3, 5),

Sorry this may well be defined as Off Topic. I apologize in advance. I am interested in performing what I think would be a probabilistic sensitivity simulation. I've done some crude ones before in excel but I'm wondering if R can help me do it more effectively? I have a set of theoretical variables for simplicity lets use (what I think) is an easier example: I have a peg and a hole