<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix"><br>
Tiago,<br>
<br>
veja uma simulação sem a correção dos dados pela formula de Abbot:<br>
<br>
<br>
conc = c(0.15,0.25,0.41,0.67,1.11,1.83,3.03,5)<br>
ni = c(53,51,48,45,50,50,48,50)<br>
im = c(7,15,28,31,42,49,48,50)<br>
<br>
mod=drm(im/ni~conc,weight=ni,fct=LL.2(),type="binomial")<br>
modelFit(mod)<br>
ED(mod,c(50),interval="delta")<br>
<br>
tratVec = with(mod, seq(min(conc), max(conc), length.out=100))<br>
Pred=predict(mod,newdata=data.frame(conc=tratVec,1),type="response",se.fit=TRUE)<br>
Li=Pred[,1]-2*Pred[,2]<br>
Ls=Pred[,1]+2*Pred[,2]<br>
<br>
plot(mod,conName="0",axes=FALSE,type="average",xlim=c(0.15,5),ylim=c(0,1),log="x")<br>
axis(1,conc)<br>
y=seq(0,1,0.1)<br>
axis(2,y)<br>
<br>
lines(tratVec, Li, lty = 2)<br>
lines(tratVec, Ls, lty = 2)<br>
<br>
<br>
<br>
Stephan<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
Em 11/03/2013 16:42, Tiago Souza Marçal escreveu:<br>
</div>
<blockquote cite="mid:BLU170-W75CA8084F33B0FFFAA7385D9E10@phx.gbl"
type="cite">
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
<div dir="ltr"> Se eu corrigir o meus dados pelo critério de
Abbott qual o modelo você me aconcelha a utilizar (LL.2()) ? <br>
<br>
CMR:<br>
<br>
dados<-structure(list(conc = c(0.15, 0.25, 0.41, 0.67, 1.11,
1.83, 3.03, <br>
5), lconc = c(-0.823908741, -0.602059991, -0.387216143,
-0.173925197, <br>
0.045322979, 0.26245109, 0.481442629, 0.698970004), ni = c(53L,
<br>
51L, 48L, 45L, 50L, 50L, 48L, 50L), im = c(7L, 15L, 28L, 31L, <br>
42L, 49L, 48L, 50L), m = c(0.132075472, 0.294117647,
0.583333333, <br>
0.688888889, 0.84, 0.98, 1, 1), mc = c(0.053686206, 0.215855573,
<br>
0.505298273, 0.610936682, 0.762166405, 0.902276295, 0.922291994,
<br>
0.922291994)), .Names = c("conc", "lconc", "ni", "im", "m", "mc"<br>
), class = "data.frame", row.names = c(NA, -8L))<br>
<br>
Onde:<br>
<br>
conc = concentrações<br>
lconc = logaritimo neperiano das concentrações.<br>
ni = número de individúos<br>
im = individúos mortos<br>
m = mortalidade<br>
mc = mortalidade corrigida<br>
<br>
Att.<br>
<br>
Tiago. <br>
<div>
<hr id="stopSpelling">Date: Mon, 11 Mar 2013 13:38:30 -0300<br>
From: <a class="moz-txt-link-abbreviated"
href="mailto:smalfitano@uol.com.br">smalfitano@uol.com.br</a><br>
To: <a class="moz-txt-link-abbreviated"
href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
Subject: Re: [R-br] Pacote drc<br>
<br>
<div class="ecxmoz-cite-prefix"><br>
Bom dia Tiago,<br>
<br>
O drc assume o controle como "dose zero", sem fazer
nenhuma correção como proposto por Abbott.<br>
<br>
Além disso é preciso verificar qual tipo de distribuição
(e parâmetros) que você está usando, pois existe uma grande
diversidade no drc quando comparado por exemplo com o Polo.<br>
<br>
Att,<br>
<br>
Stephan<br>
<br>
<br>
Em 11/03/2013 13:26, Tiago Souza Marçal escreveu:<br>
</div>
<blockquote
cite="mid:BLU170-W943D474EA7A9415361CB64D9E10@phx.gbl">
<style><!--
.ExternalClass .ecxhmmessage P
{padding:0px;}
.ExternalClass body.ecxhmmessage
{font-size:10pt;font-family:Tahoma;}
--></style>
<div dir="ltr"> Bom dia pessoal,<br>
<br>
gostaria de saber se alguem da lista utiliza o pacote
drc para estimar a CL-50. Caso tenha gostaria que me
esclarecesse as seguintes dúvidas:<br>
<br>
Ele calcula a mortalidade corrigida pelo critério de
Abbott?<br>
<br>
Porque as estimativas geradas por este pacote são
diferentes das geradas pelo SAS e o Polo?<br>
<br>
Att.<br>
<br>
Tiago.<br>
</div>
<br>
<fieldset class="ecxmimeAttachmentHeader"></fieldset>
<br>
<pre>_______________________________________________
R-br mailing list
<a moz-do-not-send="true" class="ecxmoz-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="ecxmoz-txt-link-freetext" 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>
Leia o guia de postagem (<a moz-do-not-send="true" class="ecxmoz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
<br>
<pre class="ecxmoz-signature">--
***************************************************
STEPHAN M. CARVALHO
"Chemical control and ecotoxicology"
Universidade Federal de Uberlândia
Instituto de Ciências Agrárias
Avenida Amazonas s/n – Bloco 2E
Caixa Postal 593 - Umuarama
38.400-902 Uberlândia/MG
Brazil
<a moz-do-not-send="true" class="ecxmoz-txt-link-abbreviated" href="mailto:stephan.carvalho@iciag.ufu.br">stephan.carvalho@iciag.ufu.br</a>
+55 34 3218-2225 (university)
+55 34 3218-2225 (fax)
+55 35 9103-9822 (mobile)
+55 35 9949-0707 (mobile)
Antes de imprimir, pense no ambiente!
Before printing, think in the environment!
Avant d'imprimer, pensez à l'environnement!
***************************************************
</pre>
<br>
_______________________________________________ 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.</div>
</div>
<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>
<br>
<pre class="moz-signature" cols="72">--
***************************************************
STEPHAN M. CARVALHO
"Chemical control and ecotoxicology"
Universidade Federal de Uberlândia
Instituto de Ciências Agrárias
Avenida Amazonas s/n – Bloco 2E
Caixa Postal 593 - Umuarama
38.400-902 Uberlândia/MG
Brazil
<a class="moz-txt-link-abbreviated" href="mailto:stephan.carvalho@iciag.ufu.br">stephan.carvalho@iciag.ufu.br</a>
+55 34 3218-2225 (university)
+55 34 3218-2225 (fax)
+55 35 9103-9822 (mobile)
+55 35 9949-0707 (mobile)
Antes de imprimir, pense no ambiente!
Before printing, think in the environment!
Avant d'imprimer, pensez à l'environnement!
***************************************************
</pre>
</body>
</html>