Eu encontrei a forma de calcular os quantis que precisava mas não conheço nem percebo os erros que me deram
e os resultados de "IRQ" não me aparecem.
código:
# Estatística por ano
porano<-apply(dados[,-1],2,function(x)c(summary(x),0)[1:7])
desvio padrão<-sd(dados[,-1], na.rm= TRUE)
Error: unexpected symbol in "desvio padrão"
variância<-diag(var(dados[,-1], na.rm= TRUE))
IQR) <- quantile(dados[,-1], 3/4) -
Error: unexpected ')' in "IQR)"
quantile(dados[,-1], 1/4)
Error in quantile.default(dados[, -1], 1/4) :
missing values and NaN's not allowed if 'na.rm' is FALSE
amplitude<-apply(dados[,-1],2,range, na.rm= TRUE)
amplitudediff<-c(amplitude[2,] - amplitude[1,])
rownames(porano) <- c('Minimo','Primeiro quartil','Mediana','Média','Terceiro Quartil','Máximo','NAs')
poranol<-rbind(porano,desvio padrão,variância,amplitude,amplitudediff,IQR)
Error: unexpected symbol in "poranol<-rbind(porano,desvio padrão"
colnames(porano) = 2008:2012X2008 X2009 X2010 X2011 X2012
poranol
Minimo 66.000000 65.000000 62.000000 59.000000 60.000000
Primeiro quartil 76.000000 84.000000 86.000000 88.250000 91.000000
Mediana 80.000000 87.000000 89.500000 91.500000 94.000000
Média 80.270000 86.290000 88.500000 90.500000 91.840000
Terceiro Quartil 85.000000 89.500000 92.000000 95.000000 96.000000
Máximo 94.000000 98.000000 98.000000 98.000000 99.000000
NAs 4.000000 3.000000 3.000000 3.000000 0.000000
desviopadrão 7.362374 6.254838 6.542402 7.439656 7.879345
variância 54.204545 25.746212 21.530303 25.130682 12.780303
66.000000 65.000000 62.000000 59.000000 60.000000
94.000000 98.000000 98.000000 98.000000 99.000000
amplitudediff 28.000000 33.000000 36.000000 39.000000 39.000000
porano2008 2009 2010 2011 2012
Minimo 66.00 65.00 62.0 59.00 60.00
Primeiro quartil 76.00 84.00 86.0 88.25 91.00
Mediana 80.00 87.00 89.5 91.50 94.00
Média 80.27 86.29 88.5 90.50 91.84
Terceiro Quartil 85.00 89.50 92.0 95.00 96.00
Máximo 94.00 98.00 98.0 98.00 99.00
NAs 4.00 3.00 3.0 3.00 0.00
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.