[R-br] Problema ao truncar PDF

Cristiano Strieder cstrieder em gmail.com
Quarta Junho 15 23:41:34 BRT 2011


Obrigado pelas respostas srs, mas acho q ainda n compreendi como
solucionar o problema. Mesmo se eu passar o trunc na chamada da
função, como abaixo, ainda retorna o erro, Error in -b * x :
non-numeric argument to binary operator.

Estou me basendo neste material,
http://ecologia.ib.usp.br/let/doku.php?id=tutoriais:pdf.

dricker =  function (x,b,...) {b^2*x*exp(-b*x)}
pricker = function (x, b,...){exp(-b*x)*(exp(b*x)-b*x-1)} # CDF
dricker.t = function (x, b, ...){
  p.trunc = pricker(trunc,b)
  dricker(x,b)/(1 - p.trunc)
}
dricker.t(x=10, b=0.1, trunc=1)


2011/6/15 Paulo Justiniano <paulojus em leg.ufpr.br>:
> veja trecho marcado
> a variavel "trunc" nao eh definida em lugar algum
> e portanto a funcao nao sabe de onde ler trunc
>
> note que seu codigo pode estar lendo variaveis do GlobalEnv
>
> On Wed, 15 Jun 2011, Cristiano Strieder wrote:
>
>> Olá,
>>
>> estou com problemas no seguinte código:
>>
>> dricker =  function (x,b,...) {b^2*x*exp(-b*x)}
>> pricker = function (x, b,...){exp(-b*x)*(exp(b*x)-b*x-1)} # CDF
>> dricker.t = function (x, b,...){
>>  p.trunc = pricker(trunc,b)
>
>                     ^^^^^
>>
>>        dricker(x,b)/(1 - p.trunc)}
>> dricker.t(10, b=0.1)
>>
>
>
>
>
>> Error in -b * x : non-numeric argument to binary operator
>>
>> Alguma sugestão?
>>
>> Obrigado.
>>
>> --
>> Cristiano Strieder
>> São José dos Campos - SP
>> _______________________________________________
>> R-br mailing list
>> R-br em listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>> código mínimo reproduzível.
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
> mínimo reproduzível.
>



-- 
Cristiano Strieder
São José dos Campos - SP


Mais detalhes sobre a lista de discussão R-br