<div dir="ltr"><div><div><div><div><div><div>Prezados,<br><br></div>sou biólogo e um dos objetivos da minha pesquisa é examinar os efeitos da heterogeneidade do ambiente na variabilidade genética. Para isso, eu estimei a distância genética, geodésica e ecológica e usei uma função do R - <a href="http://datadryad.org/resource/doi:10.5061/dryad.kt71r/1">http://datadryad.org/resource/doi:10.5061/dryad.kt71r/1</a> - para testar a regressão múltipla dessas matrizes.<br>
<br></div>O script para realizar o teste eu consegui fazer, entretanto, não encontrei uma função para fazer o gráfico da regressão linear sobre meus dados. O único gráfico que consigo fazer é entre duas variáveis é usando a função splom (matriz1, matriz2, lm=TRUE) e obtive uma resposta de erro (lm object not found). Mas, o que eu queria era fazer um gráfico de uma matriz1 x matriz2+matriz3.<br>
<br></div>Alguém tem uma sugestão de função que eu possa usar?<br><br></div>Agradeço a todos desde já.<br><br></div>Cordialmente,<br><br></div>Marcos Batista<br></div><div class="gmail_extra"><br clear="all"><div><div><span style="font-size:13px"><font face="verdana, sans-serif" color="#000099">Marcos Roberto Dias Batista<br>
Departamento de Genética, Evolução e Bioagentes<br>Instituto de Biologia, Unicamp<br>Caixa Postal 6109, Campinas, São Paulo<br>CEP 13083-970, Brasil<br>Telefone: (19) 3521 1151</font><br></span></div></div>
<br><br><div class="gmail_quote">2014-03-21 12:00 GMT-03:00 <span dir="ltr"><<a href="mailto:r-br-request@listas.c3sl.ufpr.br" target="_blank">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. Re: Calendário do pacote googleVis (Éder Comunello)<br>
2. Re: Grárico de controle (Edson Lira)<br>
3. Iteração com optim (jborssoi .)<br>
4. Re: Iteração com optim (Diogo Ferrari)<br>
5. Re: Grárico de controle (Éder Comunello)<br>
6. Re: Iteração com optim (walmes .)<br>
7. Re: Iteração com optim (Wagner Bonat)<br>
8. Re: Iteração com optim (jborssoi .)<br>
9. Re: Iteração com optim (Wagner Bonat)<br>
10. Re: Grárico de controle (Edson Lira)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Fri, 21 Mar 2014 07:17:25 -0400<br>
From: Éder Comunello <<a href="mailto:comunello.eder@gmail.com">comunello.eder@gmail.com</a>><br>
To: "<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>" <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: Re: [R-br] Calendário do pacote googleVis<br>
Message-ID:<br>
<CABmC8gnNgPU2mWLFczTUrWYsqOf2A9OjKv5dpF7UnNM=<a href="mailto:aTie6A@mail.gmail.com">aTie6A@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Olá, pessoal! Bom dia!<br>
<br>
Tudo indica que o problema está na API do Google e não propriamente no R. O<br>
link que segue relata o mesmo problema numa aplicação fora do R, feita<br>
também por um brasileiro, reforçando que o problema está ligado ao 'Brazilian<br>
Portuguese'.<br>
<br>
<a href="https://code.google.com/p/google-visualization-api-issues/issues/detail?id=1477" target="_blank">https://code.google.com/p/google-visualization-api-issues/issues/detail?id=1477</a><br>
<br>
Éder Comunello <c <<a href="mailto:comunello.eder@gmail.com">comunello.eder@gmail.com</a>><a href="mailto:omunello.eder@gmail.com">omunello.eder@gmail.com</a>><br>
Dourados, MS - [22 16.5'S, 54 49'W]<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/a91fe90a/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/a91fe90a/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Fri, 21 Mar 2014 05:47:41 -0700 (PDT)<br>
From: Edson Lira <<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>><br>
To: "<a href="mailto:ecomunel@gmail.com">ecomunel@gmail.com</a>" <<a href="mailto:ecomunel@gmail.com">ecomunel@gmail.com</a>>,<br>
"<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>" <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: Re: [R-br] Grárico de controle<br>
Message-ID:<br>
<<a href="mailto:1395406061.56210.YahooMailNeo@web160103.mail.bf1.yahoo.com">1395406061.56210.YahooMailNeo@web160103.mail.bf1.yahoo.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Eder, muito obrigado, ajudou muito, vi que tenho no mesmo dia vários grupos.<br>
<br>
Como faço para ter um único resultado por dia? Ou seja, para cada dia um único ponto no gráfico.<br>
<br>
<br>
[ ]'s.<br>
<br>
Edson Lira<br>
Estatístico<br>
Manaus-Amazonas<br>
<br>
<br>
<br>
Em Terça-feira, 18 de Março de 2014 12:52, Éder Comunello <<a href="mailto:comunello.eder@gmail.com">comunello.eder@gmail.com</a>> escreveu:<br>
<br>
Edson, boa tarde!<br>
<br>
Não sei se irá lhe atender, mas tentei alguma coisa com os dados de exemplo que você havia enviado (n=300). <br>
<br>
### <code r><br>
### restaurando dados de exemplo!<br>
dat <-<br>
data.frame(tempo=c(30.48,32.17,32.32,32.68,33.42,34.1,34.33,35.12,35.95,36.18,36.73,36.75,36.78,37.05,37.18,38.1,38.4,40.05,40.62,40.87,41.33,41.38,41.5,41.63,43.38,43.65,43.85,43.88,44.37,45.67,45.73,46.17,46.5,46.75,47.17,47.55,47.63,47.65,49.62,49.85,49.98,50.65,51.15,51.67,51.7,51.97,53.12,53.23,53.77,54.08,54.18,54.32,54.65,54.98,55.55,56.47,58.13,58.65,58.93,59.2,59.67,59.92,60.08,60.28,60.63,60.93,61.78,62.58,63.42,63.93,64.13,65.12,66.02,66.55,67.23,67.85,71.32,73.03,73.68,75.1,75.52,79.75,79.93,80.88,81.23,91.85,103.73,109.27,115.65,123.63,131.88,173.28,30.68,31.68,32.1,33.12,33.25,33.48,34.5,35.15,35.52,35.87,36.52,36.63,36.75,39.32,40.37,41.4,41.57,41.57,41.78,43.85,44.05,44.97,45.47,45.93,47.03,47.1,47.35,47.63,48.3,48.38,48.62,48.85,48.88,49,49.87,50.65,50.73,51.07,51.42,51.8,52.08,52.72,52.83,53.23,53.32,53.57,53.68,54,54.12,54.23,54.28,54.5,55.13,55.32,55.65,55.85,56.02,56.15,56.43,56.45,56.58,56.58,57.22,57.77,57.78,58.17,58.43,58.57,58<br>
.58,59.38,59.8,59.87,60.18,60.67,61.35,61.42,61.6,62.07,62.18,62.42,63.07,63.93,64.57,65.2,65.33,65.78,65.8,66.38,66.55,66.78,67.07,67.43,67.57,68.28,68.3,68.32,68.43,68.68,68.82,69,69.55,70.22,70.58,71.42,71.6,72.13,72.17,72.55,72.57,72.67,72.7,72.8,72.8,73.38,73.38,73.67,74.28,75.05,75.4,75.53,75.75,76.6,77.75,78.22,78.42,78.47,78.63,79.27,79.58,79.78,80.18,81.53,81.78,82.17,84.87,85.07,85.2,85.77,86.37,87.25,88.78,88.93,88.98,89.68,91,92.62,93.5,93.83,95.12,97.1,97.23,97.77,98.83,101.47,107.87,108.17,109.07,112.25,112.58,113.53,113.75,119.48,121.32,132.88,139.73,31.95,32.4,35.2,35.35,36.62,36.88,37.62,38.65,39.68,40.05,41.25,41.28,42,42.03,42.23,42.63,42.67,43.02,43.32,43.32,44.23,44.27,44.7,45.08,46.42,48.17,48.25,49.55,50.13,50.85,50.97,51.85,52.28,52.85,52.9,53.13,54.23,54.98,55.45,56.28,56.53,57.33,57.75), dia_doa=rep(c(1:3), c(92,165,43)))<br>
<br>
table(dat$dia_doa)<br>
<br>
### (r)andom (Pos)itions (dentro de cada dia)<br>
{set.seed(765); rPos <- unlist(sapply(unique(dat$dia_doa), function(x) sample(which(dat$dia_doa==x))))}<br>
<br>
### set.seed() é pra controle do CMR, considere removê-lo<br>
<br>
rDat <- dat[rPos,]; head(rDat)<br>
rDat$dia_doa ### observe que não muda ordem dos dias!<br>
<br>
dN <- unique(table(dat$dia_doa)) ### dados em cada dia!<br>
grp <- NULL; for (n in dN) grp <- c(grp, rep(1:ceiling(n/5), each=5, len=n)) ### grupos dentro de dias<br>
<br>
DAT <- cbind(rDat, grp); row.names(DAT) <- NULL; head(DAT)<br>
DAT$id <- DAT$dia_doa*1000+DAT$grp ### cria um id/label único para grupos!<br>
<br>
### Caso deseje remover dados/grupos incompletos...<br>
grp.i <- names(table(DAT$id))[table(DAT$id)<5];grp.i ### grupos incompletos<br>
dat.i <- unlist(sapply(grp.i, function(x) which(DAT$id==x))); dat.i ### dados pertencentes a grupos incompletos<br>
<br>
DAT[dat.i,] ### linhas dos grupos incompletos<br>
DATCompl <- DAT[-dat.i,] ### só grupos completos<br>
table(DATCompl$dia_doa)<br>
table(DATCompl$id)<br>
<br>
###<br>
<br>
require(qcc)<br>
qcc(qcc.groups(DAT$tempo, DAT$id)) ### todos<br>
qcc(qcc.groups(DATCompl$tempo, DATCompl$id)) ### somente completos<br>
### </code><br>
<br>
<br>
<br>
<br>
Éder Comunello <<a href="mailto:comunello.eder@gmail.com">comunello.eder@gmail.com</a>> <br>
<br>
Dourados, MS - [22 16.5'S, 54 49'W]<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/451376a0/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/451376a0/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Fri, 21 Mar 2014 10:26:15 -0300<br>
From: "jborssoi ." <<a href="mailto:joelmirb@gmail.com">joelmirb@gmail.com</a>><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
Subject: [R-br] Iteração com optim<br>
Message-ID:<br>
<<a href="mailto:CAEFzrQsPbZ_DxKqC5w3iSYP7851UZaqqCdTUUikwBTnVeSShPA@mail.gmail.com">CAEFzrQsPbZ_DxKqC5w3iSYP7851UZaqqCdTUUikwBTnVeSShPA@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Caros,<br>
tenho 100 amostras e preciso fazer a estimação de parâmetros usando<br>
a função optim. Como faço para que, a cada nova estimação, o vetor com o<br>
chute inicial ("par") seja atualizado com as estimativas obtidas no passo<br>
anterior?<br>
<br>
Grato pela ajuda<br>
Joelmir<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/85bdad8e/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/85bdad8e/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Fri, 21 Mar 2014 09:55:59 -0400<br>
From: Diogo Ferrari <<a href="mailto:diogoferrari@gmail.com">diogoferrari@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] Iteração com optim<br>
Message-ID:<br>
<<a href="mailto:CAEvjS8Lwo0C5DuQUisvNfATUJ6Zy3Qcao_pa_UxM5jA5G_2eDg@mail.gmail.com">CAEvjS8Lwo0C5DuQUisvNfATUJ6Zy3Qcao_pa_UxM5jA5G_2eDg@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Nao sei se compreendi bem sua questao, mas acredito q optim faz exatamente<br>
isso automaticamente.<br>
On Mar 21, 2014 9:26 AM, "jborssoi ." <<a href="mailto:joelmirb@gmail.com">joelmirb@gmail.com</a>> wrote:<br>
<br>
> Caros,<br>
> tenho 100 amostras e preciso fazer a estimação de parâmetros usando<br>
> a função optim. Como faço para que, a cada nova estimação, o vetor com o<br>
> chute inicial ("par") seja atualizado com as estimativas obtidas no passo<br>
> anterior?<br>
><br>
> Grato pela ajuda<br>
> Joelmir<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>
> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
> código mínimo reproduzível.<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/20140321/7fdd0eb1/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/7fdd0eb1/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Fri, 21 Mar 2014 10:09:24 -0400<br>
From: Éder Comunello <<a href="mailto:comunello.eder@gmail.com">comunello.eder@gmail.com</a>><br>
To: "<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>" <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>>, Edson Lira<br>
<<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>><br>
Subject: Re: [R-br] Grárico de controle<br>
Message-ID:<br>
<<a href="mailto:CABmC8gkK_HxJ7wGUqsUhounM84TPkQgzd6ZHcA7O0VyPWPzqAA@mail.gmail.com">CABmC8gkK_HxJ7wGUqsUhounM84TPkQgzd6ZHcA7O0VyPWPzqAA@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Edson, bom dia!<br>
<br>
Um ponto por dia simplifica bastante o procedimento e você pode dispensar<br>
praticamente todo o código anteriormente enviado.<br>
<br>
Tendo criado (ou restaurado) o objeto de dados 'dat', utilize:<br>
<br>
### qcc.groups(data, sample)<br>
qcc.groups(dat$tempo, dat$dia_doa)<br>
qcc(qcc.groups(dat$tempo, dat$dia_doa))<br>
<br>
Reforçando que serão amostras de tamanhos distintos.<br>
<br>
Éder Comunello <c <<a href="mailto:comunello.eder@gmail.com">comunello.eder@gmail.com</a>><a href="mailto:omunello.eder@gmail.com">omunello.eder@gmail.com</a>><br>
Dourados, MS - [22 16.5'S, 54 49'W]<br>
<br>
<br>
Em 21 de março de 2014 08:47, Edson Lira <<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>>escreveu:<br>
<br>
> Eder, muito obrigado, ajudou muito, vi que tenho no mesmo dia vários<br>
> grupos.<br>
><br>
> Como faço para ter um único resultado por dia? Ou seja, para cada dia um<br>
> único ponto no gráfico.<br>
><br>
><br>
> [ ]'s.<br>
> Edson Lira<br>
> Estatístico<br>
> Manaus-Amazonas<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/20140321/8124fe7c/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/8124fe7c/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 6<br>
Date: Fri, 21 Mar 2014 11:15:50 -0300<br>
From: "walmes ." <<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>" <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: Re: [R-br] Iteração com optim<br>
Message-ID:<br>
<CAFU=<a href="mailto:EkZMw34AABsMzb9r8Ey5KkbZi6rg1UFgpmBcc0u03wfzUA@mail.gmail.com">EkZMw34AABsMzb9r8Ey5KkbZi6rg1UFgpmBcc0u03wfzUA@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Faça com um laço for por exemplo<br>
<br>
chute <- ...<br>
for(...){<br>
op <- optim(chute, ...)<br>
chute <- op$par<br>
}<br>
<br>
À disposição.<br>
Walmes.<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/20140321/d078d6db/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/d078d6db/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 7<br>
Date: Fri, 21 Mar 2014 15:16:27 +0100<br>
From: Wagner Bonat <<a href="mailto:wbonat@gmail.com">wbonat@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] Iteração com optim<br>
Message-ID:<br>
<CANt=<a href="mailto:4MhG-eB4uOTxSCzmJrdm1aEbpmXa3LF6ax3c_OsaDAfgdw@mail.gmail.com">4MhG-eB4uOTxSCzmJrdm1aEbpmXa3LF6ax3c_OsaDAfgdw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
O que vc quer dizer com 100 amostras ?? Vc tem 100 conjuntos de dados<br>
diferentes e quer estimar para cada um determinado modelo ... Nao consegui<br>
entender ... Que metodo vc quer usar para inferencia ??<br>
<br>
<br>
Em 21 de março de 2014 14:55, Diogo Ferrari <<a href="mailto:diogoferrari@gmail.com">diogoferrari@gmail.com</a>>escreveu:<br>
<br>
> Nao sei se compreendi bem sua questao, mas acredito q optim faz exatamente<br>
> isso automaticamente.<br>
> On Mar 21, 2014 9:26 AM, "jborssoi ." <<a href="mailto:joelmirb@gmail.com">joelmirb@gmail.com</a>> wrote:<br>
><br>
>> Caros,<br>
>> tenho 100 amostras e preciso fazer a estimação de parâmetros usando<br>
>> a função optim. Como faço para que, a cada nova estimação, o vetor com o<br>
>> chute inicial ("par") seja atualizado com as estimativas obtidas no passo<br>
>> anterior?<br>
>><br>
>> Grato pela ajuda<br>
>> Joelmir<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>
>> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
>> código mínimo reproduzível.<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>
> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
> código mínimo reproduzível.<br>
><br>
<br>
<br>
<br>
--<br>
Wagner Hugo Bonat<br>
LEG - Laboratório de Estatística e Geoinformação<br>
UFPR - Universidade Federal do Paraná<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/a5d6012f/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/a5d6012f/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 8<br>
Date: Fri, 21 Mar 2014 11:23:04 -0300<br>
From: "jborssoi ." <<a href="mailto:joelmirb@gmail.com">joelmirb@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] Iteração com optim<br>
Message-ID:<br>
<<a href="mailto:CAEFzrQuNS%2Boj9FmamJJkmUtp9zp8rYhHXknMdawxnxcmRE151A@mail.gmail.com">CAEFzrQuNS+oj9FmamJJkmUtp9zp8rYhHXknMdawxnxcmRE151A@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Olá Wagner,<br>
é isso... tenho 100 conjuntos de dados. Eu preciso maximizar o vetor de<br>
parâmetros para a primeira amostra. Já para a segunda (e assim por diante)<br>
o optim deve utilizar o vetor obtido para a primeira amostra (e assim por<br>
diante).<br>
Imagino algo do tipo:<br>
for(i in 1:J){<br>
maxL.rep <- optim(par=???, fn, gr, method="BFGS",<br>
control=list(fnscale=-1,maxit=1000),y=W1$y[,i], u=W1$u[,i], X=d[,c(1,2)],<br>
Z=Zi, sigma.e=1)<br>
}<br>
<br>
mas não sei como fazer para que o optim atualize "par".<br>
<br>
abraço<br>
<br>
<br>
Em 21 de março de 2014 11:16, Wagner Bonat <<a href="mailto:wbonat@gmail.com">wbonat@gmail.com</a>> escreveu:<br>
<br>
> O que vc quer dizer com 100 amostras ?? Vc tem 100 conjuntos de dados<br>
> diferentes e quer estimar para cada um determinado modelo ... Nao consegui<br>
> entender ... Que metodo vc quer usar para inferencia ??<br>
><br>
><br>
> Em 21 de março de 2014 14:55, Diogo Ferrari <<a href="mailto:diogoferrari@gmail.com">diogoferrari@gmail.com</a>>escreveu:<br>
><br>
> Nao sei se compreendi bem sua questao, mas acredito q optim faz exatamente<br>
>> isso automaticamente.<br>
>> On Mar 21, 2014 9:26 AM, "jborssoi ." <<a href="mailto:joelmirb@gmail.com">joelmirb@gmail.com</a>> wrote:<br>
>><br>
>>> Caros,<br>
>>> tenho 100 amostras e preciso fazer a estimação de parâmetros usando<br>
>>> a função optim. Como faço para que, a cada nova estimação, o vetor com o<br>
>>> chute inicial ("par") seja atualizado com as estimativas obtidas no passo<br>
>>> anterior?<br>
>>><br>
>>> Grato pela ajuda<br>
>>> Joelmir<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>
>>> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
>>> código mínimo reproduzível.<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>
>> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
>> código mínimo reproduzível.<br>
>><br>
><br>
><br>
><br>
> --<br>
> Wagner Hugo Bonat<br>
> LEG - Laboratório de Estatística e Geoinformação<br>
> UFPR - Universidade Federal do Paraná<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>
> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
> código mínimo reproduzível.<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/20140321/94ab16c2/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/94ab16c2/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 9<br>
Date: Fri, 21 Mar 2014 15:36:43 +0100<br>
From: Wagner Bonat <<a href="mailto:wbonat@gmail.com">wbonat@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] Iteração com optim<br>
Message-ID:<br>
<CANt=<a href="mailto:4Mh5gG2QAodJvBDv%2B3R9vZs3_wf%2BqRnGZ9Bo0iV-NjjYhw@mail.gmail.com">4Mh5gG2QAodJvBDv+3R9vZs3_wf+qRnGZ9Bo0iV-NjjYhw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Talvez isso ajude !<br>
<br>
ll <- function(mu,x){-sum(dnorm(x,mean=mu, sd=1,log=TRUE))}<br>
ll(mu = 0, x = rnorm(10))<br>
dados <- matrix(rnorm(1000), ncol=100, nrow=10)<br>
saida <- list()<br>
saida[1][[1]] <- optim(par= 0, f = ll, x = dados[,1], method="Brent", upper<br>
= 10, lower = -10)<br>
saida[1][[1]]$par<br>
for(i in 2:100){<br>
print(saida[c(i-1)][[1]]$par)<br>
saida[i][[1]] <- optim(par= saida[c(i-1)][[1]]$par, f = ll, x =<br>
dados[,i], method="Brent", upper = 10, lower = -10)<br>
}<br>
<br>
<br>
<br>
Em 21 de março de 2014 15:23, jborssoi . <<a href="mailto:joelmirb@gmail.com">joelmirb@gmail.com</a>> escreveu:<br>
<br>
> Olá Wagner,<br>
> é isso... tenho 100 conjuntos de dados. Eu preciso maximizar o vetor de<br>
> parâmetros para a primeira amostra. Já para a segunda (e assim por diante)<br>
> o optim deve utilizar o vetor obtido para a primeira amostra (e assim por<br>
> diante).<br>
> Imagino algo do tipo:<br>
> for(i in 1:J){<br>
> maxL.rep <- optim(par=???, fn, gr, method="BFGS",<br>
> control=list(fnscale=-1,maxit=1000),y=W1$y[,i], u=W1$u[,i], X=d[,c(1,2)],<br>
> Z=Zi, sigma.e=1)<br>
> }<br>
><br>
> mas não sei como fazer para que o optim atualize "par".<br>
><br>
> abraço<br>
><br>
><br>
> Em 21 de março de 2014 11:16, Wagner Bonat <<a href="mailto:wbonat@gmail.com">wbonat@gmail.com</a>> escreveu:<br>
><br>
> O que vc quer dizer com 100 amostras ?? Vc tem 100 conjuntos de dados<br>
>> diferentes e quer estimar para cada um determinado modelo ... Nao consegui<br>
>> entender ... Que metodo vc quer usar para inferencia ??<br>
>><br>
>><br>
>> Em 21 de março de 2014 14:55, Diogo Ferrari <<a href="mailto:diogoferrari@gmail.com">diogoferrari@gmail.com</a>>escreveu:<br>
>><br>
>> Nao sei se compreendi bem sua questao, mas acredito q optim faz<br>
>>> exatamente isso automaticamente.<br>
>>> On Mar 21, 2014 9:26 AM, "jborssoi ." <<a href="mailto:joelmirb@gmail.com">joelmirb@gmail.com</a>> wrote:<br>
>>><br>
>>>> Caros,<br>
>>>> tenho 100 amostras e preciso fazer a estimação de parâmetros usando<br>
>>>> a função optim. Como faço para que, a cada nova estimação, o vetor com o<br>
>>>> chute inicial ("par") seja atualizado com as estimativas obtidas no passo<br>
>>>> anterior?<br>
>>>><br>
>>>> Grato pela ajuda<br>
>>>> Joelmir<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>
>>>> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
>>>> código mínimo reproduzível.<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>
>>> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
>>> código mínimo reproduzível.<br>
>>><br>
>><br>
>><br>
>><br>
>> --<br>
>> Wagner Hugo Bonat<br>
>> LEG - Laboratório de Estatística e Geoinformação<br>
>> UFPR - Universidade Federal do Paraná<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>
>> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
>> código mínimo reproduzível.<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>
> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
> código mínimo reproduzível.<br>
><br>
<br>
<br>
<br>
--<br>
Wagner Hugo Bonat<br>
LEG - Laboratório de Estatística e Geoinformação<br>
UFPR - Universidade Federal do Paraná<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/e0e422aa/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/e0e422aa/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 10<br>
Date: Fri, 21 Mar 2014 07:57:30 -0700 (PDT)<br>
From: Edson Lira <<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>><br>
To: "<a href="mailto:ecomunel@gmail.com">ecomunel@gmail.com</a>" <<a href="mailto:ecomunel@gmail.com">ecomunel@gmail.com</a>>,<br>
"<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>" <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: Re: [R-br] Grárico de controle<br>
Message-ID:<br>
<<a href="mailto:1395413850.17732.YahooMailNeo@web160101.mail.bf1.yahoo.com">1395413850.17732.YahooMailNeo@web160101.mail.bf1.yahoo.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Era isso mesmo que precisava.<br>
<br>
Agora tenho as seguintes perguntas:<br>
<br>
Como faço para fixar os limites em 2 sigmas, usando o comando que vi no help não funcionou, não sei por que.<br>
Quando uso size ele fixa, mas não me permite mexer no nsigma. tem algo que não estou sabendo fazer.<br>
<br>
Veja se você consegue alguma sugestão.<br>
<br>
Muito obrigado pela sua atenção Eder.<br>
<br>
<br>
[ ]'s.<br>
<br>
Edson Lira<br>
Estatístico<br>
Manaus-Amazonas<br>
<br>
<br>
<br>
Em Sexta-feira, 21 de Março de 2014 10:09, Éder Comunello <<a href="mailto:comunello.eder@gmail.com">comunello.eder@gmail.com</a>> escreveu:<br>
<br>
Edson, bom dia!<br>
<br>
Um ponto por dia simplifica bastante o procedimento e você pode dispensar praticamente todo o código anteriormente enviado.<br>
<br>
Tendo criado (ou restaurado) o objeto de dados 'dat', utilize:<br>
<br>
### qcc.groups(data, sample)<br>
qcc.groups(dat$tempo, dat$dia_doa)<br>
qcc(qcc.groups(dat$tempo, dat$dia_doa))<br>
<br>
Reforçando que serão amostras de tamanhos distintos.<br>
<br>
<br>
Éder Comunello <<a href="mailto:comunello.eder@gmail.com">comunello.eder@gmail.com</a>><br>
Dourados, MS - [22 16.5'S, 54 49'W]<br>
<br>
<br>
<br>
Em 21 de março de 2014 08:47, Edson Lira <<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>> escreveu:<br>
<br>
Eder, muito obrigado, ajudou muito, vi que tenho no mesmo dia vários grupos.<br>
><br>
>Como faço para ter um único resultado por dia? Ou seja, para cada dia um único ponto no gráfico.<br>
><br>
><br>
><br>
><br>
>[ ]'s.<br>
><br>
>Edson Lira<br>
>Estatístico<br>
>Manaus-Amazonas<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/20140321/3cd59f3c/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/3cd59f3c/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Legenda do Digest<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>
------------------------------<br>
<br>
Fim da Digest R-br, volume 39, assunto 21<br>
*****************************************<br>
</blockquote></div><br></div>