similar to: OPENSSH Support Policies

Dear listers, After activating the name space for my bioconductor package (prada) I successfully ran R CMD check. However when loading the package in R and running the examples the imported function brewer.pal from package RColorBrewer is not found. I can directly call brewer.pal from the RColorBrewer name space typing RColorBrewer::brewer.pal, but it is not imported into my prada name space. When

Dear listers, After activating the name space for my bioconductor package (prada) I successfully ran R CMD check. However when loading the package in R and running the examples the imported function brewer.pal from package RColorBrewer is not found. I can directly call brewer.pal from the RColorBrewer name space typing RColorBrewer::brewer.pal, but it is not imported into my prada name space. When

Hi, I'm trying to create a heatmap with color ranges for different values in my matrix. For example: If x > 5 , use orange gradient if x < 1.5, use red gradient ..... Right now I have the following: orgPal<-brewer.pal(3,"Oranges") bluPal<-brewer.pal(3,"Blues") redPal<-brewer.pal(3,"Reds") grad <- ifelse(randMat >

Hello, I would like to reshape my data for presentation purposes from something like this: > test <- data.frame(a=c("A","A","A","A","B","B","B"),b=c(1,1,2,2,1,1,1),c=1:7) > test a b c 1 A 1 1 2 A 1 2 3 A 2 3 4 A 2 4 5 B 1 5 6 B 1 6 7 B 1 7 to something like this: a b c 1 A 1 1,2 3 A 2 3,4 5 B 1 5,6,7 This seems

Hello, Is there a simple way to get the class type for each column of a data.frame? I am in the situation where I would like to get all the columns of a data.frame that are factors. I have tried: apply(df,2,class) but all the columns come back as class "character". Thanks Dan -- ************************************************************** Daniel Brewer, Ph.D. Institute of

Hi, I am trying to produce a plot with an xlabel that reads (x >= 5%) with the >= turned into the correct epression. I can do this up to the percentage sign by specifiing xlab=expression(x>=5). Whatever I do to include the % sign as well doesn't seem to work. xlab=bquote(x>=5.("%")) almost works but includes brackets. Anyonw know how to solve this one Dan PS I am

Hello, I would like to scrape some html tables from a web page and convert them to a data.frame so I can perform further analysis. Could anyone tell me the best way to do this? Would it be more appropriate to use an external tool first? Thanks Dan -- ************************************************************** Daniel Brewer, Ph.D. Institute of Cancer Research Molecular Carcinogenesis

Hello, I have a data.frame like this: a 1,2,3,4 b1 b 6,7 b3 And I would like to transform it to this: a 1 b1 a 2 b1 a 3 b1 a 4 b1 b 6 b3 b 7 b3 I have been looking at ddply but can't seem to work it out. ANy help would be gratefully received. Dan -- ************************************************************** Daniel Brewer, Ph.D. Institute of Cancer Research Molecular Carcinogenesis

Hello, I would like to be able to calculate the age of someone at a particular date. Both dates are date objects. Here is what I have come up with: floor(as.numeric(sampleInfo$Date.of.DIAGNOSIS-sampleInfo$Date.of.birth)/365.25) Is this the best approach? or is there an inbuilt function? I have looked at difftime but that does not seem to allow output in years. Many thanks Dan --

Hi all, I have 35 years, I have been working since I exit from University. I have returned to make a phd...I´m fortunatelly I´m working so unfortunatelly I have few time to study (I have also a baby with 16 months). I´m trying to use R but sometimes it gets very difficult. I´m trying to put this code in my polygon function colramp=colorRampPalette(brewer.pal(9,"YlOrRd") to use

Hello, I have a data.frame which I would like to sort with the primary key decreasing while the secondry key is increasing e.g. x <- data.frame(One=c(1,1,1,2,2,3,4,5),Two=c(2,3,1,2,3,3,3,3)) I would like to order it so it looks like this: One Two 8 5 3 7 4 3 6 3 3 4 2 2 5 2 3 3 1 1 1 1 2 2 1 3 i.e. primarily decreasing in the 1st column but if there is a

Hi, I have a particular barplot I would like to generate, but I am having trouble getting it to work. What I would like is in effect two barplots with stacked bars merged into one. For example, I have two samples (yoda1,yoda2) on which I measure whether two variables (var1,var2) are present or absent for a number of measurements on that sample. > var1 <- data.frame(yoda1=c(3,7),

Hello, I have a gene expression matrix with columns being samples and rows being genes. I would like to display the expression values for each gene. I have two groups which I colour differently. The aim is to see if there is any difference between the two groups consistently across genes. So the following works well:

Hello, I would like to calculate the 95% confidence intervals for the estimates of a linear model and I just wanted to check that I am doing it correct. Is it just: Estimate + 1.95996*Std.Error to Estimate - 1.95996*Std.Error or is there another approach that doesn't assume a normal distrbution? Thanks. Apologies for my naiivity Dan --

Hello, I have some data that has dates in the form 27.02.37. I convert them to a date object as follows: as.Date(data$date,format="%d.%m.%y") But this gives me years such as 2037 when I would like them to be 1937. I thought of trying to take off some time i.e. as.Date(camCD$DoB,format="%d.%m.%y") - 100*365 But that doesn't seem to work out correctly. Any ideas how to

I am trying to display data ion top of colored rectangle on the same levelplot. While using superpose the shape and symbol for levelplot changes to small circle. I am unable to figure what exactly I am doing wrong... It may be straightforward solution... .I will appreciate help to fix this issue.

Does anyone have any ideas how I could do a power calculation for a log rank test. I would like to know what the suggested sample sizes would be to pick a difference when the control to active are in a ratio of 80% to 20%. Thanks Dan -- ************************************************************** Daniel Brewer, Ph.D. Institute of Cancer Research Email: daniel.brewer at icr.ac.uk

All: Below is my code for creating a basic horizontal, stacked barchart. I'd like to label the plot in two ways: 1) place the x values in each piece and 2) place the y values above each piece (angled). I'm currently using lattice, but I'm open to suggestions using ggplot2. Questions: 1. Can this be done?...I assume yes. So, what are the best options/functions for doing this. 2.

Hello, I am finding that the melt function from the reshape library causes errors when applied to a data.frame that contains numeric and character columns. For example, melt(id.vars="ID",data.frame(ID=1:3,date=c("a","b","c"),value=c(1,4,5))) ID variable value 1 1 date a 2 2 date b 3 3 date c 4 1 value <NA> 5 2

Hi, I have a script that takes a subset of genes on a microarray and tries to fit a coxph model to the expression values for each gene. This seems to work fine but in some cases it produces warnings and/or errors. For example: Error in fitter(X, Y, strats, offset, init, control, weights = weights, : NA/NaN/Inf in foreign function call (arg 6) In addition: Warning message: In fitter(X, Y,