<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Walmes,<br>
<br>
Fiz a predição sobre um modelo linear abaixo, isto parece
correto:<br>
<br>
## Dados artificiais<br>
z1 <- rnorm(n=100, mean=0.73436, sd=0.104739)<br>
z2 <- rnorm(n=100, mean=0.60173, sd=0.054492)<br>
z <- c(z1,z2)<br>
x1<-rep("sadio",length(z)/2)<br>
x2<-rep("atacado",length(z)/2)<br>
x<-c(x1,x2)<br>
<br>
status <- rep("sadio", length(z))<br>
i <- 0.69173+c(-1,1)*0.008838306<br>
fi <- findInterval(z, i)<br>
status[fi==1] <- "atacado"<br>
<br>
<br>
#Ajuste do modelo<br>
m.1<-lm(z~x)<br>
prd<-predict(m.1, interval=c("prediction"))<br>
<br>
## Representação gráfica<br>
plot(z, col=fi+1)<br>
abline(h=i, lty=2)<br>
matlines(prd[,c("lwr","upr")], col="red")<br>
<br>
<br>
<br>
Obrigado,<br>
<br>
Alexandre<br>
<div class="moz-cite-prefix">Em 02/07/2014 12:20, ASANTOS escreveu:<br>
</div>
<blockquote cite="mid:53B43138.3070007@yahoo.com.br" type="cite">
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
Exatamente Walmes,<br>
<br>
Eu estava fazendo errado como comentou, queria
classificar novas observações para áreas não utilizadas para
calcular a média e o intervalo de confiança. Uma particularidade
dos dados, são que são oriundos de um cálculo de índice de índice
de vegetação por diferença normalizada (NDVI), portanto, qualquer
valor abaixo de 0.69173 indicam áreas atacadas e acima disto
sadias. Teria alguma sugestão de que abordagem eu poderia utilizar
para fazer a predição de áreas atacadas ou não?<br>
<br>
Obrigado,<br>
<br>
Alexandre<br>
<br>
<div class="moz-cite-prefix">Em 02/07/2014 10:24, walmes .
escreveu:<br>
</div>
<blockquote
cite="mid:CAFU=EkZjhF2Sr0OapePTc-CLyH=FzP+1TFd+CNioPAmkjQAkdA@mail.gmail.com"
type="cite">
<div dir="ltr">
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif">Pensando que você pode aumentar o número de
classes, a findInterval() é interessante.<br>
<br>
<span style="font-family:courier new,monospace">z1 <-
rnorm(n=100, mean=0.73436, sd=0.104739)<br>
z2 <- rnorm(n=100, mean=0.69173, sd=0.104492)<br>
z <- c(z1,z2)<br>
<br>
status <- rep("sadio", length(z))<br>
i <- 0.69173+c(-1,1)*0.008838306<br>
fi <- findInterval(z, i)<br>
status[fi==1] <- "atacado"<br>
<br>
plot(z, col=fi+1)<br>
abline(h=i, lty=2)<br>
</span><br>
</div>
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif">Vai uma dúvida minha que tá fora do escopo.
Você tá classificando novas observações (observações
futuras) pelo intervalo de confiança (para um parâmetro)
gerado em uma análise. Talvez isso não seja adequado. O
nível de confiança desse intervalo (95%) não significa que
95% das observações serão classificadas como "atacada". 95%
é o nível de cobertura do intervalo que representa a
probabilidade dele conter o parâmetro. Não seria mais
apropriado um intervalo de predição (para observação
futura)?<br>
<br>
À disposição.<br>
</div>
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif">Walmes.<br>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a moz-do-not-send="true" class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a moz-do-not-send="true" class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a moz-do-not-send="true" class="moz-txt-link-freetext" href="http://lattes.cnpq.br/1360403201088680">http://lattes.cnpq.br/1360403201088680</a>
======================================================================
</pre>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<a class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a class="moz-txt-link-freetext" href="http://lattes.cnpq.br/1360403201088680">http://lattes.cnpq.br/1360403201088680</a>
======================================================================
</pre>
</body>
</html>