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