similar to: read.table mystery

Hello! I have this problem: I want to create a Venn's diagram with three lists of genes'names. The first is all the genome, the second a subset of it comprising all mitochondrial genes, and the third including all genes that correlate with a given gene. This is what I do: > library(gplots) > A<-read.delim("F:/.../mito genes just names.txt") >

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

I am a relative newbie to survival analysis and R in general, but would like to use the coxme package to analyse some data I currently have. The data is relative to survival times of drosophila melanogaster populations to infection with pathogens, and has the variables: Time, Status, Treatment (4 treatments + 2 controls) Population Replicate ?and I'm currently using the following call

Dear all, Here is my code which am using to combine 5th column from different data sets. Here is the function to do my job genesymbol.append.file <-NULL gene.column <- NULL readGeneSymbol <- function(files,genesymbol.column=5){ for(i in fnames){ temp <- read.table(fnames,header=T,sep="\t",stringsAsFactors=F,quote="\"")

Dear all, we seek a talented and motivated post-doc to develop computational methods for inferring the molecular basis of genetic diseases by integration of personal omics data. Research topics include: identifying causal mutations of rare disease patients by meta-analysis; inferring disease-causing molecular pathways from genotype, phenotypes, and omics profile of patient-derived cell lines

Dear all, I have a data which is in text file and i would like to import the data to R. From the manual, i?ve found the read.table command function is the most appropriate but when i wrote the command an error had occur. It say ?Error in read.table"C:/Users/user/Documents/cfa-1.txt", header = T, sep = "\t",skip=10) :more columns than column names?. Please help me with this as

Dear all, I am currently analyzing a set of arrays hybe on the lattest affy Drosophila 2.0 GeneChip. I am trying to run simple annaffy analysis but can¹t find what is the name of the annotation file I need to use. Here is the output of the AffyBatch object I am using: > expData AffyBatch object size of arrays=732x732 features (16749 kb) cdf=Drosophila_2 (18952 affyids) number of samples=4

Hello, I've read in a data file as a frame and now I'd like the columns to be rows and the rows to be columns. What's the easiest way to do this in R? > class(d) [1] "data.frame" > rownames(d) [1] "98x101" "40x98" "30x40" "0x30" > colnames(d) [1] "H..sapiens" "C..elegans"

Hello, I've just started R, and I'm getting a bit mad using it. I've managed to produce a barplot with the labels for the ytick marks placed horizontally (perpendicular to the y-axis) usiing par(las=1). The problem is that most of my labels are in part beyond the plotting area because they are rather long (e.g. "H. sapiens", "D. melanogaster" ...). What is the

Hi everybody, I have a problem when trying to do the quality control with the packages simpleaffy and affyQCReport with the drosophila chip 2.0 At first I got the messeage, that the *.qcdef file is not there. I followed the instructions in tha manual and created the file like that: array drosophila2cdf alpha1 0.05 alpha2 0.065 spk bioB AFFX-r2-Ec-bioB-3_at spk bioC AFFX-r2-Ec-bioC-3_at spk bioD

I am running R version ??? under Redhat 5.2. It seems as though the --nsize object has no effct on the size of the allocated Ncells as determined using gc(). Yes, I have that much data.... That is if I envoke R with R --vsize 100 --nsize 5000000 then type gc() I get free total Ncells 92202 200000 Vcells 12928414 13107200 Thanks Tony Long Ecology and Evolutionary Biology Steinhaus

Hi, everybody. I am working processing EEG data from 1000 pacients. I have a specific syntax to perform the Spectral Analysis and a loop to analyse all subjects. each subject data are in separate folders (P1, P2 P3...) My question is: in some cases, some errors can appear in one subject. I want to know if is possible to jump to the next subject and perform the same syntax , exibiting an error