<html><head></head><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue-Light, Helvetica Neue Light, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:16px"><div id="yui_3_16_0_ym19_1_1487706532585_14845">Prezados, meu objetivo final é calcular valores esperados qq de uma distribuição truncada.<br></div><div id="yui_3_16_0_ym19_1_1487706532585_14847"><br></div><div id="yui_3_16_0_ym19_1_1487706532585_14901">Com a ajuda do colega Michel Montoril, fizemos da seguinte forma:</div><div id="yui_3_16_0_ym19_1_1487706532585_14902">Seja X1,...,XR a.a de F (suporte real)<br></div><div id="yui_3_16_0_ym19_1_1487706532585_14957" dir="ltr">A: suporte do truncamento, P(A)=P(X pertencer a A).</div><div dir="ltr" id="yui_3_16_0_ym19_1_1487706532585_15124">E=E(h(X)|X pertence a A)</div><div id="yui_3_16_0_ym19_1_1487706532585_13121">Seja Yi=h(Xi)*Ind(Xi pertence a A)/P(A), i=1,...,R (nº de réplicas de Monte Carlo. Na verdade, gerei round(R/P(A)))<br></div><div id="yui_3_16_0_ym19_1_1487706532585_15171">E=mean(Y)</div><div id="yui_3_16_0_ym19_1_1487706532585_15343"><br></div><div id="yui_3_16_0_ym19_1_1487706532585_15344">Na verdade, fiz Yi=h(Xi)*Ind(Xi pertence a A) e depois E=mean(Y)/P(A).</div><div id="yui_3_16_0_ym19_1_1487706532585_15530"><br></div><div id="yui_3_16_0_ym19_1_1487706532585_15543">Usei o integrate do R em um exemplo e pareceu boa a aproximação (R=50000). Mas tenho vários cálculos, por isso, gero uma vez a amostra X e faço vários MC com a mesma amostra. Com integrate, teria de ser um para cada integral.<br></div><div id="yui_3_16_0_ym19_1_1487706532585_15531"><br></div><div id="yui_3_16_0_ym19_1_1487706532585_15693">Obrigado a todos pelas dicas.</div><div id="yui_3_16_0_ym19_1_1487706532585_15755">Clécio<br></div><div id="yui_3_16_0_ym19_1_1487706532585_13102"><span></span></div> <div class="qtdSeparateBR"><br><br></div><div class="yahoo_quoted" style="display: block;"> <div style="font-family: HelveticaNeue-Light, Helvetica Neue Light, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div dir="ltr"><font size="2" face="Arial"> Em Segunda-feira, 20 de Fevereiro de 2017 11:05, Clécio da Silva Ferreira por (abe) <abe@lists.ime.usp.br> escreveu:<br></font></div> <br><br> <div class="y_msg_container"><div id="yiv3245922403"><div><div style="color:#000;background-color:#fff;font-family:HelveticaNeue-Light, Helvetica Neue Light, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:16px;"><div id="yiv3245922403yui_3_16_0_ym19_1_1487596860647_11453" dir="ltr">Prezados, alguém tem referência de metodologia de geração de amostras aleatórias em distribuições truncadas? Aceitação-rejeição seria viável? Eu sei como gerar da distribuição (suporte na reta). Preciso gerar de valores na parte positiva. Pacote 'truncdist' do R não tem a distribuição.<br></div><div id="yiv3245922403yui_3_16_0_ym19_1_1487596860647_11500"><br></div><div id="yiv3245922403yui_3_16_0_ym19_1_1487596860647_11527">Atte,</div><div id="yiv3245922403yui_3_16_0_ym19_1_1487596860647_11528">Clécio<br></div></div></div></div><br>_______________________________________________<br>abe mailing list<br><a ymailto="mailto:abe@lists.ime.usp.br" href="mailto:abe@lists.ime.usp.br">abe@lists.ime.usp.br</a><br><a href="https://lists.ime.usp.br/listinfo/abe" target="_blank">https://lists.ime.usp.br/listinfo/abe</a><br><br><br></div> </div> </div> </div></div></body></html>