<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:12pt"><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><span>Bom dia senhores!!</span></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><span><br></span></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><span>Estou encontrando um problema no ajuste de uma regressão segmentada, e não estou encontrando onde está este problema. Já fiz vários ajustes deste tipo, mais com este modelo particularmente tem alguma coisa errada que não sei onde é. No estudo dos pontos de chutes, o erro é claro, mais não estou sabendo como contornar. Deveria haver três segmentos, mais pelo estudo com janelas interativas, parece haver 5 segmentos. Se alguém puder ajudar, agradeço desde já. Segue um CMR.</span></div><div style="font-family: arial, helvetica, sans-serif;
font-size: 12pt; "><span><br></span></div><div><span><div>x <- c(0,0,0,10,10,10,20,20,20,35,35,35)</div><div>y <- c(7.77, 7.81, 7.81, 7.60, 7.67, 7.65, 7.60, NA, 7.66, 7.60, 7.42, 7.50)</div><div><br></div><div>segment_lpl <- function(x,b0,b1,b2,K,X0,X1){</div><div> (b0+b1*x)*(x<=X0)+</div><div> (K)*(X0<x & x<X1)+</div><div> (b0+b1*X0-b2*(X1-x))*(x>=X1)</div><div>}</div><div><br></div><div>teste <- nls(y ~ segment_lpl(x,b0,b1,b2,K,X0,X1),</div><div> start=list(b0=7.8,b1=-0.02,b2=-0.01,K=7.5,X0=10,X1=20)) </div><div><br></div><div>library(gWidgetsRGtk2) </div><div>limits <- list(b0=c(7,8), </div><div> b1=c(-0.05,0), </div><div> b2=c(-0.04,0), </div><div>
K=c(7,8), </div><div> X0=c(8,10),</div><div> X1=c(18,25))</div><div><br></div><div>start <- list() # lista com os valores para chute</div><div><br></div><div>#-----------------------------------------------------------------------------</div><div>#função que será atualizada a cada movimento do deslizador</div><div>#parâmetros dentro de svalue() são controlados, nomes igual aos da lista</div><div><br></div><div>plot.chute <- function(...){</div><div> #faz o gráfico de dispersão</div><div> plot(y ~ x)</div><div> #sobrepõe a curva com os valores dos deslizadores</div><div> curve(segment_lpl(x, svalue(b0), svalue(b1), svalue(b2), svalue(K),svalue(X0), svalue(X1)),</div><div>
add=TRUE, col=2)</div><div> #reescreve o start com os valores dos delizadores, para usar na nls()</div><div> start <- list(b0=svalue(b0), b1=svalue(b1), b2=svalue(b2), K=svalue(K), X0=svalue(X0), X1=svalue(X1))</div><div>}</div><div><br></div><div>#-----------------------------------------------------------------------------</div><div>#criação da janela com deslizadores</div><div>#na primeira chamada escolher uma das opções (sempre escolho a 1)</div><div>#Select a GUI toolkit</div><div>#essa função pode estar num arquivo fn.R e carregada com source("fn.R")</div><div><br></div><div>w <- gwindow("Caixa com deslizadores para controlar parâmetros")</div><div>tbl <- glayout(cont=w)</div><div>for(i in 1:length(limits)){</div><div> tbl[i,1] <- paste("Controle", names(limits)[i])</div><div> tbl[i,2, expand=TRUE] <- (assign(names(limits)[i],</div><div>
gslider(from=limits[[i]][1],</div><div> to=limits[[i]][2],</div><div> by=diff(limits[[i]])/20,</div><div> value=mean(limits[[i]]),</div><div> container=tbl,
handler=plot.chute)))</div><div>}</div><div># </div><div>#-----------------------------------------------------------------------------</div><div>#agora com a caixa criada, basta chamar a função e mover os deslizadores</div><div>plot.chute() </div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><br></div></span></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; ">(S,P,f)</div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; ">Allaman</div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; background-color: transparent; " align="center"><font style="background-color:transparent;" face="comic sans ms" size="2"><font size="1"><b><br></b></font></font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left;
background-color: transparent; " align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif" size="2">\begin{signature}</font></font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif" size="2"><<>>=</font></font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Prof. Dr. Ivan Bezerra Allaman</font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif"
size="2">Universidade Estadual de Santa Cruz</font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Departamento de Ciências Exatas e Tecnológicas</font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Ilhéus/BA - Brasil</font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Fone: +55 73 3680-5076</font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color:
transparent; " align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com</font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif" size="2">@</font></font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif" size="2">\end{signature}</font></font></div></div></body></html>