Re: [R-br] Residuos do Teste T de Student

[Graciliano via R-br]

Olá Cesar, 

Então, não falam nada de resíduo. Pensei que o resíduo da analise seria algo 
relacionado a dispersão das observações em relação à média de cada tratamento. 

Na verdade minha preocupação é com a premissa de normalidade. Como geralmente é 
testado o resíduo da análise, fiquei sem rumo no caso do teste T.
  No stackexchange.com falam sobre testar normalidade de Y, mas isso não faz 
muito sentido se devemos ver o comportamento de Y sobre efeito de X. 

Espero estar sendo claro.

-Mensagem Original-
De: "Cesar Rabak" 
Enviada em: ‎19/‎11/‎2016 15:20
Para: "Graciliano Galdino" ; "a lista Brasileira oficial de 
discussão do programa R." 
Assunto: Re: [R-br] Residuos do Teste T de Student

Gracialiano,


Por que você esperaria resíduos de um teste t de Student?


Você conhece alguma referência ou obra que trate desse teste e fale sobre 
resíduos calculados devido a esse teste?




2016-11-19 15:55 GMT-02:00 Graciliano Galdino via R-br 
:

Olá a todos,


Quando faço uma análise qualquer usando teste T (t.test) não aparecem os 
resíduos na lista dos resultados da análise, como aparece na ANOVA. Como faço 
para extraí-los? 



Abraço
-- 

Graciliano Galdino A. dos Santos
Biólogo
Doutorando em Ciências Florestais - PPGCF
Universidade Federal Rural da Amazônia - UFRA

___
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.___
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.

Re: [R-br] Ajuste de modelo não linear

[Graciliano via R-br]

Olá Dr. Walmes, obrigado por compartilhar esse script do rpanel, muito legal! 
Mas ó, não consegui achar a função  as.lm no caminho indicado. Daria para 
anexa-la aqui?

Att.

-Mensagem Original-
De: "Walmes Zeviani via R-br" 
Enviada em: ‎17/‎10/‎2016 08:33
Para: "Andre Oliveira" ; "a lista Brasileira oficial 
de discussão do programa R." 
Assunto: Re: [R-br] Ajuste de modelo não linear

André,

Apenas para complementar a solução do Marcus, encaminho o código que faz uso de 
uma função com recursos do pacote rpanel para ajuste de modelo não linear de 
forma interativa. Essa função (na versão antiga) gerou o GIF que fica na capa 
da página da lista no Nabble (rs). O uso da interface não elimina o problema do 
chute inicial, só diminui porque ao invés de passar um valor você passa um 
intervalo e calibra visualmente. Além disso, tem a parte de análise de resíduos.

Brasil <- c(13.63, 13.93, 14.09, 15.12, 15.98, 16.08, 16.45, 17.96,
18.83, 19.28, 19.18, 20.89, 21.87, 21.85, 23.43, 25.38,
25.64, 26.99, 28.68, 31.92, 32.92, 33.32, 34.37, 34.99,
36.08, 37.72, 40.41, 42.81, 44.18, 44.91, 48.97, 51.59,
53.44, 56.49, 62.2, 66.49, 67.01, 68.41, 74.14, 77.58,
80.39, 85.71, 89.58, 96.3, 100.9, 102.8, 107.69, 111.51,
111.42, 116.59, 118.85, 121.42, 127.53)
anos <- 1961:2013

anos2 <- anos - min(anos) + 1

# Calibrando o chute inicial (por tentativa-erro faça aproximar).
start <- list(a = 12, b = 1.049)
plot(Brasil ~ anos2)
with(start,
 curve(a * (b^(anos2)),
   xname = "anos2",
   add = TRUE,
   col = 2))

fit <- nls(Brasil ~ a * (b^(anos2)),
   start = start)

# Verificando o ajuste.
plot(Brasil ~ anos2)
with(as.list(coef(fit)),
 curve(a * (b^(anos2)),
   xname = "anos2",
   add = TRUE,
   col = 4))

# rp.nls: função que usa recursos no pacote rpanel para ajuste de
# modelos de regressão não linear.
source(paste0("https://raw.githubusercontent.com/walmes/;,
  "wzRfun/master/R/rp.nls.R"))
library(rpanel)

fit <- rp.nls(Brasil ~ a * (b^(anos2)),
  data = data.frame(Brasil = Brasil, anos2 = anos2),
  start = list(a = c(5, 20),
   b = c(0.95, 1.5)))
summary(fit)

# Função que era do pacote nls2() mas foi removida (guardei o fonte).
source("~/Dropbox/CursoR/GeneticaEsalq/as.lm.R")
library(nls2)

mit <- as.lm(fit)

# Diagnóstico do modelo.
par(mfrow = c(2, 2))
plot(mit)
layout(1)

# ! Relação média variância !

# Linear!
plot(log(Brasil) ~ anos2)

À disposição.

Walmes.



​___
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.

Re: [R-br] boxplot comparativos

[Graciliano via R-br]

Estou no celular e não consigo o script, mas no ggplot2 dá facinho. Bota o mês 
no eixo x, temperatura no y e o comando "fill" por cidade. 

-Mensagem Original-
De: "Elisa Henning via R-br" 
Enviada em: ‎14/‎10/‎2016 22:06
Para: "a lista Brasileira oficial de discussão do programa R." 

Assunto: [R-br] boxplot comparativos

Boa noite


Gostaria de construir boxplots múltiplos comparativos.
Exemplificando: tenho dois arquivos com dados climáticos referentes a duas 
cidades.
Gostaria de comparar, mês a mês, por meio de boxplots, as temperaturas destas 
duas cidades.


Alguém  poderia  indicar algum material que auxilie a fazer estes gráficos?


Desde já agradeço


Elisa___
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.