similar to: keep track of selected observations over time

hello every one, How to function more than one loop in R? I have the following problem to be solved with the a method of three loops, can you help me please? The data is attached with this message. The data is composed of two parts, cleaved (denoted by ?cleaved?) and non cleaved (denoted by ?noncleaved?). ? to access to the ith peptide, you can use X$Peptide[i] ? to access to the ith label,

Hello everyone, I would like to parse very large xml files from MS/MS experiments and create R objects from their content. (By very large, I mean going up to 5-10Gb, although I am using a 'small' 40M file to test my code.) My first attempt at parsing the 40M file, using the XML package, took more than 2200 seconds and left me quite disappointed. I managed to cut that down to around 40

Hi, I'm fairly new to R but I have a large dataset (300000 obs) containing patient material. Some patients came 2-9 times during the three year observation period. The patients are identified by a unique idnr, the sessions can be distinguished using the session date. How can I recode the date of the session to a session id (1-9). This would be necessary to obtain information and do some

Dear R-users! I am using R 1.7.0, under Windows XP. Having some hospital discharge data (admission date and discharge date for each patient), I want to get the number of patients in the hospital on a given date. My data look like (simple example): > x <- data.frame(patid=c("pat1", "pat2"), adm.date = c("15.03.2002", "16.03.2002"),

the content of th character vector (of length 1) is as follows: a <- "something2 ....pat1 name1 pat2 something2....pat1 name2 pat2....pat1 name3 pat2 " I would like to extract the character bewteen pat1 and pat2. That's to say, I would like to get a vecter of c("name1", "name2","name3"). What I did is use strsplit() twise. But I wonder if there

How can I delete both rows and columns that do not meet a particular cut off value. Example: > d <- rbind(c(0, 1, 6, 4), + c(2, 5, 7, 5), + c(3, 6, 1, 6), + c(4, 4, 4, 4)) > f <- as.matrix(d) > f [,1] [,2] [,3] [,4] [1,] 0 1 6 4 [2,] 2 5 7 5 [3,] 3 6 1 6 [4,] 4 4 4 4 I would like to delete all

Hi !  Facing problem with " getSequence" commend .  when only biomaRt package loaded the following example working well  >mart <- useMart("ensembl",dataset="hsapiens_gene_ensembl") >seq = getSequence(id="BRCA1", type="hgnc_symbol", seqType="peptide", mart = mart) show(seq) but when i have loaded the seqinr, i got problem

I have following codes for ggplots. The legends are given in the plot do not match with the values specified in the codes given below. Your helps highly appreciated. Greg library(ggplot2) p <- ggplot(a,aes(x=NO_BMI_FI_beta ,y=FI_beta ,color= Super.Pathway))+ theme_bw() +theme(panel.border=element_blank()) + geom_point(size=3) p2<-p+scale_color_manual(name="Super.Pathway",

Hi there, I am a beginner. I would like to change the error bars in the bargraph.CI function from the default (se) to (sd). The help file says ci.fun= function(x) c(fun(x)-se(x), fun(x)+se(x)) Is there a simple way of telling the function what (x) precisely is - I already define in in the of the bargraph.CI function and assume that is should be able to use that information. cheers, Herwig --

Hello, Everyone, I am a student an a new learner of R and I am trying to do my homework in R. I have 10 files need to be read and process seperately. I really want to write the codes into something like "macro" to save the lines instead of repeating 10 times of similar work. The following is part of my codes and I only extracted three lines for each repeating section. data.1 <-

I have been trying to download the climstats package: https://r-forge.r-project.org/R/?group_id=861 but it doesn't seem to run on R 3.0.2 or 3.0.1 and the zipfile is empty. Does anyone know the status of this package or where I can download it. Thanks ****************** Jenny Williams Spatial Information Scientist, GIS Unit Herbarium, Library, Art & Archives Directorate Royal Botanic

Loris, thank you thank you so much!!! do you know which CRAN mirror i select? There are loads - all related to various countries.. so appreciate the help Jenny Jenny Coetzee (PhD Candidate) Head of Prevention in Key Populations & Head of the Soweto Sex Worker Project Perinatal HIV Research Unit (PHRU) New Nurses Home Chris Hani Baragwanath Academic Hospital Soweto B: +27 11 989 9821 M:

Dear R-help, I have a loop, which is set to take about 26 hours to run at the rate it's going - this is ridiculous and I really need your help to find a more efficient way of loading up my array gpcc.array: #My data is stored in a table format with all the data in one long column #running though every longitute, for every latitude, for every year. The #original data is sotred as

Dear R-help community, I need to find out how to make the black axis line around my plot thicker in levelplot, I cannot find the correct command in ?levelplot or ?par - or if it is there I cannot get it to work - any help would be more that appreciated! Here's an example script - I would want the code for making the axis bordering the plot thicker: x<- seq(pi/4, 5*pi,length=100)

Here is the situation that I would like to create: Call comes in Receptionist sees that the caller ID is Jenny <8675309> Receptionist picks up phone and transfers call to Batman Batman looks at his phone and sees that the caller ID is Jenny <8675309> I can't seem to figure out how to forward the caller ID. Is this possible with Asterisk?

Hi Jenny, Jenny Coetzee <coetzeeJ at phru.co.za> writes: > Hello > > I am new to R, and am running a programme called RDS Analyst which > uses R. When trying to open RDS i am getting the following message: > JGR not installed. Attempting to install from CRAN? > > I am not sure if the issue is a Java one or an R one? > > I am working off a Mac using High Sierra

Dear R-help, Which random number generator function would you recommend for simply picking 15 random numbers from the sequence 0-42? I want to use replacement (so that the same number could potentially be picked more than once). I have read the R-help archives and the statistics and computing book on modern Applied statistics with S but the advice seems to be for much form complicated

Dear R-community, I am looking for some simple advice - I have a matrix (therefore 2 dimensional) of global temperature. Having read R-help I think that when I ask R to image() or levelplot() my matrix will it actually appear upside down - I think I therefore need to use the line: > levelplot(temperature.matrix[,ncol(output.temp):1], ........) to get it looking like it was on the globe due

Dear R-help community I am trying to write an R object (data.out) to a file in order to re-access it later and not have to re-load up the array with data every time. Here is the form of data.out > data.out <- list(lats=seq(88.542, -88.542, length=94), lons=seq(0, 360-1.875, length=192), date=vector(length=nyr*12), data=array(NA, c(nyr*12, 94*192)) ) I tried >

Dear R-help community, Does anybody know of a stats function in R that tells you which distribution best fits your data? I have tried look through the archives but have only found functions that tell you if it's normal or log etc. specifically - I am looking for a function that tells you (given a timeseries) what the distribution is. Any help/advice will be greatly appreciated, All the