Modelos segmentados de regressão não linear...

Ajuste o modelo com a nls() por exemplo e peça o confint() do modelo ajustado. Forneça um CMR. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

Bem, você não passou CMR portanto eu não sei como anda sua programação, mas o linear platô seria algo assim ym <- 0.5; # altura do platô xm <- 6; # ponto do cotôvelo ou troca de regime/taxa cm <- -0.1; # taxa de decaimento curve(ym+cm*(x-xm)^1*(x<=xm), 0, 10, main="linear platô") x <- seq(0,10,by=0.5) y <- ym+cm*(x-xm)^1*(x<=xm)+rnorm(x,0,0.01) plot(y~x) n0 <- nls(y~ym+cm*(x-xm)*(x<=xm), start=list(ym=0.5, cm=-0.1, xm=5)) summary(n0) confint(n0) confint.default(n0) À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

Olá Walmes, aproveitando a resposta que você deu ao André, tem como plotar o intervalo de confiança no gráfico. No meu caso estou com o modelo logístico y=A/(1-B*exp(-C*x)). Att. Tiago. Date: Wed, 24 Oct 2012 22:19:53 -0200 From: walmeszeviani@gmail.com To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Modelos segmentados de regressão não linear... Bem, você não passou CMR portanto eu não sei como anda sua programação, mas o linear platô seria algo assim ym <- 0.5; # altura do platôxm <- 6; # ponto do cotôvelo ou troca de regime/taxacm <- -0.1; # taxa de decaimento curve(ym+cm*(x-xm)^1*(x<=xm), 0, 10, main="linear platô") x <- seq(0,10,by=0.5) y <- ym+cm*(x-xm)^1*(x<=xm)+rnorm(x,0,0.01)plot(y~x) n0 <- nls(y~ym+cm*(x-xm)*(x<=xm), start=list(ym=0.5, cm=-0.1, xm=5))summary(n0)confint(n0) confint.default(n0) À disposição. Walmes.========================================================================== Walmes Marques ZevianiLEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paranáfone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173e-mail: walmes@ufpr.br skype: walmeszevianitwitter: @walmeszevianihomepage: http://www.leg.ufpr.br/~walmes linux user number: 531218========================================================================== _______________________________________________ R-br mailing list R-br@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.

Seu exemplo não é reproduzível, logo na primeira linha dá erro!
cv=c(x1,x2,x3,x4,x6,x8,x9,x12,x16,x18,x24,x32,x36,x48,x64,x72,x96,x144,x192,x288,x384,x576);cv Erro: objeto 'x1' não encontrado Para plotar o IC, é só fazer aux <- confint(modelo) plot(aux) Se quiser um gráfico com intervalos do tipo ponto médio + extremos, pode usar a função latticeExtra::segplot(), consulte a sessão exemplos da documentação da função para ver. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

Apesar de você me mandar os dados colados no e-mail isso ainda não é diretamente reproduzível. Vou ter que jogar num arquivo e importar e tals e isso toma 40x mais tempo do que você, principal envolvido/interessado, colar o resultado de um dput(), recorde a saída do comando dput() e cole na mensagem. A lógica é: facilite a vida de quem tá disposto a te ajudar. Leia o guia de postagem. dput(dados) À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

André, Sinto informar mas não facilitou ainda! Fácil é resolver com control+c e control+v. Adote um dos dois procedimentos: 1. textConnection(); 2. dput(); Veja como fazer aqui: http://www.leg.ufpr.br/doku.php/ridiculas?grafico_com_dois_eixos_coordenad... Um CMR significa Código Mínimo **Reproduzível**, e reproduzível é daquele tipo "Recorte o código da mensagem, cole numa sessão R e veja a mágica acontecer. Ou os problemas!". À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (3)
-
andrebvs
-
tiago souza marçal
-
Walmes Zeviani