[R-es] calculo de otros percentiles (0.95, 0.99)

[fcortes en inidep.edu.ar]

Hola Alex, no revisé tus datos pero quizás esto ayude,

>mis.quantiles<-function (x) {quantile(x,probs = c(0.95,0.99))}

>aggregate(datos$TSM, list(DIA,MES),mis.quantiles,na.rm=T))


Cambiando el argumento probs vas a obtener los quantiles que quieras.

Saludos

Federico





Dr. Federico Cortés
fcortes en inidep.edu.ar
Programa Pesquerías de Condrictios
Instituto Nacional de Investigación y Desarrollo Pesquero (INIDEP)
Tel: +54(223)4862586 int.282 Fax: +54(223)4861830
Paseo Victoria Ocampo Nº 1, (7600) Mar del Plata,  Argentina.


Estimados junto con saludarlos,
Tengo una serie tiempo de datos de temperatura del mar horaria desde el 
2004. Y necesito obtener el percentil 0.95 y 0.99 diario. Sin embargo, 
solo he logrado sacar los percentiles tipicos 0.25,0.5,0.75 usando la 
funcion aggregate:

datos<-read.table(file.choose(), header=T)
attach(datos)
aggregate(datos$TSM, list(DIA,MES),quantile,na.rm=T))


he intentando con la funcion quantile pedir estos percentiles (0.99 y 
0.95) pero no me los entrega :( .... Alguno podria ayudarme a salir del 
entuerto..... subo el txt de los datos.... muchas muchas gracias de 
Chile



--------------------------------------------------------------------------------


_______________________________________________
R-help-es mailing list
R-help-es en r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es