<div dir="ltr">Olá, <div>você diz eu colar aqui o dput ou mandar em anexo o data-set?</div></div><div class="gmail_extra"><br><div class="gmail_quote">Em 8 de outubro de 2015 23:04, sznelwar <span dir="ltr"><<a href="mailto:sznelwar@uol.com.br" target="_blank">sznelwar@uol.com.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><span style="font-size:medium;color:#000080">Tem o data-set para rodar?</span><br> <br><br></div>
<hr style="border-top:1px solid #ccc">
<div><br><br></div>
<div dir="ltr"><div><div class="h5">Caros,boa noite,Â
<div>Â </div>
<div>estou tentando fazer um fit em um conjunto de dados mas ele me retorna o erro:
<div>
<div>Â </div>
<div>Erro em nls(y ~ func(x, c, a, b), data = dados2, start = guess, trace = TRUE) :Â </div>
<div>Â fator de passos 0.000488281 reduzido abaixo de 'minFactor' de 0.000976562</div>
</div>
</div>
<div>Â </div>
<div>Alguém, por favor, saberia me dizer como posso concertar isso?</div>
<div>O programa segue abaixo:</div>
</div></div><div><div><div class="h5">
<div>> require(lattice)</div>
<div>> l<-1</div>
<div>> file<-read.table(paste0("/Users/bau/ProjetoU-Shape/AA.N_Novo_MomentosEstatÃsticosSemestre",l,".txt"), header=TRUE)</div>
<div>> media<-file[3:389,2]</div>
<div>> media</div>
<div>Â [1] 24117.500 13168.333 13861.667 18031.667 19581.667 18797.500 15147.500 16699.167 17911.667</div>
<div>Â [10] 16430.000 15874.167 15470.000 16360.833 18240.000 15018.333 18155.833 14793.333 17348.333</div>
<div>Â [19] 19142.500 21510.833 17417.500 17450.833 19335.833 14439.167 16995.000 14624.167 18152.500</div>
<div>Â [28] 18240.000 16124.167 18489.167 14679.167 17181.667 18115.833 15302.500 13445.833 15159.167</div>
<div>Â [37] 20161.667 15321.667 14252.500 15850.833 13745.000 13915.833 15810.833 13280.833 14680.833</div>
<div>Â [46] 12119.167 12100.833 16263.333 13210.833 15323.333 12921.667 14374.167 12805.833 13955.000</div>
<div>Â [55] 12920.833 12016.667 13795.000 14751.667 12735.833 12267.500 12896.667 11947.500 12723.333</div>
<div>Â [64] 14574.167 12818.333 14187.500 11534.167 12134.167 10910.000 12596.667 11503.333 12515.833</div>
<div>Â [73] 12709.167 14239.167 11495.000 13163.333 11990.000 10506.667 11984.167 10225.000 11447.500</div>
<div>Â [82] 11966.667 12115.833 11273.333 10478.333 10375.000 12240.833 11178.333 10699.167 11474.167</div>
<div>Â [91] 10075.000 12563.333 14194.167 13675.833 12290.000 11357.500 10956.667 12120.000 11663.333</div>
<div>[100] 11078.333 11371.667 10549.167 12766.667 14523.333 10948.333 10809.167 10772.500 11190.000</div>
<div>[109] Â 9748.333 14768.333 11070.000 10362.500 10942.500 Â 9190.833 12433.333 10302.500 10838.333</div>
<div>[118] 12184.167 Â 8899.167 10073.333 10393.333 Â 9760.000 10390.833 11309.167 Â 9319.167 Â 8183.333</div>
<div>[127] Â 9531.667 Â 9843.333 Â 8924.167 10070.833 Â 8730.000 10109.167 10507.500 11991.667 Â 8540.833</div>
<div>[136] 11396.667 Â 9405.000 Â 9284.167 Â 9146.667 Â 9640.000 Â 8991.667 Â 8352.500 Â 7969.167 10366.667</div>
<div>[145] Â 9319.167 Â 8555.000 Â 8578.333 Â 7897.500 Â 8881.667 Â 8842.500 Â 9579.167 Â 9022.500 Â 9914.167</div>
<div>[154] Â 9933.333 10987.500 Â 8835.000 Â 8473.333 Â 9029.167 Â 8619.167 Â 8782.500 Â 9901.667 Â 8566.667</div>
<div>[163] Â 8490.833 Â 7475.833 Â 7510.833 Â 8990.000 Â 6865.833 Â 6564.167 Â 7418.333 Â 7356.667 Â 7485.000</div>
<div>[172] Â 8482.500 Â 9325.000 Â 7989.167 Â 7370.833 Â 7921.667 Â 8243.333 Â 8005.833 Â 8551.667 Â 9119.167</div>
<div>[181] Â 8760.000 Â 7193.333 Â 9693.333 Â 7742.500 11112.500 Â 7497.500 Â 6685.000 Â 8441.667 Â 7150.000</div>
<div>[190] Â 6189.167 Â 6760.833 Â 6255.833 Â 8953.333 Â 8413.333 Â 7731.667 Â 6546.667 Â 7179.167 Â 8060.000</div>
<div>[199] Â 7023.333 Â 9560.000 Â 7071.667 Â 7770.000 Â 6640.000 Â 5910.833 Â 8896.667 Â 7938.333 Â 7098.333</div>
<div>[208] Â 7282.500 Â 7291.667 Â 8557.500 Â 7326.667 Â 8196.667 12226.667 Â 8380.000 Â 6513.333 Â 9100.833</div>
<div>[217] Â 6303.333 Â 7525.000 Â 7297.500 Â 7375.833 Â 7224.167 Â 7029.167 Â 9325.833 Â 5923.333 Â 7674.167</div>
<div>[226] Â 5882.500 Â 8318.333 Â 7263.333 Â 7574.167 Â 7268.333 Â 7530.000 Â 9355.833 Â 8654.167 Â 7426.667</div>
<div>[235] Â 7794.167 Â 7008.333 Â 7974.167 10410.000 Â 7181.667 Â 7440.000 Â 6345.000 Â 5503.333 Â 6674.167</div>
<div>[244] Â 9780.000 Â 7740.000 Â 6981.667 Â 5926.667 Â 9352.500 10264.167 Â 8406.667 Â 6532.500 Â 6478.333</div>
<div>[253] Â 8561.667 Â 8072.500 Â 7700.000 Â 7189.167 Â 6443.333 Â 8120.000 Â 7190.833 Â 6520.833 Â 7539.167</div>
<div>[262] Â 8013.333 Â 7191.667 Â 7660.000 Â 7276.667 Â 8233.333 Â 7500.833 Â 8582.500 Â 8019.167 Â 7231.667</div>
<div>[271] 10179.167 Â 8445.000 Â 9302.500 Â 7680.833 Â 9114.167 Â 8432.500 Â 7475.000 Â 7779.167 Â 8895.000</div>
<div>[280] Â 9171.667 Â 9753.333 Â 7490.000 10211.667 11380.000 Â 8524.167 Â 8077.500 10155.000 Â 9406.667</div>
<div>[289] 11199.167 Â 8286.667 Â 9850.000 10032.500 Â 9740.833 Â 7321.667 Â 8494.167 10023.333 Â 9450.000</div>
<div>[298] 10995.833 10166.667 10881.667 Â 8832.500 10015.000 11526.667 11015.833 10857.500 Â 9627.500</div>
<div>[307] 11716.667 10213.333 Â 9765.000 Â 8673.333 Â 8560.833 10153.333 12676.667 11633.333 10933.333</div>
<div>[316] 10492.500 10300.000 Â 9853.333 10602.500 Â 9778.333 Â 9030.000 12785.833 10950.000 11523.333</div>
<div>[325] 12445.000 10896.667 10875.833 11619.167 13154.167 11693.333 12561.667 11002.500 11017.500</div>
<div>[334] 10700.833 14557.500 12236.667 10875.833 11762.500 12705.000 13452.500 11402.500 11232.500</div>
<div>[343] 11739.167 13043.333 11583.333 11468.333 11772.500 13278.333 14787.500 13462.500 13837.500</div>
<div>[352] 13491.667 12455.000 14412.500 14211.667 15135.000 13923.333 15175.000 16048.333 16182.500</div>
<div>[361] 17939.167 17191.667 18522.500 18923.333 17106.667 15255.833 14752.500 15868.333 21636.667</div>
<div>[370] 17807.500 19451.667 17064.167 20163.333 16959.167 17077.500 19470.833 17424.167 18996.667</div>
<div>[379] 24262.500 21454.167 21115.833 18455.000 22215.000 27894.167 28049.167 30302.500 37730.000</div>
<div>> y<-media</div>
<div>> tempo<-c(-193:193/193)</div>
<div>> tempo</div>
<div>Â [1] -1.000000000 -0.994818653 -0.989637306 -0.984455959 -0.979274611 -0.974093264 -0.968911917</div>
<div>Â [8] -0.963730570 -0.958549223 -0.953367876 -0.948186528 -0.943005181 -0.937823834 -0.932642487</div>
<div>Â [15] -0.927461140 -0.922279793 -0.917098446 -0.911917098 -0.906735751 -0.901554404 -0.896373057</div>
<div>Â [22] -0.891191710 -0.886010363 -0.880829016 -0.875647668 -0.870466321 -0.865284974 -0.860103627</div>
<div>Â [29] -0.854922280 -0.849740933 -0.844559585 -0.839378238 -0.834196891 -0.829015544 -0.823834197</div>
<div>Â [36] -0.818652850 -0.813471503 -0.808290155 -0.803108808 -0.797927461 -0.792746114 -0.787564767</div>
<div>Â [43] -0.782383420 -0.777202073 -0.772020725 -0.766839378 -0.761658031 -0.756476684 -0.751295337</div>
<div>Â [50] -0.746113990 -0.740932642 -0.735751295 -0.730569948 -0.725388601 -0.720207254 -0.715025907</div>
<div>Â [57] -0.709844560 -0.704663212 -0.699481865 -0.694300518 -0.689119171 -0.683937824 -0.678756477</div>
<div>Â [64] -0.673575130 -0.668393782 -0.663212435 -0.658031088 -0.652849741 -0.647668394 -0.642487047</div>
<div>Â [71] -0.637305699 -0.632124352 -0.626943005 -0.621761658 -0.616580311 -0.611398964 -0.606217617</div>
<div>Â [78] -0.601036269 -0.595854922 -0.590673575 -0.585492228 -0.580310881 -0.575129534 -0.569948187</div>
<div>Â [85] -0.564766839 -0.559585492 -0.554404145 -0.549222798 -0.544041451 -0.538860104 -0.533678756</div>
<div>Â [92] -0.528497409 -0.523316062 -0.518134715 -0.512953368 -0.507772021 -0.502590674 -0.497409326</div>
<div>Â [99] -0.492227979 -0.487046632 -0.481865285 -0.476683938 -0.471502591 -0.466321244 -0.461139896</div>
<div>[106] -0.455958549 -0.450777202 -0.445595855 -0.440414508 -0.435233161 -0.430051813 -0.424870466</div>
<div>[113] -0.419689119 -0.414507772 -0.409326425 -0.404145078 -0.398963731 -0.393782383 -0.388601036</div>
<div>[120] -0.383419689 -0.378238342 -0.373056995 -0.367875648 -0.362694301 -0.357512953 -0.352331606</div>
<div>[127] -0.347150259 -0.341968912 -0.336787565 -0.331606218 -0.326424870 -0.321243523 -0.316062176</div>
<div>[134] -0.310880829 -0.305699482 -0.300518135 -0.295336788 -0.290155440 -0.284974093 -0.279792746</div>
<div>[141] -0.274611399 -0.269430052 -0.264248705 -0.259067358 -0.253886010 -0.248704663 -0.243523316</div>
<div>[148] -0.238341969 -0.233160622 -0.227979275 -0.222797927 -0.217616580 -0.212435233 -0.207253886</div>
<div>[155] -0.202072539 -0.196891192 -0.191709845 -0.186528497 -0.181347150 -0.176165803 -0.170984456</div>
<div>[162] -0.165803109 -0.160621762 -0.155440415 -0.150259067 -0.145077720 -0.139896373 -0.134715026</div>
<div>[169] -0.129533679 -0.124352332 -0.119170984 -0.113989637 -0.108808290 -0.103626943 -0.098445596</div>
<div>[176] -0.093264249 -0.088082902 -0.082901554 -0.077720207 -0.072538860 -0.067357513 -0.062176166</div>
<div>[183] -0.056994819 -0.051813472 -0.046632124 -0.041450777 -0.036269430 -0.031088083 -0.025906736</div>
<div>[190] -0.020725389 -0.015544041 -0.010362694 -0.005181347 Â 0.000000000 Â 0.005181347 Â 0.010362694</div>
<div>[197] Â 0.015544041 Â 0.020725389 Â 0.025906736 Â 0.031088083 Â 0.036269430 Â 0.041450777 Â 0.046632124</div>
<div>[204] Â 0.051813472 Â 0.056994819 Â 0.062176166 Â 0.067357513 Â 0.072538860 Â 0.077720207 Â 0.082901554</div>
<div>[211] Â 0.088082902 Â 0.093264249 Â 0.098445596 Â 0.103626943 Â 0.108808290 Â 0.113989637 Â 0.119170984</div>
<div>[218] Â 0.124352332 Â 0.129533679 Â 0.134715026 Â 0.139896373 Â 0.145077720 Â 0.150259067 Â 0.155440415</div>
<div>[225] Â 0.160621762 Â 0.165803109 Â 0.170984456 Â 0.176165803 Â 0.181347150 Â 0.186528497 Â 0.191709845</div>
<div>[232] Â 0.196891192 Â 0.202072539 Â 0.207253886 Â 0.212435233 Â 0.217616580 Â 0.222797927 Â 0.227979275</div>
<div>[239] Â 0.233160622 Â 0.238341969 Â 0.243523316 Â 0.248704663 Â 0.253886010 Â 0.259067358 Â 0.264248705</div>
<div>[246] Â 0.269430052 Â 0.274611399 Â 0.279792746 Â 0.284974093 Â 0.290155440 Â 0.295336788 Â 0.300518135</div>
<div>[253] Â 0.305699482 Â 0.310880829 Â 0.316062176 Â 0.321243523 Â 0.326424870 Â 0.331606218 Â 0.336787565</div>
<div>[260] Â 0.341968912 Â 0.347150259 Â 0.352331606 Â 0.357512953 Â 0.362694301 Â 0.367875648 Â 0.373056995</div>
<div>[267] Â 0.378238342 Â 0.383419689 Â 0.388601036 Â 0.393782383 Â 0.398963731 Â 0.404145078 Â 0.409326425</div>
<div>[274] Â 0.414507772 Â 0.419689119 Â 0.424870466 Â 0.430051813 Â 0.435233161 Â 0.440414508 Â 0.445595855</div>
<div>[281] Â 0.450777202 Â 0.455958549 Â 0.461139896 Â 0.466321244 Â 0.471502591 Â 0.476683938 Â 0.481865285</div>
<div>[288] Â 0.487046632 Â 0.492227979 Â 0.497409326 Â 0.502590674 Â 0.507772021 Â 0.512953368 Â 0.518134715</div>
<div>[295] Â 0.523316062 Â 0.528497409 Â 0.533678756 Â 0.538860104 Â 0.544041451 Â 0.549222798 Â 0.554404145</div>
<div>[302] Â 0.559585492 Â 0.564766839 Â 0.569948187 Â 0.575129534 Â 0.580310881 Â 0.585492228 Â 0.590673575</div>
<div>[309] Â 0.595854922 Â 0.601036269 Â 0.606217617 Â 0.611398964 Â 0.616580311 Â 0.621761658 Â 0.626943005</div>
<div>[316] Â 0.632124352 Â 0.637305699 Â 0.642487047 Â 0.647668394 Â 0.652849741 Â 0.658031088 Â 0.663212435</div>
<div>[323] Â 0.668393782 Â 0.673575130 Â 0.678756477 Â 0.683937824 Â 0.689119171 Â 0.694300518 Â 0.699481865</div>
<div>[330] Â 0.704663212 Â 0.709844560 Â 0.715025907 Â 0.720207254 Â 0.725388601 Â 0.730569948 Â 0.735751295</div>
<div>[337] Â 0.740932642 Â 0.746113990 Â 0.751295337 Â 0.756476684 Â 0.761658031 Â 0.766839378 Â 0.772020725</div>
<div>[344] Â 0.777202073 Â 0.782383420 Â 0.787564767 Â 0.792746114 Â 0.797927461 Â 0.803108808 Â 0.808290155</div>
<div>[351] Â 0.813471503 Â 0.818652850 Â 0.823834197 Â 0.829015544 Â 0.834196891 Â 0.839378238 Â 0.844559585</div>
<div>[358] Â 0.849740933 Â 0.854922280 Â 0.860103627 Â 0.865284974 Â 0.870466321 Â 0.875647668 Â 0.880829016</div>
<div>[365] Â 0.886010363 Â 0.891191710 Â 0.896373057 Â 0.901554404 Â 0.906735751 Â 0.911917098 Â 0.917098446</div>
<div>[372] Â 0.922279793 Â 0.927461140 Â 0.932642487 Â 0.937823834 Â 0.943005181 Â 0.948186528 Â 0.953367876</div>
<div>[379] Â 0.958549223 Â 0.963730570 Â 0.968911917 Â 0.974093264 Â 0.979274611 Â 0.984455959 Â 0.989637306</div>
<div>[386] Â 0.994818653 Â 1.000000000</div>
<div>> x<-tempo</div>
<div>>Â </div>
</div></div><div>> dados
<div>> dados2
<span class=""><div>>Â </div>
<div>> xyplot(y~x, data=dados, type=c("p","smooth"))</div>
<div>>Â </div>
<div>> ## Valores iniciais.</div>
<div>>Â </div>
</span><div>> guess
<div>>Â </div>
<div>> ## Modelo não linear.</div>
<div>> func
<div>+ Â y
<span class=""><div>+ Â return(y)</div>
<div>+ }</div>
<div>> with(guess,</div>
<div>+ Â Â Â curve(func(x, c,a,b),min(dados$x), max(dados$x)))</div>
<div>Mensagens de aviso perdidas:</div>
<div>In atanh((sqrt((x + (b * (x^2)))^2)/c)^2) : NaNs produzidos</div>
<div>> abline(v=194, lty=2)</div>
<div>>Â </div>
<div>>Â </div>
<div>> ## Sobrepondo dados e função a partir dos valores iniciais.</div>
<div>> plot(y~x, data=dados, xlab="tempo", ylab="volume/média")</div>
<div>>Â </div>
<div>>Â </div>
</span><div>> f expression(atanh((sqrt((x+(b*(x^2)))^2)/c)^2) + a)</div>
<div>>Â </div>
<div>> fit
<span class=""><div>55389072390 : Â 1.000 1.000 0.001</div>
<div>6020970037 : Â -224.49929 10765.82603 Â Â 49.04659</div>
<div>Erro em nls(y ~ func(x, c, a, b), data = dados2, start = guess, trace = TRUE) :Â </div>
<div>Â fator de passos 0.000488281 reduzido abaixo de 'minFactor' de 0.000976562</div>
<div>>Â </div>
<div>Â </div>
<div>Muito Obrigada,</div>
</span></div>
<div>Â </div>
<div>Michelle</div>
<div>Â </div>
</div>
_______________________________________________<br>R-br mailing list<br><a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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 c�digo m�nimo reproduz�vel.</div>
</div>
</div>
</div>
</div>
</div><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" rel="noreferrer" 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" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mÃnimo reproduzÃvel.<br></blockquote></div><br></div>