similar to: applying cbind (or any function) across all components in a list

Hello list: I am using "cut" and "table" to obtain a frequency table from a numeric sample vector.  The idea is to calculate mean and standard deviation on grouped data.  However, I can't extract the midpoints of the class intervals, which seem to be strings treated as factors.  How do i extract the midpoint? Thanks, jose loreto [[alternative HTML version deleted]]

Hello: I'm having quite a bit of difficulty adjusting the x-labels using bargraphCI. I've tried using text and srt=45 to rotate the labels or mtext for 2 lines to break up the labels. However, using either method, I cannot line up the labels with the midpoints of the bars (they line up with some sort of tick mark that is off the midpoint of the bars). Any suggestions would be greatly

I am using dgc.genetics to perform TDT analysis on SNP data from a cohort of trios. I now have a file with about 6008 variables. The first few variables related to the pedigree data such as the pedigree ID the person ID etc. Thereafter each variable is a specific locus or marker. The variables are named by a pattern such as "Genotype.nnnnn" with nnnnn corresponding to a number which

Dear R-users, I need some help concerning the plotting of dendrograms for hierarchical agglomerative clustering. The agglomeration niveau of each step should be displayed at the branches of the dendrogram. For this I need the x/y coordinates of the branch-agglomerations of the dendrogram. The y-values are known (the heights of the agglomeration), but how can I get the x-values? > mydata

Hello r-help, I have a lengthy vector of data (with values anywhere from 1-200), and another index vector of 'groups' representing values 0-2, 3-5, 6-8, ... of length 67. The index vector has the structure (1, 4, 7, ... , 196, 199), where each value is the midpoint of each respective group. I'm trying to convert the data vector such that values falling into each group are changed to

SYSLINUX 3.11-pre5 is hereby a bug fix release candidate. If I don't receive any problem reports, I will release it as-is on Tuesday, August 30. http://www.kernel.org/pub/linux/utils/boot/syslinux/Testing/ -hpa

Hi, I need to construct a formula programaticly, and pass it to a function such as the linear mixed model lme. The help says it requires "a two-sided linear formula object describing the fixed-effects part of the model" but I do not know how to create this formula. I have tried various things using formula(x, ...), as.formula(object, env = parent.frame()) and as.Formula(x, ...)

Hello, I am working with a matrix of multilocus genotypes for ~180 individual snail samples, with substantial missing data. I am trying to calculate the pairwise genetic distance between individuals using the stats package 'dist' function, using euclidean distance. I took a subset of this dataset (3 samples x 3 loci) to test how euclidean distance is calculated: 3x3 subset used

Bottom Line Up Front: How does one reshape genetic data from long to wide? I currently have a lot of data. About 180 individuals (some probands/patients, some parents, rare siblings) and SNP data from 6000 loci on each. The standard formats seem to be something along the lines of Famid, pid, fatid, motid, affected, sex, locus1Allele1, locus1Allele2, locus2Allele1, locus2Allele2, etc In other

Hello again I'm trying to put error bars onto a barplot. I've tried something that Bill Simpson suggested a while ago, ie: x<-c(1,2,3,4,5) y<-c(1.1, 2.3, 3.0, 3.9, 5.1) ucl<-c(1.3, 2.4, 3.5, 4.1, 5.3) lcl<-c(.9, 1.8, 2.7, 3.8, 5.0) plot(x,y, ylim=range(c(lcl,ucl))) arrows(x,ucl,x,lcl,length=.05,angle=90,code=3) #or segments(x,ucl,x,lcl) but I can't get it to work on a

#Hello, #I have a question about obtaining results from a loop I have written. #Below is a sample of individual genotypes from a genetic question I am working on called "P.genotype.sample ". P.genotype.sample<-matrix(10,10,10) P.genotype.sample[,1]<-c(2,2,1,5,1,1,5,6,1,3) P.genotype.sample[,2]<-c(6,3,3,6,8,1,6,7,2,3) P.genotype.sample[,3]<-c(2,2,2,3,3,2,2,2,3,3)

Hi dear R-listers, I'm trying to fit a 3-level model using lme in R. My sample size is about 2965 and 3 factors: year (5 levels), ssize (4 levels), condition (2 levels). When I issue the following command: > lme(var~year*ssize*condition,random=~ssize+condition|subject,data=smp,method ="ML") I got the following error: Error in logLik.lmeStructInt(lmeSt, lmePars) :

Hi , is it possible to calculate ld-measures D, D', r and perhaps corresponding p-values with r IF THE PHASE IS KNOWN? The genetics - package provides the LD function only for ambigious phase. Thank you very much Bettina Kulle

Hello, I have a data-frame in which one-column is a factor: > str(data); `data.frame': 194 obs. of 8 variables: $ Type : Factor w/ 3 levels "Nuclear-Rec..",..: 1 2 2 2 2 2 2 2 2 2 ... $ Locus : num 0.000571 0.004000 0.001429 0.004857 0.007429 ... And I'd like to make a boxplot of the data$Locus values, where each level of the factor gets its own

Hi, I'm trying to make a chromosomal map in R by using the locus. I have a list of genes and their locus, and I want to visualise that so you can see if there are multiple genes on a specific place on a chromosome. A example of what I more or less want is below: http://www.nabble.com/file/p21474206/untitled.JPG untitled.JPG The genes and locus are here:

Hello everyone: I'm a new member of this group. I have a question about "oneway.test". When I use "anova(lm(....))" to analysis the ANOVA, I can get the information about Sum Sq and Mean Sq. (The R code and the results are as follows.)

Hi All, I would like to ask the following. I have an array of data in an objetct, let's say X. I need to use a for loop on the elements of one or more columns of X and I am having a debate with a colleague about the best memory management. I believe that if I do: col1 = X[,1] col2 = X[,2] ... colx = X[,x] and then for(i in whatever){ do something using col1[i], col2[i] ... colx[i] }

I know that there are quite a few packages out that there for cluster analysis. The problem that I am facing is finding a package that will not incorporate all my samples into clusters but just the samples that fit a threshold (that I have not set yet and may need help finding the right level) for genotyping. It should be able to "no call" samples outside the clusters. It also needs to

Hello I have a vector gd of 279 values of diameter of limpets. Using 'hist(gd, breaks=seq(10,60,5),plot=F)' I can obtain a vector of 10 centers of size-class ($mids). But I'd like to work on a vector of 279 values being the mids (centers of size-classes) corresponding to the respective 279 values in gd. I looked in the different manuals and in archives of list but didn't find

Dear R-helpers, I am developing a Mixed-Effects model for a study of immunoassays using 'lme4' library. The study design is as follows: 10 samples were run using 7 different immunoassays, 3 times each, in duplicates. Data attached. I have developed the following model: c.lme <- lmer(Result~SPL + (SPL|Assay/Run) -1, data=data) This model has excellent predictions - the Rsquared of