Elias,
Você precisará ajustar o script de acordo com sua conveniência. Sei que
kurtosis() e skewness() do pacote {moments} aceitam o argumento na.rm,
tendo FALSE por padrão.
================================================
Éder Comunello
Agronomist (UEM), MSc in Environ. Sciences (UEM)
DSc in Agricultural Systems Engineering (USP/Esalq)
Brazilian Agricultural Research Corporation (Embrapa)
Dourados, MS, Brazil |<O>|
================================================
GEO, -22.2752, -54.8182, 408m
UTC-04:00 / DST: UTC-03:00
Em 11 de maio de 2016 09:35, Elias Carvalho <[email protected]> escreveu:
> Fantástico codigo Fernando vou usar com certeza.
>
> Me tire uma dúvida: meu data set tem bastante Missing, eu devo removê-los
> antes de verificar a normalidade ou a função em si já trata isso ?
>
> Em 6 de maio de 2016 09:12, Elias Carvalho <[email protected]> escreveu:
>
>> Bom dia Pessoal
>>
>> Eu tenho o seguinte data frame:
>>
>> 'data.frame': 1999 obs. of 14 variables:
>> $ AGE: int 21 31 44 46 49 50 52 64 23 33 ...
>> $ PHA: Factor w/ 2 levels "NO","YES": 1 1 1 1 1 1 1 1 2 2 ...
>> $ SMK: Factor w/ 2 levels "NO","YES": 1 1 1 1 1 1 1 1 1 1 ...
>> $ ALC: Factor w/ 2 levels "NO","YES": 1 1 1 1 1 NA 2 NA 1 NA ...
>>
>>
>>
>>
>> E preciso rodar uma rotina para verificar a normalidade de cada variável
>> de uma só vez, então fiz assim (ps: me perdoem o 'for' eu sou programador
>> das antigas e ainda estou aprendendo a user os ...apply):
>>
>> for (x in 1:length(names(data)))* # vai de 1 a 4 passando por todas as
>> variáveis *
>> {
>> * # Insere o nome da variável em uma matriz*
>> commandP <- paste("nmatrix[",x,",1] <- names(data)[x]", sep="")
>> eval(parse(text=commandP))
>> * # aplicar apenas para variaveis continuas*
>> * if (!is.factor(data$names(data)[x])) *
>> {
>> *# na coluna 2 da matriz insere valores referente a skewness*
>> commandP <- paste("nmatrix[",x,",2] <- round(skewness(data$",
>> names(data)[x],"),2)", sep="")
>> eval(parse(text=commandP))
>> *# na coluna3 da matriz insere valores referente a kurtosis*
>> commandP <- paste("nmatrix[",x,",3] <- round(kurtosis(data$",
>> names(data)[x],"),2)", sep="")
>> eval(parse(text=commandP))
>> * # gera um histograma para cada variável* já com par(mfrow...)
>> calculado de acordo com o
>> # numero de variáveis
>> commandP <- paste("hist(data$", names(data)[x],")", sep="")
>> eval(parse(text=commandP))
>> *# gera um gráfico de quartil para cada variável já com
>> par(mfrow...) calculado de *
>> * # acordo **com a quantidade de variáveis*
>> commandP <- paste("qqnorm(data$", names(data)[x],")", sep="")
>> eval(parse(text=commandP))
>> * # Insere a linha no gráfico *
>> commandP <- paste("abline(0,1)", sep="")
>> eval(parse(text=commandP))
>> }
>>
>> O Resultado que quero parecido com o abaixo, mas para todas as variáveis:
>>
>> nmatrix [,1] [,2] [,3]
>> [1,] "AGE" "-0.13" "2.28"
>>
>>
>>
>>
>>
>>
>> O meu problema é na linha if (!is.factor(data$names(data)[x])) que
>> identifica todas as variáveis como caracter, o que não está errado, pois
>> pega por exemplo "AGE" e então sempre retorna falso e na verdade eu queria
>> algo como "!is.factor(data$AGE), !is.factor(data$SMK)....
>>
>> ou se alguem conhecer um pacote que faça isso de maneira mais fácil
>> agradeço a sugestão
>>
>> Obrigado
>> -
>>
>> Elias
>> <http://lattes.cnpq.br/4248328961021251>
>>
>
>
>
> --
> Best regards... 8^)
>
> “The mind that is open to new ideas never come back
> to its original size” *Albert Einstein*
>
>
> --
> Obrigado
> Elias
>
> _______________________________________________
> R-br mailing list
> [email protected]
> 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.
>
_______________________________________________
R-br mailing list
[email protected]
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.
