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

"Marcuzzi, Javier Rubén" javier.ruben.marcuzzi en gmail.com
Mie Sep 24 00:02:45 CEST 2014


Estimado Alex Fabian

Vecino, ya que saludas desde Chile, yo estoy en Argentina..

Copio y pego un código, no llega a la solución porque me da un mensaje 
de error, pero se aproxima, en este momento no puedo dedicarle más 
tiempo como para solucionarlo pero seguramete usted o algún otro miembro 
de la lista terminará la corrección.

#datos<-read.table(file.choose(), header=T)
datos <- read.table("C:/Users/Javier/Downloads/TSM_AN_2004.txt", 
header=TRUE, quote="")
attach(datos)

aggregate(datos$TSM, by=list(DIA,MES),quantile(datos,c(.95, 
.99)),na.rm=FALSE)

duration = faithful$eruptions
quantile(duration, c(.32, .57, .98))

quantile(datos$MES, c(.32, .57, .98))

Saludos.
Javier Marcuzzi

El 23/09/2014 06:37 p.m., ALEX FABIAN Mellado escribió:
> 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


	[[alternative HTML version deleted]]



Más información sobre la lista de distribución R-help-es