fcortes en inidep.edu.ar
2014-Sep-23 22:08 UTC
[R-es] calculo de otros percentiles (0.95, 0.99)
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