Caentes?<br>Prezado Eder,<br>Se for regressão linear simples como você descreveu, o delineamento mais informativo (como o Walmes chamou a atenção) é o com níveis extremos. Um nível intermediário é recomendável para testar a falta de ajuste do modelo de regressão.<br>
No caso, suspeito que seu problema pode ser mais sutil do que isto.<br>Qual o interesse em tomar medidas em tempos diferentes? é por causa do custo de armazenamento de dados? ou por causa dos reagentes?<br>Qual é o recurso mais escasso?<br>
[]s<br>Júlio<br><br><br><div class="gmail_quote">2011/9/15 <span dir="ltr"><<a href="mailto:r-br-request@listas.c3sl.ufpr.br">r-br-request@listas.c3sl.ufpr.br</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Enviar submissões para a lista de discussão R-br para<br>
<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
<br>
Para se cadastrar ou descadastrar via WWW, visite o endereço<br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
ou, via email, envie uma mensagem com a palavra 'help' no assunto ou<br>
corpo da mensagem para<br>
<a href="mailto:r-br-request@listas.c3sl.ufpr.br">r-br-request@listas.c3sl.ufpr.br</a><br>
<br>
Você poderá entrar em contato com a pessoa que gerencia a lista pelo<br>
endereço<br>
<a href="mailto:r-br-owner@listas.c3sl.ufpr.br">r-br-owner@listas.c3sl.ufpr.br</a><br>
<br>
Quando responder, por favor edite sua linha Assunto assim ela será<br>
mais específica que "Re: Contents of R-br digest..."<br>
<br>
<br>
Tópicos de Hoje:<br>
<br>
1. União de poligonos (Jose Vilton Costa)<br>
2. Seleção Níveis Regressão linear (Eder David Borges da Silva)<br>
3. RES: Seleção Níveis Regressão linear (Alan Rodrigo Panosso)<br>
4. Re: RES: Seleção Níveis Regressão linear (Walmes Zeviani)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Wed, 14 Sep 2011 12:17:40 -0300<br>
From: Jose Vilton Costa <<a href="mailto:josevilton@hotmail.com">josevilton@hotmail.com</a>><br>
To: <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: [R-br] União de poligonos<br>
Message-ID: <SNT113-W5836662C704F025AAC51DFDA040@phx.gbl><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
<br>
Professor Paulo,Eu já estava tentando usar a função unionSpatialPolygons() , mas não estava conseguindo informar os IDs. Agora usando a biblioteca aRT foi possível fazer isto. Mauro, tente fazer um unpdate no pacote Maptools para ver se não funciona a função readShapePolys(). Obrigado,José Vilton. Estava vendo vossa discussão, tentei rodar o script e não encontra esta função readShapePolys(), e já tenho estas bibliotécas instaladas!<br>
<br>
> require(maptools)<br>
> readShapePolys()<br>
Erro: não foi possível encontrar a função "readShapePolys"<br>
<br>
<br>
Uma possibilidade é o uso das funcionalidades<br>
do pacote maptools que voce já está usando<br>
<br>
As funcoes relevantes seriam<br>
<br>
require(maptools)<br>
readShapePolys()<br>
gpclibPermit() # pode ser evitado se voce tiver a biblioteca rgeos instalada<br>
unionSpatialPolygons()<br>
<br>
O aRT (<a href="http://www.leg.ufpr.br/aRT" target="_blank">www.leg.ufpr.br/aRT</a>) tb tem funcoes para isto<br>
<br>
A ideia básica é que na funcao de uniao voce passa os ID's que quer unir,<br>
ele vai gerar um novo cunjunto de IDs no objeto resultante<br>
<br>
<br>
<br>
<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110914/37d99bd9/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110914/37d99bd9/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Wed, 14 Sep 2011 22:43:51 -0300<br>
From: Eder David Borges da Silva <<a href="mailto:eder@leg.ufpr.br">eder@leg.ufpr.br</a>><br>
To: <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
Subject: [R-br] Seleção Níveis Regressão linear<br>
Message-ID:<br>
<<a href="mailto:CALKkXXorwhyjN_-wu-Fe_QdHUwK6EfXNKoe6zHfwvLNTREorXg@mail.gmail.com">CALKkXXorwhyjN_-wu-Fe_QdHUwK6EfXNKoe6zHfwvLNTREorXg@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Boa noite pessoal,<br>
Gostaria de discutir com vocês um experimento.<br>
Estamos planejando um experimento de medição de gases do solo via câmaras,<br>
em síntese colocar uma câmara sob o solo e a cada intervalo de tempo medi-se<br>
a concentração de gases (CO2, N2O, CH4) retirando uma amostra e fazendo<br>
cromatografia dela.<br>
nosso experimento visa otimizar os tempos para coleta destas amostra hoje se<br>
faz 0, 15, 30, 45 min.<br>
Vamos fazer um experimento com 13 tempos de 0 a 120 min.<br>
Nosso interesse neste experimento é determinar a inclinação da reta que da a<br>
taxa de emissão (ppm / min)<br>
Minha sugestão de analise:<br>
1) Ajustar um modelo com os 13 tempos. sendo esta estimativa de inclinação a<br>
mais próxima da realidade. (M0)<br>
2) ajustar modelos com tempos distintos por exemplo 0, 5, 10 , 20 (Mi) e<br>
comparar esta estimativas com o (M0)<br>
repetir o passo 2 varias vezes<br>
3) computar os modelos que são "iguais" a M0, de todos os "iguais" o melhor<br>
é aquele que demore menos tempo.<br>
Em resumo se avaliar nos tempos 0,5,10,20 e igual a avaliar em 13 tempos de<br>
0 a 120???<br>
OBS: o tempo 0 sempre terá de ser medido<br>
Um Codigo possivel seria isso:<br>
require(MCMCpack)<br>
require(plyr)<br>
require(latticeExtra)<br>
<br>
dados <-<br>
expand.grid(local=c('Agri','Pastagem'),epoca=factor(c(-1,1,5,10,20)),camara=c('C1','C2','C3'),tempo=c(seq(0,45,by=5),60,90,120),bloco=factor(1:4))<br>
dim(dados)<br>
dados$resp <- rnorm(nrow(dados),10,10)<br>
head(dados)<br>
<br>
m0 <- MCMCregress(resp~tempo,data=dados)<br>
summary(m0)<br>
plot(m0)<br>
###Distribuiçao pretendida<br>
<br>
#### Simulação<br>
simu <- list(A=c(0,5,10,60),<br>
B=c(0,5,15,45),<br>
C=c(0,10,20,45),<br>
D=c(0,10,45,60))### colocar mais<br>
<br>
resS <- data.frame(V1='m0',V2=m0[,2])<br>
<br>
for(i in 1:length(simu)){<br>
m0 <- MCMCregress(resp~tempo,data=subset(dados,dados$tempo %in%<br>
simu[[i]]),burnin=1000,mcmc=5000)<br>
m0 <- as.matrix(m0)<br>
resS <- rbind(resS,data.frame(V1=names(simu)[i],V2=m0[,2]))<br>
}<br>
class(resS)<br>
<br>
### Aproximaçao para Slope<br>
marginal.plot(resS$V2, data=resS,groups = resS$V1,auto.key=TRUE)<br>
<br>
### Diferenças<br>
mS <- aov(V2~V1,resS)<br>
summary(mS)<br>
par(mfrow=c(2,2))<br>
plot(mS)<br>
<br>
Creio que esta ANOVA no final não seja a coisa mais interessante de ser<br>
feita, diante disso gostaria de sugestões para o problema.<br>
Acho que problemas de heterogenidade de variancias vão surgir....<br>
OBS1: As outras covariaveis dos dados estão sendo discutidas se serão<br>
feitas, por enquanto desconsiderar elas<br>
Atenciosamente<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110914/6f27d5ef/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110914/6f27d5ef/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Thu, 15 Sep 2011 07:27:37 -0500<br>
From: "Alan Rodrigo Panosso" <<a href="mailto:arpanosso@yahoo.com.br">arpanosso@yahoo.com.br</a>><br>
To: <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: [R-br] RES: Seleção Níveis Regressão linear<br>
Message-ID: <000001cc73a2$dc435c40$94ca14c0$@<a href="http://yahoo.com.br" target="_blank">yahoo.com.br</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Eder, não sei se você já viu esse trabalho sobre emissão de metano.<br>
<br>
COSTA, F. S. ; BAYER, Cimélio ; LIMA, Magda Aparecida de ; FRIGHETTO, Rosa T<br>
S ; MACEDO, Vera R M ; MARCOLIN, Elio . Variação diária da emissão de metano<br>
em solo cultivado com arroz irrigado no Sul do Brasil.. Ciência Rural, v.<br>
38, p. 2049-2053, 2008.<br>
<br>
Estamos utilizando o mesmo esquema para o CO2<br>
<br>
<br>
<br>
Atenciosamente<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
De: <a href="mailto:r-br-bounces@listas.c3sl.ufpr.br">r-br-bounces@listas.c3sl.ufpr.br</a><br>
[mailto:<a href="mailto:r-br-bounces@listas.c3sl.ufpr.br">r-br-bounces@listas.c3sl.ufpr.br</a>] Em nome de Eder David Borges da<br>
Silva<br>
Enviada em: quarta-feira, 14 de setembro de 2011 20:44<br>
Para: <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
Assunto: [R-br] Seleção Níveis Regressão linear<br>
<br>
<br>
<br>
Boa noite pessoal,<br>
<br>
Gostaria de discutir com vocês um experimento.<br>
<br>
Estamos planejando um experimento de medição de gases do solo via câmaras,<br>
em síntese colocar uma câmara sob o solo e a cada intervalo de tempo medi-se<br>
a concentração de gases (CO2, N2O, CH4) retirando uma amostra e fazendo<br>
cromatografia dela.<br>
<br>
nosso experimento visa otimizar os tempos para coleta destas amostra hoje se<br>
faz 0, 15, 30, 45 min.<br>
<br>
Vamos fazer um experimento com 13 tempos de 0 a 120 min.<br>
<br>
Nosso interesse neste experimento é determinar a inclinação da reta que da a<br>
taxa de emissão (ppm / min)<br>
<br>
Minha sugestão de analise:<br>
<br>
1) Ajustar um modelo com os 13 tempos. sendo esta estimativa de inclinação a<br>
mais próxima da realidade. (M0)<br>
<br>
2) ajustar modelos com tempos distintos por exemplo 0, 5, 10 , 20 (Mi) e<br>
comparar esta estimativas com o (M0)<br>
<br>
repetir o passo 2 varias vezes<br>
<br>
3) computar os modelos que são "iguais" a M0, de todos os "iguais" o melhor<br>
é aquele que demore menos tempo.<br>
<br>
Em resumo se avaliar nos tempos 0,5,10,20 e igual a avaliar em 13 tempos de<br>
0 a 120???<br>
<br>
OBS: o tempo 0 sempre terá de ser medido<br>
<br>
Um Codigo possivel seria isso:<br>
<br>
require(MCMCpack)<br>
<br>
require(plyr)<br>
<br>
require(latticeExtra)<br>
<br>
<br>
<br>
dados <-<br>
expand.grid(local=c('Agri','Pastagem'),epoca=factor(c(-1,1,5,10,20)),camara=<br>
c('C1','C2','C3'),tempo=c(seq(0,45,by=5),60,90,120),bloco=factor(1:4))<br>
<br>
dim(dados)<br>
<br>
dados$resp <- rnorm(nrow(dados),10,10)<br>
<br>
head(dados)<br>
<br>
<br>
<br>
m0 <- MCMCregress(resp~tempo,data=dados)<br>
<br>
summary(m0)<br>
<br>
plot(m0)<br>
<br>
###Distribuiçao pretendida<br>
<br>
<br>
<br>
#### Simulação<br>
<br>
simu <- list(A=c(0,5,10,60),<br>
<br>
B=c(0,5,15,45),<br>
<br>
C=c(0,10,20,45),<br>
<br>
D=c(0,10,45,60))### colocar mais<br>
<br>
<br>
<br>
resS <- data.frame(V1='m0',V2=m0[,2])<br>
<br>
<br>
<br>
for(i in 1:length(simu)){<br>
<br>
m0 <- MCMCregress(resp~tempo,data=subset(dados,dados$tempo %in%<br>
simu[[i]]),burnin=1000,mcmc=5000)<br>
<br>
m0 <- as.matrix(m0)<br>
<br>
resS <- rbind(resS,data.frame(V1=names(simu)[i],V2=m0[,2]))<br>
<br>
}<br>
<br>
class(resS)<br>
<br>
<br>
<br>
### Aproximaçao para Slope<br>
<br>
marginal.plot(resS$V2, data=resS,groups = resS$V1,auto.key=TRUE)<br>
<br>
<br>
<br>
### Diferenças<br>
<br>
mS <- aov(V2~V1,resS)<br>
<br>
summary(mS)<br>
<br>
par(mfrow=c(2,2))<br>
<br>
plot(mS)<br>
<br>
<br>
<br>
Creio que esta ANOVA no final não seja a coisa mais interessante de ser<br>
feita, diante disso gostaria de sugestões para o problema.<br>
<br>
Acho que problemas de heterogenidade de variancias vão surgir....<br>
<br>
OBS1: As outras covariaveis dos dados estão sendo discutidas se serão<br>
feitas, por enquanto desconsiderar elas<br>
<br>
Atenciosamente<br>
<br>
<br>
<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110915/9501bf3e/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110915/9501bf3e/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Thu, 15 Sep 2011 11:32:19 -0300<br>
From: Walmes Zeviani <<a href="mailto:walmeszeviani@gmail.com">walmeszeviani@gmail.com</a>><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
Subject: Re: [R-br] RES: Seleção Níveis Regressão linear<br>
Message-ID:<br>
<CAFU=<a href="mailto:EkY37yZx2ARCrwuj83Lbf_0HZbRCVK0gZ4PBr4WCwNAKYg@mail.gmail.com">EkY37yZx2ARCrwuj83Lbf_0HZbRCVK0gZ4PBr4WCwNAKYg@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Éder,<br>
<br>
Isso me parece uma questão de delineamento ótimo, que no caso linear<br>
apresenta solução simples. Como você já sabe qual o modelo que representa os<br>
dados, é só otimizar a escolha dos tempos para obter um estimador com menor<br>
erro-padrão. Seria o caso então de estudar os critérios de otimabilidade que<br>
existem para modelos lineares.<br>
<br>
À disposição.<br>
Walmes.<br>
<br>
==========================================================================<br>
Walmes Marques Zeviani<br>
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)<br>
Departamento de Estatística - Universidade Federal do Paraná<br>
fone: <a href="tel:%28%2B55%29%2041%203361%203573" value="+554133613573">(+55) 41 3361 3573</a><br>
VoIP: (3361 3600) 1053 1173<br>
e-mail: <a href="mailto:walmes@ufpr.br">walmes@ufpr.br</a><br>
twitter: @walmeszeviani<br>
homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a><br>
linux user number: 531218<br>
==========================================================================<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110915/6781ee46/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110915/6781ee46/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
<br>
<br>
Fim da Digest R-br, volume 7, assunto 16<br>
****************************************<br>
</blockquote></div><br><br clear="all"><br>-- <br>Júlio Sílvio de Sousa Bueno Filho<br><br>Departamento de Ciências Exatas <br>Universidade Federal de Lavras<br>Caixa Postal 3037<br>Lavras, MG - Brasil - 37200 000<br>tel: (+55) 35 3829 1369<br>
fax: (+55) 35 3829 1371<br><br> "Natureza da gente não cabe em nenhuma certeza."<br>João G. Rosa - Riobaldo (Grande Sertão: Veredas)<br><br>"Homo sum; humani nihil a me alienum puto."<br>Terêncio - Chremes (Heautontimoroumenos)<br>