Re: [R] QQ plot
Hi Michael, Thank you so much for that valuable idea! I will try first to clump or remove SNPs in LD and maybe the situation would improve. But this procedure of yours is definitely something that would come handy in future! Cheers, Ana On Wed, Nov 13, 2019 at 5:47 AM Michael Dewey wrote: > > Dear Ana > > As others have commented this is getting a bit off-topic but here are > some hints. > > It is helpful to distinguish two sorts of plot: archival plots and > impact plots. If you want to have an impact plot which gives you a > picture but possibly at the cost of completeness and accuracy then why not: > > 1 - plot a sample of your 5 million drawn at random > 2 - bin the data and plot median p-value against median expected > 3 - deal with overlap by choosing a graphical device which supports > transparency and plot points in very light grey so the overlap is more > visible. > > Michael > > On 12/11/2019 22:04, Ana Marija wrote: > > why I selected only those with P<0.003 to put on QQ plot is because > > the original data set contains 5556249 points and when I extract only > > P<0.001 I am getting 3713 points. Is there is a way to plot the whole > > data set, or choose only the representative points? > > > > On Tue, Nov 12, 2019 at 3:42 PM Ana Marija > > wrote: > >> > >> the smallest p value in my dataset goes to 9.89e-08. How do I make > >> that known on the new QQ plot with multiplied with 1000 values > >> > >> On Tue, Nov 12, 2019 at 3:37 PM Ana Marija > >> wrote: > >>> > >>> Just do I need to change the axis when I multiply with 1000 and what > >>> should I put on my axis? > >>> > >>> On Tue, Nov 12, 2019 at 3:07 PM Ana Marija > >>> wrote: > > Hi Duncan, > > yes I choose for QQ plot only P<1e-3 and multiplying everything with > 1000 works great! > This should not in my understanding influence the interpretation of > the plot, it is only changing the scale of axis. > > Thank you so much, > Ana > > On Tue, Nov 12, 2019 at 2:51 PM Duncan Murdoch > wrote: > > > > On 12/11/2019 2:56 p.m., Jim Lemon wrote: > >> I thought about this and did a little study of GWAS and the use of > >> p-values to assess significant associations. As Ana's plot begins at > >> values of about 0.001, this seems to imply that almost everything in > >> the genome is associated to some degree. One expects that most SNPs > >> will not be associated with a particular condition (p~1), so perhaps > >> something is going wrong in the calculations that produce the > >> p-values. > > > > I may be misunderstanding your last sentence, but if there is no > > association, the p-value would usually have a uniform distribution from > > 0 to 1, it wouldn't be near 1. > > > > I'd guess we're not seeing the p values from every test, only those that > > are less than 0.001. If that's true, and there are no effects, it makes > > sense to multiply all of them by 1000 to get U(0,1) values. On the > > plot, that would correspond to subtracting 3 from -log10(p), or adding 3 > > to the reference line, as Ana requested. > > > > Or just multiply them by 1000 and pass them to qq(): > > > > qq(dd$P*1000, main = "Q-Q plot of small GWAS p-values") > > > > As far as I can see, there's no way to tell qqman::qq to move the > > reference line. > > > > Duncan Murdoch > > > >> > >> Jim > >> > >> On Wed, Nov 13, 2019 at 12:28 AM Patrick (Malone Quantitative) > >> wrote: > >>> > >>> I agree with Abby. That would defeat the purpose of a QQ plot. > >>> > >>> On Mon, Nov 11, 2019, 9:54 PM Abby Spurdle > >>> wrote: > >>> > Hi > > I'm not familiar with the qqman package, or GWAS studies. > However, my guess would be that you're *not* supposed to change the > position of the line. > > On Tue, Nov 12, 2019 at 11:48 AM Ana Marija > > wrote: > > > > Hi, > > > > I was using this library, qqman > > https://cran.r-project.org/web/packages/qqman/vignettes/qqman.html > > > > to create QQ plot, attached. How would I change this default abline > > to > > start from the beginning of my QQ line? > > > > This is my code: > > qq(dd$P, main = "Q-Q plot of GWAS p-values") > > > > Thanks > > Ana > > __ > > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > > https://stat.ethz.ch/mailman/listinfo/r-help > > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > > and provide commented, minimal, self-contained, reproducible code. > > __ >
Re: [R] data load from excel files
Hello,
Try which.max?
Hope this helps,
Rui Barradas
Às 14:10 de 13/11/19, ani jaya escreveu:
Thank you very much Mr. Rui, but for delete the duplicated row I use:
...
library(tidyverse)
alldata<-data.frame(Reduce(rbind, pon1))
c<-(which(duplicated(alldata$Tanggal))) #duplicate
alldata<-alldata[-c,]
attach(alldata)
because not every last row from every df is bad one.
Another problem is I want to know when the max value is occurred. So
basically I have maximum value every month (maxi, n=360, from 1986 to
2015) and I want to find annual_maxima.
...
maxi<-lapply(pon1, function(x) max(x$RR,na.rm=T))
maxi<-data.frame(Reduce(rbind, maxi))
names(maxi)<-"maxi"
annual_maxima <- rep(NA,30)
date <- rep(NA,30)
for(i in 1:30){
annual_maxima[i] <- max(maxi$maxi[(i*12-11):(i*12)])
date[i]<-Tanggal[which(RR==annual_maxima[i])]
}
Here "alldata" contain "Tanggal" in this case is date and rainfall
("RR"). What I get is error stated that:
In date[i] <- Tanggal[which(RR == annual_maxima[i])] : number of items
to replace is not a multiple of replacement length
Maybe you have some idea where the problem is, I would be thankful.
Best,
Ani
On Wed, Nov 13, 2019 at 5:49 PM Rui Barradas > wrote:
Hello,
Maybe the following will get you close to what you want.
# remove the last row from every df
pon1 <- lapply(pon1, function(DF){
DF[[1]] <- as.Date(DF[["Tanggal"]], "%d-%m-%Y")
DF[-nrow(DF), ]
})
# order the list by year-month
inx_ym <- sapply(pon1, function(DF){
format(DF[["Tanggal"]][1], "%Y-%m")
})
pon1 <- pon1[order(inx_ym)]
# get the minimum and maximum of every "RR"
min.RR <- sapply(pon1, function(DF) min(DF[["RR"]], na.rm = TRUE))
max.RR <- sapply(pon1, function(DF) max(DF[["RR"]], na.rm = TRUE))
Hope this helps,
Rui Barradas
Às 07:50 de 13/11/19, ani jaya escreveu:
> Dear R-Help,
>
> I have 30 of year-based excel files and each file contain month
sheets. I
> have some problem here. My data is daily rainfall but there is
extra 1 day
> (first date of next month) for several sheets. My main goal is to
get the
> minimum value for every month.
>
> First, how to extract those data to list of data frame based on
year and
> delete every overlapping date?
> Second, how to sort it based on date with ascending order (old to
new)?
> Third, how to get the maximum together with the date?
>
> I did this one,
>
> ...
> file.list <- list.files(pattern='*.xlsx')
> file.list<-mixedsort(file.list)
>
> #
>
https://stackoverflow.com/questions/12945687/read-all-worksheets-in-an-excel-workbook-into-an-r-list-with-data-frames
>
> read_excel_allsheets <- function(filename, tibble = FALSE) {
> sheets <- readxl::excel_sheets(filename)
> x <- lapply(sheets, function(X) read.xlsx(filename, sheet=X,
rows=9:40,
> cols=1:2))
> if(!tibble) x <- lapply(x, as.data.frame)
> names(x) <- sheets
> x
> }
>
> pon<-lapply(file.list, function(i) read_excel_allsheets(i))
> pon1<-do.call("rbind",pon)
> names(pon1) <- paste0("M.", 1:360)
> pon1 <-lapply(pon1,function(x){x$RR[x$RR==] <- NA; x})
> pon1 <-lapply(pon1,function(x){x$RR[x$RR==""] <- NA; x})
> maxi<-lapply(pon1, function(x) max(x$RR,na.rm=T))
> maxi<-data.frame(Reduce(rbind, maxi))
> names(maxi)<-"maxi"
>
>
> but the list start from January for every year, and move to
February and so
> on. And there is no date in "maxi". Here some sample what I get
from my
> simple code.
>
>> pon1[256:258]$M.256
> Tanggal RR
> 1 01-09-2001 5.2
> 2 02-09-2001 0.3
> 3 03-09-2001 29.0
> 4 04-09-2001 0.7
> 5 05-09-2001 9.6
> 6 06-09-2001 0.7
> 7 07-09-2001 NA
> 8 08-09-2001 13.2
> 9 09-09-2001 NA
> 10 10-09-2001 NA
> 11 11-09-2001 0.0
> 12 12-09-2001 66.0
> 13 13-09-2001 0.0
> 14 14-09-2001 57.6
> 15 15-09-2001 18.0
> 16 16-09-2001 29.2
> 17 17-09-2001 52.2
> 18 18-09-2001 7.0
> 19 19-09-2001 NA
> 20 20-09-2001 74.5
> 21 21-09-2001 20.3
> 22 22-09-2001 49.6
> 23 23-09-2001 0.0
> 24 24-09-2001 1.3
> 25 25-09-2001 0.0
> 26 26-09-2001 1.0
> 27 27-09-2001 0.1
> 28 28-09-2001 1.9
> 29 29-09-2001 9.5
> 30 30-09-2001 3.3
> 31 01-10-2001 0.0
>
> $M.257
> Tanggal RR
> 1 01-09-2002 0.0
> 2 02-09-2002 0.0
> 3 03-09-2002 0.0
> 4 04-09-2002 12.8
> 5 05-09-2002 1.0
> 6 06-09-2002 0.0
> 7 07-09-2002 NA
> 8 08-09-2002 22.2
> 9 09-09-2002 NA
> 10 10-09-2002 NA
> 11 11-09-2002 0.0
> 12 12-09-2002 0.0
Re: [R] data load from excel files
Thank you very much Mr. Rui, but for delete the duplicated row I use:
...
library(tidyverse)
alldata<-data.frame(Reduce(rbind, pon1))
c<-(which(duplicated(alldata$Tanggal))) #duplicate
alldata<-alldata[-c,]
attach(alldata)
because not every last row from every df is bad one.
Another problem is I want to know when the max value is occurred. So
basically I have maximum value every month (maxi, n=360, from 1986 to 2015)
and I want to find annual_maxima.
...
maxi<-lapply(pon1, function(x) max(x$RR,na.rm=T))
maxi<-data.frame(Reduce(rbind, maxi))
names(maxi)<-"maxi"
annual_maxima <- rep(NA,30)
date <- rep(NA,30)
for(i in 1:30){
annual_maxima[i] <- max(maxi$maxi[(i*12-11):(i*12)])
date[i]<-Tanggal[which(RR==annual_maxima[i])]
}
Here "alldata" contain "Tanggal" in this case is date and rainfall ("RR").
What I get is error stated that:
In date[i] <- Tanggal[which(RR == annual_maxima[i])] : number of
items to replace is not a multiple of replacement length
Maybe you have some idea where the problem is, I would be thankful.
Best,
Ani
On Wed, Nov 13, 2019 at 5:49 PM Rui Barradas wrote:
> Hello,
>
> Maybe the following will get you close to what you want.
>
>
> # remove the last row from every df
> pon1 <- lapply(pon1, function(DF){
>DF[[1]] <- as.Date(DF[["Tanggal"]], "%d-%m-%Y")
>DF[-nrow(DF), ]
> })
>
>
> # order the list by year-month
> inx_ym <- sapply(pon1, function(DF){
>format(DF[["Tanggal"]][1], "%Y-%m")
> })
> pon1 <- pon1[order(inx_ym)]
>
>
> # get the minimum and maximum of every "RR"
> min.RR <- sapply(pon1, function(DF) min(DF[["RR"]], na.rm = TRUE))
> max.RR <- sapply(pon1, function(DF) max(DF[["RR"]], na.rm = TRUE))
>
>
> Hope this helps,
>
> Rui Barradas
>
>
>
> Às 07:50 de 13/11/19, ani jaya escreveu:
> > Dear R-Help,
> >
> > I have 30 of year-based excel files and each file contain month sheets. I
> > have some problem here. My data is daily rainfall but there is extra 1
> day
> > (first date of next month) for several sheets. My main goal is to get the
> > minimum value for every month.
> >
> > First, how to extract those data to list of data frame based on year and
> > delete every overlapping date?
> > Second, how to sort it based on date with ascending order (old to new)?
> > Third, how to get the maximum together with the date?
> >
> > I did this one,
> >
> > ...
> > file.list <- list.files(pattern='*.xlsx')
> > file.list<-mixedsort(file.list)
> >
> > #
> >
> https://stackoverflow.com/questions/12945687/read-all-worksheets-in-an-excel-workbook-into-an-r-list-with-data-frames
> >
> > read_excel_allsheets <- function(filename, tibble = FALSE) {
> >sheets <- readxl::excel_sheets(filename)
> >x <- lapply(sheets, function(X) read.xlsx(filename, sheet=X,
> rows=9:40,
> > cols=1:2))
> >if(!tibble) x <- lapply(x, as.data.frame)
> >names(x) <- sheets
> >x
> > }
> >
> > pon<-lapply(file.list, function(i) read_excel_allsheets(i))
> > pon1<-do.call("rbind",pon)
> > names(pon1) <- paste0("M.", 1:360)
> > pon1 <-lapply(pon1,function(x){x$RR[x$RR==] <- NA; x})
> > pon1 <-lapply(pon1,function(x){x$RR[x$RR==""] <- NA; x})
> > maxi<-lapply(pon1, function(x) max(x$RR,na.rm=T))
> > maxi<-data.frame(Reduce(rbind, maxi))
> > names(maxi)<-"maxi"
> >
> >
> > but the list start from January for every year, and move to February and
> so
> > on. And there is no date in "maxi". Here some sample what I get from my
> > simple code.
> >
> >> pon1[256:258]$M.256
> >Tanggal RR
> > 1 01-09-2001 5.2
> > 2 02-09-2001 0.3
> > 3 03-09-2001 29.0
> > 4 04-09-2001 0.7
> > 5 05-09-2001 9.6
> > 6 06-09-2001 0.7
> > 7 07-09-2001 NA
> > 8 08-09-2001 13.2
> > 9 09-09-2001 NA
> > 10 10-09-2001 NA
> > 11 11-09-2001 0.0
> > 12 12-09-2001 66.0
> > 13 13-09-2001 0.0
> > 14 14-09-2001 57.6
> > 15 15-09-2001 18.0
> > 16 16-09-2001 29.2
> > 17 17-09-2001 52.2
> > 18 18-09-2001 7.0
> > 19 19-09-2001 NA
> > 20 20-09-2001 74.5
> > 21 21-09-2001 20.3
> > 22 22-09-2001 49.6
> > 23 23-09-2001 0.0
> > 24 24-09-2001 1.3
> > 25 25-09-2001 0.0
> > 26 26-09-2001 1.0
> > 27 27-09-2001 0.1
> > 28 28-09-2001 1.9
> > 29 29-09-2001 9.5
> > 30 30-09-2001 3.3
> > 31 01-10-2001 0.0
> >
> > $M.257
> >Tanggal RR
> > 1 01-09-2002 0.0
> > 2 02-09-2002 0.0
> > 3 03-09-2002 0.0
> > 4 04-09-2002 12.8
> > 5 05-09-2002 1.0
> > 6 06-09-2002 0.0
> > 7 07-09-2002 NA
> > 8 08-09-2002 22.2
> > 9 09-09-2002 NA
> > 10 10-09-2002 NA
> > 11 11-09-2002 0.0
> > 12 12-09-2002 0.0
> > 13 13-09-2002 0.0
> > 14 14-09-2002 NA
> > 15 15-09-2002 0.0
> > 16 16-09-2002 0.0
> > 17 17-09-2002 0.0
> > 18 18-09-2002 13.3
> > 19 19-09-2002 0.0
> > 20 20-09-2002 0.0
> > 21 21-09-2002 0.0
> > 22 22-09-2002 0.0
> > 23 23-09-2002 0.0
> > 24 24-09-2002 0.0
> > 25 25-09-2002 0.0
> > 26 26-09-2002 0.5
> > 27 27-09-2002 2.1
> > 28 28-09-2002 NA
> > 29 29-09-2002 18.5
> > 30 30-09-2002 0.0
> > 31 01-10-2002 NA
> >
> > $M.258
> >
[R] [Rd] R 3.6.2 scheduled for December 12
Full schedule is available on developer.r-project.org. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Office: A 4.23 Email: pd@cbs.dk Priv: pda...@gmail.com __ r-de...@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel ___ r-annou...@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-announce __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] QQ plot
Dear Ana As others have commented this is getting a bit off-topic but here are some hints. It is helpful to distinguish two sorts of plot: archival plots and impact plots. If you want to have an impact plot which gives you a picture but possibly at the cost of completeness and accuracy then why not: 1 - plot a sample of your 5 million drawn at random 2 - bin the data and plot median p-value against median expected 3 - deal with overlap by choosing a graphical device which supports transparency and plot points in very light grey so the overlap is more visible. Michael On 12/11/2019 22:04, Ana Marija wrote: why I selected only those with P<0.003 to put on QQ plot is because the original data set contains 5556249 points and when I extract only P<0.001 I am getting 3713 points. Is there is a way to plot the whole data set, or choose only the representative points? On Tue, Nov 12, 2019 at 3:42 PM Ana Marija wrote: the smallest p value in my dataset goes to 9.89e-08. How do I make that known on the new QQ plot with multiplied with 1000 values On Tue, Nov 12, 2019 at 3:37 PM Ana Marija wrote: Just do I need to change the axis when I multiply with 1000 and what should I put on my axis? On Tue, Nov 12, 2019 at 3:07 PM Ana Marija wrote: Hi Duncan, yes I choose for QQ plot only P<1e-3 and multiplying everything with 1000 works great! This should not in my understanding influence the interpretation of the plot, it is only changing the scale of axis. Thank you so much, Ana On Tue, Nov 12, 2019 at 2:51 PM Duncan Murdoch wrote: On 12/11/2019 2:56 p.m., Jim Lemon wrote: I thought about this and did a little study of GWAS and the use of p-values to assess significant associations. As Ana's plot begins at values of about 0.001, this seems to imply that almost everything in the genome is associated to some degree. One expects that most SNPs will not be associated with a particular condition (p~1), so perhaps something is going wrong in the calculations that produce the p-values. I may be misunderstanding your last sentence, but if there is no association, the p-value would usually have a uniform distribution from 0 to 1, it wouldn't be near 1. I'd guess we're not seeing the p values from every test, only those that are less than 0.001. If that's true, and there are no effects, it makes sense to multiply all of them by 1000 to get U(0,1) values. On the plot, that would correspond to subtracting 3 from -log10(p), or adding 3 to the reference line, as Ana requested. Or just multiply them by 1000 and pass them to qq(): qq(dd$P*1000, main = "Q-Q plot of small GWAS p-values") As far as I can see, there's no way to tell qqman::qq to move the reference line. Duncan Murdoch Jim On Wed, Nov 13, 2019 at 12:28 AM Patrick (Malone Quantitative) wrote: I agree with Abby. That would defeat the purpose of a QQ plot. On Mon, Nov 11, 2019, 9:54 PM Abby Spurdle wrote: Hi I'm not familiar with the qqman package, or GWAS studies. However, my guess would be that you're *not* supposed to change the position of the line. On Tue, Nov 12, 2019 at 11:48 AM Ana Marija wrote: Hi, I was using this library, qqman https://cran.r-project.org/web/packages/qqman/vignettes/qqman.html to create QQ plot, attached. How would I change this default abline to start from the beginning of my QQ line? This is my code: qq(dd$P, main = "Q-Q plot of GWAS p-values") Thanks Ana __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained,
Re: [R] data load from excel files
Hello,
Maybe the following will get you close to what you want.
# remove the last row from every df
pon1 <- lapply(pon1, function(DF){
DF[[1]] <- as.Date(DF[["Tanggal"]], "%d-%m-%Y")
DF[-nrow(DF), ]
})
# order the list by year-month
inx_ym <- sapply(pon1, function(DF){
format(DF[["Tanggal"]][1], "%Y-%m")
})
pon1 <- pon1[order(inx_ym)]
# get the minimum and maximum of every "RR"
min.RR <- sapply(pon1, function(DF) min(DF[["RR"]], na.rm = TRUE))
max.RR <- sapply(pon1, function(DF) max(DF[["RR"]], na.rm = TRUE))
Hope this helps,
Rui Barradas
Às 07:50 de 13/11/19, ani jaya escreveu:
Dear R-Help,
I have 30 of year-based excel files and each file contain month sheets. I
have some problem here. My data is daily rainfall but there is extra 1 day
(first date of next month) for several sheets. My main goal is to get the
minimum value for every month.
First, how to extract those data to list of data frame based on year and
delete every overlapping date?
Second, how to sort it based on date with ascending order (old to new)?
Third, how to get the maximum together with the date?
I did this one,
...
file.list <- list.files(pattern='*.xlsx')
file.list<-mixedsort(file.list)
#
https://stackoverflow.com/questions/12945687/read-all-worksheets-in-an-excel-workbook-into-an-r-list-with-data-frames
read_excel_allsheets <- function(filename, tibble = FALSE) {
sheets <- readxl::excel_sheets(filename)
x <- lapply(sheets, function(X) read.xlsx(filename, sheet=X, rows=9:40,
cols=1:2))
if(!tibble) x <- lapply(x, as.data.frame)
names(x) <- sheets
x
}
pon<-lapply(file.list, function(i) read_excel_allsheets(i))
pon1<-do.call("rbind",pon)
names(pon1) <- paste0("M.", 1:360)
pon1 <-lapply(pon1,function(x){x$RR[x$RR==] <- NA; x})
pon1 <-lapply(pon1,function(x){x$RR[x$RR==""] <- NA; x})
maxi<-lapply(pon1, function(x) max(x$RR,na.rm=T))
maxi<-data.frame(Reduce(rbind, maxi))
names(maxi)<-"maxi"
but the list start from January for every year, and move to February and so
on. And there is no date in "maxi". Here some sample what I get from my
simple code.
pon1[256:258]$M.256
Tanggal RR
1 01-09-2001 5.2
2 02-09-2001 0.3
3 03-09-2001 29.0
4 04-09-2001 0.7
5 05-09-2001 9.6
6 06-09-2001 0.7
7 07-09-2001 NA
8 08-09-2001 13.2
9 09-09-2001 NA
10 10-09-2001 NA
11 11-09-2001 0.0
12 12-09-2001 66.0
13 13-09-2001 0.0
14 14-09-2001 57.6
15 15-09-2001 18.0
16 16-09-2001 29.2
17 17-09-2001 52.2
18 18-09-2001 7.0
19 19-09-2001 NA
20 20-09-2001 74.5
21 21-09-2001 20.3
22 22-09-2001 49.6
23 23-09-2001 0.0
24 24-09-2001 1.3
25 25-09-2001 0.0
26 26-09-2001 1.0
27 27-09-2001 0.1
28 28-09-2001 1.9
29 29-09-2001 9.5
30 30-09-2001 3.3
31 01-10-2001 0.0
$M.257
Tanggal RR
1 01-09-2002 0.0
2 02-09-2002 0.0
3 03-09-2002 0.0
4 04-09-2002 12.8
5 05-09-2002 1.0
6 06-09-2002 0.0
7 07-09-2002 NA
8 08-09-2002 22.2
9 09-09-2002 NA
10 10-09-2002 NA
11 11-09-2002 0.0
12 12-09-2002 0.0
13 13-09-2002 0.0
14 14-09-2002 NA
15 15-09-2002 0.0
16 16-09-2002 0.0
17 17-09-2002 0.0
18 18-09-2002 13.3
19 19-09-2002 0.0
20 20-09-2002 0.0
21 21-09-2002 0.0
22 22-09-2002 0.0
23 23-09-2002 0.0
24 24-09-2002 0.0
25 25-09-2002 0.0
26 26-09-2002 0.5
27 27-09-2002 2.1
28 28-09-2002 NA
29 29-09-2002 18.5
30 30-09-2002 0.0
31 01-10-2002 NA
$M.258
Tanggal RR
1 01-09-2003 0.0
2 02-09-2003 0.0
3 03-09-2003 0.0
4 04-09-2003 4.0
5 05-09-2003 0.3
6 06-09-2003 0.0
7 07-09-2003 NA
8 08-09-2003 0.0
9 09-09-2003 0.0
10 10-09-2003 0.0
11 11-09-2003 NA
12 12-09-2003 1.0
13 13-09-2003 0.0
14 14-09-2003 60.0
15 15-09-2003 4.5
16 16-09-2003 0.1
17 17-09-2003 2.1
18 18-09-2003 NA
19 19-09-2003 0.0
20 20-09-2003 NA
21 21-09-2003 NA
22 22-09-2003 31.5
23 23-09-2003 42.0
24 24-09-2003 43.3
25 25-09-2003 2.8
26 26-09-2003 21.4
27 27-09-2003 0.8
28 28-09-2003 42.3
29 29-09-2003 5.3
30 30-09-2003 17.3
31 01-10-2003 0.0
Any lead or help is very appreciate.
Best,
Ani
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
