Here´s my problem,
i have developed the function kurtosis (using R-book as a guide) with this
commands:
kurtosis-function(x) {
m4-sum((x-mean(x))^4)/length(x)
s4-var(x)^2
m4/s4 - 3 }
Then create the object fem, which is the difference between the count of a
trait in the left side of the body
kurtosis(fem[!is.na(fem)])
Sent from my iPad
On May 29, 2011, at 4:54, gaiarrido gaiarr...@usal.es wrote:
Here´s my problem,
i have developed the function kurtosis (using R-book as a guide) with this
commands:
kurtosis-function(x) {
m4-sum((x-mean(x))^4)/length(x)
s4-var(x)^2
m4/s4 - 3 }
Or just include is.na=TRUE in the definition of kurtosis():
kurtosis-function(x) {
m4-sum((x-mean(x,na.rm=TRUE))^4,na.rm=TRUE)/length(x)
s4-var(x,na.rm=TRUE)^2
m4/s4 - 3 }
HTH
Stephan
Am 29.05.2011 11:34, schrieb Jim Holtman:
kurtosis(fem[!is.na(fem)])
Sent from my iPad
On May 29,
That's it.
thanks very much, I save hours if work qith your reply
-
Mario Garrido Escudero
PhD student
Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola
Universidad de Salamanca
--
View this message in context:
4 matches
Mail list logo