similar to: mhtplot?

Hello All, I am drawing a graph having 18 small graphs inside using par(mfrow = c(6,3)) command. My problem is how to specify the margins of the whole 18 graphs. I used par(mar=c(6.5, 6.5, 1.5, 1.5)) for each graph separately already but it does not left any margins for the 'mtext()' for the margins of the whole 18 graphs. Any comments please. Thanks, MJK [[alternative HTML version

Hello all, I am trying to make a genome association plot for p-values related to SNPs and was fortunate to find that R contains a package that produces Manhattan plots which is what's preferred for my current project. The function mhtplot() is found in the 'gap' package which I installed in R 2.11.1 on Windows. I thought I'd test out the function first with the examples they

Hello all, I have this simple barplot code: ifn <- "id.dat" dat <- read.table(ifn) ofn <- "id.png" bitmap(ofn, type = "png256", width = 30, height = 30, pointsize = 30, bg = "white",res=50) par(mar=c(5, 5, 3, 2),lwd=5) par(cex.main=1.6,cex.lab=1.6,cex.axis=1.6) names(dat)<-c("NumberOfPeople","Average")

Hi R-ers, I am struggling with my x-axis in a association plot. What I would like is to place the labels of the x-axis between the tick markers and normally the labels are printed at the place where the tick marker is placed. I don???t want to move the tick marker (it gives the switch between one chromosome and the next) but I just want to put the chromosome number in between the

Hello All: If I do have: x = (2, 4, 5, 5, 6, 4, 5, 2, 1) y = (9, 11.5, 12.5, 13, 14, 19, 20, 21, 22) I wanted to find a simple function in R which calculates the averages of X and Y for every 4 unit increase of Y. The results should look like: x = (4, 6, 3) #where (2+ 4+ 5+ 5)/ 4 = 4 and ... y = (12, 14, 20.5) #where (9+11.5+12.5+13)/ 4 = 12 and ... Thanks, Mohsen

Dear all, I want to design a histogram and I need to have the frequency at certain points. For example I have the following 2 columns: *X Y* 0.1 25 0.4 22 0.45 11 0.55 21 I want the chart to have 4 columns. First column is from 0.0-0.1 (on X) and frequency is 25. Next colum is wider and form 0.1-0.4 with 22 frequency. Next column is narrow with 11 frequency and the last column is

Hello All: I am having the following data file named "data.dat" Number of people Number of pets Age of trees ifn <- "data.dat" dat <- read.table(ifn) colnames(dat)<-c("Variables") I want R to read all these in a string but when I ask R to read these, it gives me error because there more than one word in each line. Any comments to solve this problem?

Hello All: Can anybody tell me what is the problem with my program please. I have an error message as appears below. My program is: ifn <- "Jul08_09.LM" data <- read.table(ifn) ofn <- "Jul.png" bitmap(ofn, type = "png256", width = 30, height = 30, pointsize = 30, bg = "white",res=50) par(mar=c(5, 5, 3, 2),lwd=5)

Hello everyone, I am using the following program to get the p-value of some numbers (column 'LR' of the data.dat file). I want to write the 1st and 2nd column of the output file (data.out) as an integer while the program change them to double. Could anybody please tell me how I can write the code which writes the values of the first two columns as integer? Thanks library ('MASS')

Hello All: I am having the following error message when I increase the resolution ("res=30" to "res=350") of my graph: Error in plot.xy(xy, type, ...) : ignoring SIGPIPE signal Calls: plot -> plot.default -> plot.xy Execution halted Here is the code: ofn <- "MAY11.png" bitmap(ofn, type = "png256", width = 21, height = 30, pointsize = 30, bg =

Hello everyone, I have this problem, when I launch findFn from RStudio ( but not from R console) the table of results does not pop in the browser and I got a warning. findFn("test") found 9914 matches; retrieving 20 pages, 400 matches. 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 Warning message: In parseHTML(href) : NAs introduced by coercion Thanks H. [[alternative HTML

Hello all, I am doing a Likelihood Ratio (LR) test in my simulation and I have a vector LR values (each with 1 degree of freedom) at the end of my simulation. Can anybody tell me how I can write a 'R' code which gives me the p-value for each of those LR values. Thanks [[alternative HTML version deleted]]

Issue with sos package I am now using R 2.14 and the sos package did work with 2.13 > library(sos) Loading required package: brew Attaching package: 'sos' The following object(s) are masked from 'package:JGR': installPackages The following object(s) are masked from 'package:utils': > library("sos"); findFn("random forest") found 253

Hello, I have used the following data to draw my barplot: BL LR Q 36.35 1.00 1.92 36.91 4.00 0.00 25.70 6.00 0.00 34.38 3.00 1.92 05.32 0.50 0.00 BL<-c(36.35, 36.91, 25.70, 34.38, 05.32) LR<-c(1.00, 4.00, 6.00, 3.00, 0.50) Q<-<(1.92, 0.00, 0.00, 1.92, 0.00) barplot(dt$LR, main='LR Value', col='orange', border='black',

Hi all, We know that Hermite polynomial is for Gaussian, Laguerre polynomial for Exponential distribution, Legendre polynomial for uniform distribution, Jacobi polynomial for Beta distribution. Does anyone know which kind of polynomial deals with the log-normal, Student抯 t, Inverse gamma and Fisher抯 F distribution? Thank you in advance! David [[alternative HTML version deleted]]

Hello everyone, I am using the following program to get the p-value of some numbers (column 'LR' of the data.dat file). I want to write the 1st and 2nd column of the output file (data.out) as an integer while the program change them. Could anybody please tell me how I can write the code which writes the values of the first two columns as integer? For your convenience, I have attached the

All Is there a package or library that will, given a nucleotide sequence 1. calculate the extinction coefficient at 260 nm for (Beer-Lambert's law) 2. calculate molecular weight 3. return it's complementary sequence I was able to find several packages that can do similar calculations for an amino acid sequence for proteins but none for nucleic acids. Any pointers, etc. would be

I have the data set like this: A B C D E A 2 235 327 278 2569 B 244 5 122 248 259 C 324 2235 13 98 69 D 22 23 2 7 6569 E 23 2 2 278 9 This is the interaction data between A, B, C, D, and E. How to construct a tree including A-E? -- View

library(sos) (mp <- findFn('{molecular properties}')) ????? ** found 7 matches in 4 packages and opened two web pages in my default browser with (a) the 7 matches and (b) the 4 packages. The first function was something for amino acids, like you suggested.? Two others returned compound and substance information from PubChem. ????? Does this help? ????? Spencer On

Don't do anything yet. I may have found the problem by accident. I tried to use the computer from something else, and it was being drastically slowed down by some leftover processes, which turned out to be xlhtml. That is something that converts Excel files. Apparently, some excel files got into the libraries, and they were causing the indexing to hang completely. I am now running everything