similar to: list.dir() function

Hi Thanks Gabor for your suggestion. I am posting the code that worked for me. dataframe1 = data.frame(cbind(Src = c(1,1,1,2,3), Target1 = c('aaa','bbb','ccc','aaa','ddd'))); #must be data frame dataframe2 = data.frame(cbind(Src = c(2,3,4,4,4), Target2 = c('aaaa','dddd','bbbb','eeee','ffff'))); dataframe3 =

Hi R users, I try to read a data file (tab delimited format) in which some of the entries in a particular field are missing. Is it possible to fill the unavailable data with 'UnAvailable' string while performing read.table() Something like df = read.table(DataFile, header=FALSE, fill_missing_entry = 'unAvailable') 1 thanks -- ------------- Mary Kindall Yorktown Heights,

I have a value a=300 observation (x) = sample(1:50) How to find a p-value from this. I need to show that "a" is different fom mean(x). Thanks -- ------------- Mary Kindall Yorktown Heights, NY USA [[alternative HTML version deleted]]

Hi I have a data frame with values in following format. param case1 1 a 2 b 2 c 2 d 3 e 4 f how to use aggregate so that it I only one row for each 'param' value. the output for the above input should be param case1 1 a 2 b,c,d 3 e 4 f Thanks M -- ------------- Mary

> > I have two dataframes and want to perform cbind and then write into a > file. The number of entries are more than a million in both frames. R is > taking a lot of time performing this operation. > > Is there any alternate way to perform cbind? > > x = table1[1:1000000,1:4] > y = table2[1:1000000,3:6] > > z = cbind(x,y) //hanging the machine > >

In a frequency distribution table (bell shaped), how can we find the most frequent range? for example: x = c(1,2, 4,4,4,4, 5,5,5,6,6,5,5,5,5,5,6,6,6,13, 17,17,30,100,300) barplot(table(x)) In the code above, which function do we use to find that the most frequent value range from 4 to 6. Thanks. -- ------------- Mary Kindall Yorktown Heights, NY USA [[alternative HTML version deleted]]

I installed some downloaded packages in R. I always do $sudo R CMD INSTALL <anRpackage.tar.gz> By default it is storing these packages into my directory /home/mary/R/x86_64-pc-linux-gnu-library/2.13/. However I want them to be systemwide into /usr/local/lib/R/site-library/ folder. I tried $sudo R R> install.packages("anRpackage", dep=TRUE) I did not succeed into getting

I installed some downloaded packages in R. I always do $sudo R CMD INSTALL <anRpackage.tar.gz> By default it is storing these packages into my directory /home/mary/R/x86_64-pc-linux-gnu-library/2.13/. However I want them to be systemwide into /usr/local/lib/R/site-library/ folder. I tried $sudo R R> install.packages("anRpackage", dep=TRUE) I did not succeed into getting

Hi R users, I am new to R and am trying to merge data frames in the following way. Suppose I have n data frames each with two fields. Field 1 is common among data frames but may have different entries. Field 2 is different. Data frame 1: Src Target1 1 aaa 1 bbb 1 ccc 2 aaa 3 ddd Data frame 2: Src Target2 2 aaaa 3 dddd 4 bbbb 4

I am downloading say 100 files from ucsc website and storing it into dest folder. download.file function create a file in destination folder even if the file is not present which is something I dont want. So I wrote if condition to remove the file if the download function has non zero value. Now it exits when there is an error or file not present. How can I use "try" and "if"

I am planning to download a large number of files from some website. I am using the following script. files2down = c('aaa', 'bbb', ................) for (i in 1: len) { print(paste('downloading file', i, ' of total ', len)); url = paste(urlPrefix, files2down[i], sep='') destfile = paste (dest, 'inDir', files2down[i], sep='/' )

Hi I have a data frame in the following form. There are two groups and for each 'width' relative frequency for group1 and group2 is given. How to plot this in R using ggplot or other package. Width relativeFrequency1 relativeFrequency2 1 100 0.0006388783 0.02265428 2 200 0.0022677303 0.02948625 3 300 0.0061182673 0.01739936 4 400 0.0152237225 0.02569902 5 500

I have a list of txt files that I want to convert into .rdata R data object. filenames 1. "./file1.txt" 2. "./file2.txt" 3. "./file3.txt" 4. "./file4.txt" 5. "./file5.txt" 6. "./file6.txt" 7. "./file7.txt" 8. "./file8.txt" 9. "./file9.txt" 10. "./file10.txt" I saved these files as for ( i in

I have two one dimensional list of elements and want to perform cbind and then write into a file. The number of entries are more than a million in both lists. R is taking a lot of time performing this operation. Is there any alternate way to perform cbind? x = table1[1:1000000,1] y = table2[1:1000000,5] z = cbind(x,y) //hanging the machine write.table(z,'out.txt) -- -------------

Hello, Everyone, how to use "loop" to make the process automatic and fast? When compute each sample, the script type in R almost the same, just the input and output file's name is changed(chr1 change to chr2, chr3,chr4...). The first sample's script like this: >chr1=MEDIPS.readAlignedSeqences(BSgenome="hg19", file="chr1",numrows= )

Hi, I was wondering I'm going about this in the correct way. I need to test if there are coding sequences or exons in hg19 which match a string of 100bp "D" i.e. [A,G or T]. However I'm getting a strange result. I get a hit on chr7, using the 100bp search however when I search with 60bp sequence of "D" I don't get any hits. library("BSgenome")

Hi, I was installing hundreds of packages on a machine with a single call to install.packages() and after a long time the call to install.packages() finally returned with the following warnings and errors: Warning messages: 1: packages ?hgu133aprobe?, ?hgu95av2.db?, ?BSgenome.Celegans.UCSC.ce2?, ?BSgenome.Mmusculus.UCSC.mm10?, ?BSgenome.Dmelanogaster.UCSC.dm3.masked?,

Folks, The registration site for the Xen Summit is now online, please see details below. Although the official deadline for presentation proposals has passed its still possible to consider them -- please send a paragraph describing your proposal to summit-talks <at> xensource.com as soon as possible. Hope to see you at the summit! Best, Ian We are pleased to announce that IBM has

Hao Yu Commercial Scale Out IBM T.J. Watson Research Center 1101 Kitchawan Rd/Route 134, Rm 36-019 Yorktown Heights, NY 10598-0218 914-945-1854, Email: yuh@us.ibm.com ----- Forwarded by Hao Yu/Watson/IBM on 10/31/2006 09:52 AM ----- Jimi

Hao Yu Commercial Scale Out IBM T.J. Watson Research Center 1101 Kitchawan Rd/Route 134, Rm 36-019 Yorktown Heights, NY 10598-0218 914-945-1854, Email: yuh@us.ibm.com ----- Forwarded by Hao Yu/Watson/IBM on 10/31/2006 09:52 AM ----- Jimi