R help es - Mar 2016 - RV: heplots

De: Dr. José A. Betancourt Bethencourt [mailto:jbetancourt en iscmc.cmw.sld.cu]

Enviado el: lunes, 7 de marzo de 2016 05:40
Para: 'r-help-es en r-project.org' <r-help-es en r-project.org>
Asunto: heplots

 

Estimados, quisiéramos que el script que utilizamos brindara información de
la salida summary (soils.mod)  tales como R2=0.8543;p-value: < 2.2e-16 IMC
*** grasas_C *** aguas_c  ***

en el propio gráfico que se genera con el comando heplot(soils.mod) parecido
a lo del SjPlot

saludos

José

 

 

#SCRIPT

 

rm(list = ls())

#setwd("D:/Public/Documents/R/r_epidemiología/") 

#data<-#read.csv(paste("D:/Public/Documents/R/r_epidemiología/DATA/duda.csv"
),header=T#RUE, sep=";", dec=",")

attach(data)

soils.mod <- lm(cbind(Z50,  R50) ~IMC+grasas_C+aguas_c, data=data)

summary(soils.mod)

library(heplots)

heplot(soils.mod)

 



--
Este mensaje le ha llegado mediante el servicio de correo electronico que ofrece
Infomed para respaldar el cumplimiento de las misiones del Sistema Nacional de
Salud. La persona que envia este correo asume el compromiso de usar el servicio
a tales fines y cumplir con las regulaciones establecidas

Infomed: http://www.sld.cu/



------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20160307/c2cf3475/attachment.html>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: duda.csv
Type: application/vnd.ms-excel
Size: 2117 bytes
Desc: no disponible
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20160307/c2cf3475/attachment.xlb>

Hola,

Te adjunto como incluir una parte de la que pides.
Las que faltan, puedes seguir el ejemplo ...

#---------------------------
datIn <- read.table("duda.csv", dec=",", header=T,
sep=";")
attach(datIn)
soils.mod <- lm(cbind(Z50,  R50) ~ IMC + grasas_C + aguas_c, data=datIn)
summary(soils.mod)

*r2_z50 <- round(summary(soils.mod)$`Response Z50`$r.squared, 2)*
*r2_r50 <- round(summary(soils.mod)$`Response R50`$r.squared, 2)*


library(heplots)
heplot(soils.mod)
*text(700,400, bquote(R^2 ~~ .(paste("Z50: ",r2_z50))), cex=0.7,
col="tomato", font=2)*
*text(700,400*0.95, bquote(R^2 ~~ .(paste("R50: ",r2_r50))), cex=0.7,
col="tomato", font=2)*

#---------------------------

Saludos,
Carlos Ortega
www.qualityexcellence.es

El 7 de marzo de 2016, 11:43, Dr. José A. Betancourt Bethencourt <
josebetancourt.cmw en infomed.sld.cu> escribió:
>
>
>
>
> *De:* Dr. José A. Betancourt Bethencourt [mailto:
> jbetancourt en iscmc.cmw.sld.cu]
> *Enviado el:* lunes, 7 de marzo de 2016 05:40
> *Para:* 'r-help-es en r-project.org' <r-help-es en
r-project.org>
> *Asunto:* heplots
>
>
>
> Estimados, quisiéramos que el script que utilizamos brindara información
> de la salida summary (soils.mod)  tales como R2=0.8543;p-value: <
2.2e-16
> IMC  *** grasas_C *** aguas_c  ***
>
> en el propio gráfico que se genera con el comando heplot(soils.mod)
> parecido a lo del SjPlot
>
> saludos
>
> José
>
>
>
>
>
> #SCRIPT
>
>
>
> rm(list = ls())
>
> #setwd("D:/Public/Documents/R/r_epidemiología/")
>
>
#data<-#read.csv(paste("D:/Public/Documents/R/r_epidemiología/DATA/duda.csv"),header=T#RUE,
> sep=";", dec=",")
>
> attach(data)
>
> soils.mod <- lm(cbind(Z50,  R50) ~IMC+grasas_C+aguas_c, data=data)
>
> summary(soils.mod)
>
> library(heplots)
>
> heplot(soils.mod)
>
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>


-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]