similar to: question about extreme value distribution

HI, this is my problem I want to subset this file df, using only unique df$exon printing the line once even if df$exon appear several times: unique(df$exon) will show me the unique exons If I try to print only the unique exon lines with df[unique(df$exon),] -this doesn't print only the unique ones :( could you help? thanks Nat exon size chr start

Hi R friends, I've been attempting to create plots with multiple alpha values using Cairo to save them on a windows (32b XP) platform as it doesn't support more than 3 alpha values. This worked well until I wanted a postscript file (unsupported) and as a attempted work around I installed RGtk2. So far so good, however now when I try to use a >CairoPDF("alpha.pdf", 6, 6,

Hi, I do not get how par works, help please. Let's say I have a simple plot: plot(1:10) I want to change the font size for the x axis... how do I do that? Thank you, Stephane -- The Wellcome Trust Sanger Institute is operated by Genome Research Limited, a charity registered in England with number 1021457 and a compa ny registered in England with number 2742969,

Hi everybody, I've been looking around, but can't seem to find the answer. I get a list of names (which vary, so I never know them in advance), and transform them into variables, each of which is a dataframe, using assign: polyList <- c("rs123", "rs124", "rs555", "rs000") numPoly <- length(polyList) for (k in 1:numPoly) {

Hello, The 179th and 180th elements of my list of lists look like this: [[179]] [[179]]$desc [1] ">ipi|IPI00646510|IPI00646510.2 ISOFORM P60-HCK OF TYROSINE-PROTEIN KINASE HCK." [[179]]$seq [1] "MGGRSSCEDPGCPRDEERAPRMGCMKSKFLQVGGNTFSKTETSASPHCPVYVPDPTSTIKPGPNSHNSNTP GIREAGSEDIIVVALYDYEAIHHEDLSFQKGDQMVVLEESGEWWKARSLATRKEGYIPSNYVARVDSLETEE

Hi - If I have a hypothetical distribution as an array distribution<-c(0,1,2,3,4,5,6,7,8,9) and I want to find the probability there is a value smaller than a new value. new_value<-4 (such that I'd get this type of output) new_value p-value 4 0.5 3.4 0.4 3 0.4 0 0.1 -1 0.0 Thanks for the help, I bet this is really easy... :/ Stephen -- The Wellcome Trust Sanger Institute is

Hi, I would like some extra information on the 'cgh' package in R. I noticed that there isn't much activity regarding this package on the R and BioC mailing list (I googled it). I started using this package and I have few questions: 1/ As I have a custom tiling like array @8um features resolution (affy), I have a lot of probes to work with. I'm assuming it is correct to

Hi, I'm trying to get a plot ready for publication, which involves getting it to look nice at a rather small size and to fine-tune positioning all the labels and sizes of the margins. I realise that I may not be doing this the right way and I welcome any comments about better approaches to do this. What I have done so far is open an X11 device with the size I want for the final output and I

To whom it may concern: My OS: CentOS 5.5 R version: R-2.11.1 My questions are as follows: First, I inputted commands: ----------------------------- # ./configure --enable-R-shlib # make # make check ----------------------------- Everything goes well until "make check". screen log: ************************************** make[1]: Entering directory `/tmp/R-2.11.1/tests' make[2]:

Hello, I need a specific version of R installed for consistency reasons. I do the standard setup steps: echo "deb https://cran.ma.imperial.ac.uk/bin/linux/ubuntu trusty/" | sudo tee -a /etc/apt/sources.list sudo apt-key adv --keyserver keyserver.ubuntu.com --recv-keys E084DAB9 sudo apt-get update And then a simple install call, which used to work just fine some time ago, tosses an

I am running dovecot 2.2.2 with tcp based replication, and experiencing some duplicated emails. `doveconf -n` output is below. I have narrowed it down to the following scenario: An email arrives, and is successfully replicated to both nodes. It is in INBOX/new/ at this point on both servers. Connect with a mail client, and delete the message - without delayed expunge. So, for example, mutt

I am running two servers with Dovecot v2.2.12 on CentOS x86_64 (5.10 and 6.5 respectively) and users are virtual over ldap. I have setup our main internal server (vmail.example.com) with dsync replication according to the first part of http://wiki2.dovecot.org/Replication. The second one (vmail1.example.com) will be the failover server which we want to be a real-time mirror (but can be

Hello List I have a strange problem here which i try to analyse, but i'm stuck. Maybe someone has a hint? What happened: A few weeks ago one of the LDAPS Servers which is not maintained by us has crashed. From that moment on, users could still login to check their emails, but they were not able to send any email through postfix (which uses smtpd_sasl_type = dovecot) What i do not

Hi - I am interested in solving variance components for the data below with respect to the response variable, Expression within R. However, the covariates aren't independent and they also have a class (of which the total variance explained by covariates in that class I am most interested in). Very naively, I have tried to look at each individual covariates variance like this >

First, apologies for no example data but I don't think it's needed in this case, Q: Can (and if so how ) the line end style be changed for 'cloud' plots? I've tried par(lend=2), trellis.par.set(add.line = list(lend=2)) and much googling but to no avail Thanks in advance Dan P.S. the reason for this is that the round end looks bad at lwd=3 or more Daniel Alcock Malaria

Hi I hope someone can shed some light on this: For some reason when I read.table("bfx.txt") R decides to only give back the first character from each column in each row as one single column. Like this: V1 1 ÿþr 2 \n 3 r 4 1 5 0 6 A 7 G 8 \n 9 r 10 1 11 0 12 T 13 C 14 \n The data should be:

hi, I have 2 files example 1 and example 2 and would like to know what is in example2 and not in example1 (attached) V1 contain data which could be in duplicated which I am using as identifiers I used setdiff(example2$V1,example1$V1) to find the identifiers which are specific to example2: [1] "rs2276598" "rs17253672" I am looking for a way to get an output with all

Hi, I'm using read.table in a loop, to read in multiple files. The problem is that when a file is missing there is an error message and the loop is broken; what I'd like to do is to test for the error and simply do "next" instead of breaking the loop. Anybody knows how to do that? Example: filelist <- c("file1.txt", "file2.txt",

Dear all, I have a question about the 'sdev' value returned by the princomp function (which does principal components analysis). On the help page for princomp it says 'sdev' is 'the standard deviations of the principal components'. However, when I calculate the principal components for the USArrests data set, I don't find this to be the case: Here is how I

Hello, I'm a newcomer to R so please forgive me if this is a silly question. It's that I have a linear regression: fm <- lm (x ~ y) and I want to test whether the slope of the regression is significantly less than 1. How can I do this in R? I'm also interested in comparing the slopes of two regressions: fm1 <- lm (x ~ y) fm2 <- lm (a ~ b) and asking if the slope of fm1 is