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 minus the right side. It?s a numeric variable. But i?ve got no data for some individuals, so in the variable "fem" i got numbers and smoe "NA". When i try to apply kurtosis to "fem", the NA, don?t allow me to do that, that`s what i get:> kurtosis(fem)[1] NA I tried with simple tricks, but i?ve got no results:> kurtosis(fem,na.rm=T)Error en kurtosis(fem, na.rm = T) : el argumento(s) no fue utilizado(s) (na.rm = T) What can i do? the values of fem vary from -4 to 9. Thanks very much> kurtosis(fem[fem!="NA"])[1] NA ----- 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: http://r.789695.n4.nabble.com/Applying-a-function-to-a-subset-tp3558567p3558567.html Sent from the R help mailing list archive at Nabble.com.
kurtosis(fem[!is.na(fem)]) Sent from my iPad On May 29, 2011, at 4:54, gaiarrido <gaiarrido at 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 } > > Then create the object fem, which is the difference between the count of a > trait in the left side of the body minus the right side. It?s a numeric > variable. But i?ve got no data for some individuals, so in the variable > "fem" i got numbers and smoe "NA". > When i try to apply kurtosis to "fem", the NA, don?t allow me to do that, > that`s what i get: > >> kurtosis(fem) > [1] NA > > I tried with simple tricks, but i?ve got no results: >> kurtosis(fem,na.rm=T) > Error en kurtosis(fem, na.rm = T) : > el argumento(s) no fue utilizado(s) (na.rm = T) > > What can i do? the values of fem vary from -4 to 9. > > Thanks very much >> kurtosis(fem[fem!="NA"]) > [1] NA > > > ----- > 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: http://r.789695.n4.nabble.com/Applying-a-function-to-a-subset-tp3558567p3558567.html > Sent from the R help mailing list archive at Nabble.com. > > ______________________________________________ > R-help at r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code.
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, 2011, at 4:54, gaiarrido<gaiarrido at 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 } >> >> Then create the object fem, which is the difference between the count of a >> trait in the left side of the body minus the right side. It?s a numeric >> variable. But i?ve got no data for some individuals, so in the variable >> "fem" i got numbers and smoe "NA". >> When i try to apply kurtosis to "fem", the NA, don?t allow me to do that, >> that`s what i get: >> >>> kurtosis(fem) >> [1] NA >> >> I tried with simple tricks, but i?ve got no results: >>> kurtosis(fem,na.rm=T) >> Error en kurtosis(fem, na.rm = T) : >> el argumento(s) no fue utilizado(s) (na.rm = T) >> >> What can i do? the values of fem vary from -4 to 9. >> >> Thanks very much >>> kurtosis(fem[fem!="NA"]) >> [1] NA >> >> >> ----- >> 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: http://r.789695.n4.nabble.com/Applying-a-function-to-a-subset-tp3558567p3558567.html >> Sent from the R help mailing list archive at Nabble.com. >> >> ______________________________________________ >> R-help at r-project.org mailing list >> https://stat.ethz.ch/mailman/listinfo/r-help >> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html >> and provide commented, minimal, self-contained, reproducible code. > > ______________________________________________ > R-help at r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code.