Displaying 1 result from an estimated 1 matches for "heavyprec".
2010 Mar 26
0
CDF calculation from kernel density estimates for a 324X 15 matrix
Hi,
I have a 324X15 matrix (No of years vs. heavy precipitation days) and I want
to calculate the cdf at 5 different data points for each row. I tried by the
following codes but it's not working.
heavyprec <- read.csv (file="heavyprecdays_CSV.csv",header=TRUE,sep=",")
a <- heavyprec
pdf <- density (a, bw="SJ", kern= "gaussian")
f <- approxfun (pdf$x, pdf$y)
cdf <- integrate (f, -Inf, 25.25) # I also have to calculate cdf at 30,
33, 37 heavy p...