[R-br] Elipsóides de Confiança

[Fabio Mathias Corrêa]

----- Mensagem encaminhada ----
De: walmes zeviani <walmeszeviani em yahoo.com.br>
Para: R_STAT em yahoogrupos.com.br
Enviadas: Segunda-feira, 20 de Outubro de 2008 19:08:34
Assunto: Res: [R_STAT] Elipsóides de Confiança


Saudações Willian,

Sobre o que você relatou eu consegui a abordagem que está apresentada na 
seguência de comandos a seguir.

#====== Contruir contornos de confiança ======#
#======= para dois modelos de regressão ======#

library(ellipse) #para lidar com contornos
data(mtcars)     #conjunto de dados

fit <- lm(mpg ~ am + vs , mtcars)    #modelo 1
fit2 <- lm(mpg ~ drat + vs , mtcars) #modelo 2

# Plotando os contornos para o 'modelo 1'
 plot(ellipse( fit, which = c('am',  'vs'),
      level = 0.90), type = 'l',
      xlim=c(2,9), ylim=c(2,11) ,
      xlab="drat e am", ylab="vs",
      col=2, lwd=2,
      main="Contornos de confiança")
 par(new=T) #permite a sobreposição de gráficos
# Sobrepondo os contornos do 'modelo 2'
 plot(ellipse( fit2, which = c('drat', 'vs'),
      level = 0.90), type = 'l',
      xlim=c(2,9), ylim=c(2,11) ,
      xlab="", ylab="", col=3,
      lwd=2, main="")
# Adicionando informações úteis (acessórios)
 points(fit$coeffici ents['am' ],
        fit$coefficients[ 'vs'])
 points(fit2$ coefficients[ 'drat'],
         fit2$coefficients[ 'vs'])
 text(7,10,"mpg~ am+vs")
 text(3.5,3," mpg~drat+ vs")

#=========== ========= ========= ========= =======#


Tomando o cuidado de delarar o mesmo xlim ylim para os dois gráficos é o que 
importa. Procure delarar xlab ylab apenas no primeiro gráfico. Espero ter 
ajudado  você.

Uma ótima semana.
Walmes.



----- Mensagem original ----
De: William Marciano <fwpm2004 em yahoo. com.br>
Para: R_STAT em yahoogrupos. com.br
Enviadas: Segunda-feira, 20 de Outubro de 2008 16:54:34
Assunto: [R_STAT] Elipsóides de Confiança


Por favor, gostaria de saber se alguém já construiu em um mesmo gráfico 
elipsóides de confiança para  modelos de regressão padrão diferentes, já tentei 
com a função confidence.ellipse( ), mas não consegui colocar em um mesmo gráfico 
dois modelos que eu ajustei para compará-los. 

Desde já agradeço a atenção de todos!
 
____________ _________ _________ _________ _________ __
Fale com seus amigos  de graça com o novo Yahoo! Messenger 
http://br.messenger .yahoo.com/ 
____________ _________ _________ _________ _________ __
Fale com seus amigos  de graça com o novo Yahoo! Messenger 
http://br.messenger .yahoo.com/ __._,_.___ 
Mensagens neste tópico  (1)  Responder  (através da web)  | Adicionar um novo 
tópico  

Mensagens | Arquivos | Fotos | Links | Banco de dados | Enquetes 
 
Alterar configurações via web (Requer Yahoo! ID) 
Alterar configurações via e-mail: Alterar recebimento para lista diária de 
mensagens | Alterar formato para o tradicional 

Visite seu Grupo  | Termos de uso do Yahoo! Grupos  | Sair do grupo  
Atividade nos últimos dias
	*  8
Novos usuáriosVisite seu Grupo  
Yahoo! Mail
Conecte-se ao mundo
Proteção anti-spam
Muito mais espaço
Yahoo! Barra
Instale grátis
Buscar sites na web
Checar seus e-mails .
Yahoo! Grupos
Crie seu próprio grupo
A melhor forma de comunicação
. 
__,_._,___    


      
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110317/18f3d43a/attachment-0001.html>