Scatter plot de matriz
Prezados, Estou tentando inserir no "scatter plot de matriz" (veja o conjunto dados de exemplo, abaixo) as linhas de correlação entre as variáveis, entretanto não estou conseguindo de um modo mais direto. Seria para como usar para cada um dos gráficos o seguinte exemplo: #------- z <- lm(dist ~ speed, data = cars) plot(cars) abline(z) # equivalent to abline(reg = z) or abline(coef = coef(z)) #--------- Entretanto eu preciso colocar uma reta para cada plot abaixo #DADOS DE EXEMPLO var_1<-c(0.492077114759013 , 0.836016362765804 , 0.481298430124298 , 0.386469675228 , 0.777098514838144 , 0.617958343820646 , 0.219213913194835 , 0.678557483945042 , 0.423462828854099 , 0.829378235153854 , 0.806068887002766 , 0.337134317029268 , 0.661742311669514 , 0.619118098868057 , 0.788589890813455 , 0.279048688709736 , 0.820936637930572 , 0.247459845617414 , 0.752754796063527 , 0.934190914966166) var_2<-c(0.658414233683847 , 0.605210445823305 , 0.801341124716519 , 1.011018222403 , 0.572025150910094 , 0.836177895660953 , 1.03193629108894 , 0.656809315729027 , 0.949461899587905 , 0.553121786467357 , 0.817865236941146 , 0.925669452041572 , 0.578475899228501 , 0.476615221718288 , 0.571342170684858 , 0.867686931852481 , 0.615533868096807 , 0.926396173911956 , 0.667834629397337 , 0.65446214987638) var_3<-c(0.805190744576976 , 0.00357919884845614 , 0.296429808018729 , 0.147386643337086 , 0.224931408418342 , 0.967352513456717 , 0.88843358354643 , 0.991240894189104 , 0.906759183388203 , 0.635862796101719 , 0.90428285067901 , 0.755293083842844 , 0.484490071190521 , 0.0830997943412513 , 0.66727612237446 , 0.314308791654184 , 0.478725346503779 , 0.885625332361087 , 0.567443900275975 , 0.175529750995338) var_4<-c(1.51875793231285 , 2.08332774568708 , 1.85743044603302 , 1.94252530517918 , 1.70102115019736 , 1.4099218798547 , 1.37410255920602 , 1.63332942855045 , 1.59918963127413 , 1.54217901497363 , 1.39090041665641 , 1.56455567202422 , 1.70518441073333 , 1.96772731408659 , 1.76799247833419 , 1.85562083102216 , 1.81231983166867 , 1.58934006350458 , 1.77421167457209 , 2.01928129901513) data<-data.frame(var_1,var_2,var_3,var_4) plot(data,pch=16) Certo de vossa atenção, Atenciosamente, Luciano -- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732
Olha se eu entendi, vc quer colocar uma linha em cada plot. Veja a documentação de ?pairs e os exemplos la que podem te dar varias ideias. Mas acho que é isso que você precisa #faça uma função panel, que vai lidar com os pontos de cada quadradinho #aqui eu ploto os pontos e a linha. panel.lm <- function(x, y) { points(y~x) abline(lm(y~x)) } #dai é so entrar com a função, no caso eu escolhi faze isso na diagonal inferior, mas voce pode setar como quiser. pairs(data,lower.panel=panel.lm) Lembrando que plot e pairs são a mesma coisa nesse caso, com vc entrando com um data.frame. E evite postar em multiplas listas. Espero ter ajudado. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.wordpress.com/ Lattes: http://lattes.cnpq.br/7355685961127056
Augusto, Muito obrigado... era exatamente isso que eu estava procurando. Eu não vejo problema em fazer isso, visto que eu não fiz cross-post. O motivo de eu ter postado nas duas listas é que o publico compartilhado pelas duas é extremamente pequeno (são menos de 10%). Além disso, as mensagens podem ser úteis para outros usuários que não participam das duas. Abraços Luciano Em 13 de dezembro de 2012 17:56, Augusto Ribas <ribas.aca@gmail.com>escreveu:
Olha se eu entendi, vc quer colocar uma linha em cada plot.
Veja a documentação de ?pairs e os exemplos la que podem te dar varias ideias. Mas acho que é isso que você precisa
#faça uma função panel, que vai lidar com os pontos de cada quadradinho #aqui eu ploto os pontos e a linha.
panel.lm <- function(x, y) { points(y~x) abline(lm(y~x)) }
#dai é so entrar com a função, no caso eu escolhi faze isso na diagonal inferior, mas voce pode setar como quiser. pairs(data,lower.panel=panel.lm)
Lembrando que plot e pairs são a mesma coisa nesse caso, com vc entrando com um data.frame.
E evite postar em multiplas listas.
Espero ter ajudado. -- Grato Augusto C. A. Ribas
Site Pessoal: http://recologia.wordpress.com/ Lattes: http://lattes.cnpq.br/7355685961127056 _______________________________________________ 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.
-- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732
participantes (2)
-
Augusto Ribas -
Luciano F. Sgarbi