Auxilio com Funções de Estatistica
Bom dia! Meu nome é Yury Duarte e sou pós graduando em Engenharia de Biossistemas. Estou com uma necessidade de aplicar alguns índices estatísticos a uma serie de dados, entretanto, não consegui encontrar uma função no R que as faça. Os índices são: - R² - r - c - D - RMSE - Erro médio - Erro médio absoluto Dessa forma, gostaria de saber se já existem comandos para calcular esses índices ou se terei que desenvolver o código para cada um deles. Desde já, agradeço pela colaboração! Att Yury Duarte Engenheiro Agrônomo - ESALQ/USP
Se você fizer um modelo com a função lm por exemplo, da para extrair a maioria desses indices com summary ou anova, apesar que confesso que não sei o que são algumas dessas letras que você colocou. Segue um exemplo. set.seed(123) x <- runif(30, 2, 5) y <- rnorm(30, x, 0.3) modelo <- lm(y ~ x) summary(modelo) anova(modelo) ###RSME sqrt(sum(resid(modelo)^2)/30) ###RSME sqrt(mean((y-fitted(modelo))^2)) Espero ter ajudado. Em 9 de novembro de 2015 08:53, Yury Duarte <yurynepomuceno@gmail.com> escreveu:
Bom dia! Meu nome é Yury Duarte e sou pós graduando em Engenharia de Biossistemas. Estou com uma necessidade de aplicar alguns índices estatísticos a uma serie de dados, entretanto, não consegui encontrar uma função no R que as faça. Os índices são:
- R² - r - c - D - RMSE - Erro médio - Erro médio absoluto
Dessa forma, gostaria de saber se já existem comandos para calcular esses índices ou se terei que desenvolver o código para cada um deles. Desde já, agradeço pela colaboração!
Att
Yury Duarte Engenheiro Agrônomo - ESALQ/USP
_______________________________________________ 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.
-- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056
Obrigado, Augusto! Irei tentar esses comandos que vc sugeriu. Os outros índices que vc comentou, que não coloquei o nome, são: D: índice de concordância de Willmott c: índice de confiança ou desempenho Mais uma vez, agradeço pela ajuda! Att Yury Duarte Engenheiro Agrônomo - ESALQ/USP Em 9 de novembro de 2015 10:46, Augusto Ribas <ribas.aca@gmail.com> escreveu:
Se você fizer um modelo com a função lm por exemplo, da para extrair a maioria desses indices com summary ou anova, apesar que confesso que não sei o que são algumas dessas letras que você colocou. Segue um exemplo.
set.seed(123) x <- runif(30, 2, 5) y <- rnorm(30, x, 0.3) modelo <- lm(y ~ x)
summary(modelo) anova(modelo)
###RSME sqrt(sum(resid(modelo)^2)/30)
###RSME sqrt(mean((y-fitted(modelo))^2))
Espero ter ajudado.
Em 9 de novembro de 2015 08:53, Yury Duarte <yurynepomuceno@gmail.com> escreveu:
Bom dia! Meu nome é Yury Duarte e sou pós graduando em Engenharia de Biossistemas. Estou com uma necessidade de aplicar alguns índices estatísticos a uma serie de dados, entretanto, não consegui encontrar uma função no R que as faça. Os índices são:
- R² - r - c - D - RMSE - Erro médio - Erro médio absoluto
Dessa forma, gostaria de saber se já existem comandos para calcular esses índices ou se terei que desenvolver o código para cada um deles. Desde já, agradeço pela colaboração!
Att
Yury Duarte Engenheiro Agrônomo - ESALQ/USP
_______________________________________________ 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.
-- Grato Augusto C. A. Ribas
Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056
_______________________________________________ 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.
Use ## sum(resid(modelo)^2): use deviance(). deviance(modelo) ## estimativa de sigma, use: summary(modelo)$sigma ## R² summary(modelo)$r.squared summary(modelo)$adj.r.squared Os demais você calcula a partir dos resíduos que estraí com residuals(modelo). Pode ser que o pacote car (companium of applied regression) tenha funções diretas para essas medidas. À disposição. Walmes.
participantes (3)
-
Augusto Ribas -
Walmes Zeviani -
Yury Duarte