[R-es] RV: gráfico dispersión

2017-07-07 Thread Dr . José A . Betancourt Bethencourt
 

 

De: Dr. Jos� A. Betancourt Bethencourt [mailto:jbetanco...@iscmc.cmw.sld.cu]

Enviado el: viernes, 7 de julio de 2017 12:05
Para: 'r-help-es@r-project.org' 
Asunto: RV: gr�fico dispersi�n

 

 

 

De: Dr. Jos� A. Betancourt Bethencourt [mailto:jbetanco...@iscmc.cmw.sld.cu]

Enviado el: viernes, 7 de julio de 2017 12:05
Para: 'r-help-es@r-project.org'  >
CC: 'josebetancourt@infomed.sld.cu'  >
Asunto: gr�fico dispersi�n

 

Estimados

 

He buscado en ggplot2 y no encuentro un script que me permita plotear una
variable cuantitativa continua (representada con una   l�nea) contra cuatro
variables cuantitativas ( representadas con puntos de diferentes colores)

 

Mis saludos

Jos�



--
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/




[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] RV: pregunta

2016-10-27 Thread Dr. José A. Betancourt Bethencourt
 

 

De: Dr. Jos� A. Betancourt Bethencourt [mailto:jbetanco...@iscmc.cmw.sld.cu]

Enviado el: jueves, 27 de octubre de 2016 05:56
Para: 'r-help-es@r-project.org' 
Asunto: pregunta

 

Estimados, el paquete EPICALC ya no funciona, �existe alguna manera de
utilizarlo sin tener que regresar a una vieja versi�n de r? Es que para la
docencia es muy �til

Saludos

Jos�



--
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/




[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] de pdf a csv

2016-09-10 Thread Dr. José A. Betancourt Bethencourt
Estimados

 

En ocasionas hay informaciones epidemiológicas en reportes pdf semanales
como el que adjunto que quisiéramos llevar a csv o txt  USANDO R para poder
analizarlas estadísticamente. Apreciaríamos su ayuda si nos diesen un
script, el paquete pdftable no me resultó.

Saludos

José



--
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/





08-2016.pdf
Description: Adobe PDF document
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] RV: zyka RD

2016-08-16 Thread Dr. José A. Betancourt Bethencourt
 

Estimados

 

Quisiera saber cómo combinar varios archivos csv  pero no de manera manual,
sino con un script

 

Saludos

José



--
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/





Epidemiological_Bulletin-2016-06-04.csv
Description: MS-Excel spreadsheet


Epidemiological_Bulletin-2016-01-09.csv
Description: MS-Excel spreadsheet


Epidemiological_Bulletin-2016-01-16.csv
Description: MS-Excel spreadsheet


Epidemiological_Bulletin-2016-01-23.csv
Description: MS-Excel spreadsheet


Epidemiological_Bulletin-2016-01-30.csv
Description: MS-Excel spreadsheet
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] RV: zyka RD

2016-08-16 Thread Dr. José A. Betancourt Bethencourt
 

 

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

Enviado el: martes, 16 de agosto de 2016 06:18
Para: 'r-help-es@r-project.org' <r-help-es@r-project.org>
Asunto: zyka RD

 

Estimados

 

Quisiera saber cómo combinar varios archivos csv  pero no de manera manual,
sino con un script

 

Saludos

José



--
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/





Epidemiological_Bulletin-2016-06-04.csv
Description: MS-Excel spreadsheet


Epidemiological_Bulletin-2016-01-09.csv
Description: MS-Excel spreadsheet


Epidemiological_Bulletin-2016-01-16.csv
Description: MS-Excel spreadsheet


Epidemiological_Bulletin-2016-01-23.csv
Description: MS-Excel spreadsheet


Epidemiological_Bulletin-2016-01-30.csv
Description: MS-Excel spreadsheet
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] pregunta de fechas 1

2016-07-17 Thread Dr. José A. Betancourt Bethencourt
Estimados

 

en algunas tablas meteorológicas aparece este formato 1-1-2015

 

las funciones que tengo solo leen el formato   2010-1-1

 

?Como puedo convertir el primer formato en el segundo? Adjunto las dos bases
de datos

 

Saludos

José

 



--
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/



CST,Max TemperatureF,Mean TemperatureF,Min TemperatureF,Max Dew PointF,MeanDew 
PointF,Min DewpointF,Max Humidity, Mean Humidity, Min Humidity, Max Sea Level 
PressureIn, Mean Sea Level PressureIn, Min Sea Level PressureIn, Max 
VisibilityMiles, Mean VisibilityMiles, Min VisibilityMiles, Max Wind SpeedMPH, 
Mean Wind SpeedMPH, Max Gust SpeedMPH,PrecipitationIn, CloudCover, Events, 
WindDirDegrees
1-1-2015,85,76,67,70,66,61,94,74,39,30.21,30.13,30.06,6,6,6,18,10,,0,2,,63
1-2-2015,86,77,68,70,68,64,94,77,45,30.22,30.14,30.09,6,6,3,20,10,,0,2,Fog,76
1-3-2015,86,77,68,70,67,64,88,74,44,30.22,30.16,30.12,6,6,5,22,11,,0,2,,81
1-4-2015,88,78,68,70,67,59,94,73,30,30.25,30.17,30.09,6,6,6,18,10,,0,2,,81
1-5-2015,86,78,69,70,67,64,94,74,43,30.26,30.2,30.15,6,6,4,21,11,,0,1,Fog,70
1-6-2015,85,77,69,70,67,63,94,74,42,30.26,30.18,30.12,6,6,6,20,11,29,0,3,,63
1-7-2015,86,76,66,68,64,52,94,70,22,30.18,30.11,30.04,6,6,4,16,8,,0,2,Fog,61
1-8-2015,83,74,65,68,64,60,94,73,36,30.18,30.13,30.09,6,6,6,20,10,,0,3,,38
1-9-2015,85,76,66,66,64,61,88,71,39,30.19,30.13,30.06,6,6,5,18,11,,0,1,,61
1-10-2015,84,76,67,66,65,63,88,71,40,30.23,30.17,30.12,6,6,6,22,11,,0,3,,55
1-11-2015,83,76,68,72,67,64,94,78,54,30.19,30.14,30.09,6,6,4,21,12,,0.28,3,Rain,64
1-12-2015,85,76,68,70,68,66,94,78,51,30.16,30.09,30.04,6,6,6,20,9,,0,4,,76
1-13-2015,88,78,67,70,68,64,94,76,40,30.12,30.05,30.01,6,6,4,9,6,,0,2,Fog,84
1-14-2015,87,76,65,70,67,63,94,77,39,30.09,30.02,29.98,6,4,0,12,4,,0,2,Fog,10
1-15-2015,87,76,64,70,65,61,94,75,36,30.12,30.04,29.98,6,5,2,14,4,,0,1,Fog,32
1-16-2015,84,73,62,68,63,59,94,73,32,30.09,30.04,30.01,6,4,0,13,4,,0,2,Fog,323
1-17-2015,82,73,64,68,66,63,94,82,48,30.13,30.07,30.04,6,5,4,15,5,,0.08,4,Fog-Rain,53
1-18-2015,86,76,66,68,66,61,94,74,40,30.14,30.08,30.01,6,6,6,14,7,,0,3,,70
1-19-2015,78,72,66,68,65,61,94,80,54,30.15,30.08,30.04,6,6,3,16,7,,0.2,4,Rain,36
1-20-2015,82,74,65,66,64,61,94,75,47,30.09,30.04,29.98,6,6,6,16,8,,0,3,,71
1-21-2015,87,76,64,68,64,59,94,74,34,30.09,30.04,29.98,6,5,1,14,6,,0,1,Fog,60
1-22-2015,85,76,67,70,65,60,94,73,35,30.16,30.1,30.04,6,6,4,17,7,,0,2,Fog,61
1-23-2015,87,76,65,66,64,59,94,71,31,30.12,30.06,29.98,6,6,5,16,7,,0,2,,76
1-24-2015,88,76,65,68,63,55,94,69,27,30.09,30.02,29.95,6,6,3,14,5,21,0,2,Fog,233
1-25-2015,84,74,65,68,65,62,94,77,42,30.08,30.02,29.95,6,5,0,13,4,,0,2,Fog,16
1-26-2015,82,73,64,66,62,57,89,71,44,30.06,30.01,29.95,6,6,5,16,4,,0,2,Fog,305
1-27-2015,82,70,59,59,55,48,88,63,27,30.1,30.04,29.98,6,6,5,14,5,26,0,1,,327
1-28-2015,79,66,54,57,54,51,88,65,28,30.15,30.09,30.04,6,6,6,13,5,,0,2,,14
1-29-2015,78,68,59,59,55,53,82,63,34,30.17,30.12,30.06,6,6,6,21,10,,0,2,,41
1-30-2015,80,71,62,59,56,52,83,61,31,30.17,30.11,30.06,6,6,6,15,9,,0,2,,38
1-31-2015,78,68,58,59,56,50,83,64,31,30.18,30.12,30.09,6,6,6,21,9,,0,2,,46
2-1-2015,81,72,62,64,62,61,94,72,41,30.16,30.1,30.06,6,6,6,21,10,,0,3,,72
2-2-2015,89,78,66,68,63,57,94,68,28,30.13,30.06,30.01,6,6,6,13,6,,0,2,,84
2-3-2015,87,76,65,66,64,61,94,70,32,30.15,30.1,30.06,6,6,4,20,6,,0,2,Fog,48
2-4-2015,86,76,66,68,65,62,89,72,34,30.13,30.05,29.98,6,6,3,18,10,,0,2,,70
2-5-2015,88,77,66,68,63,55,94,66,27,30.07,30,29.95,6,6,6,15,7,,0,1,,181
2-6-2015,87,76,64,68,64,59,94,71,32,30.09,30.03,29.98,6,6,4,18,4,,0,2,Fog,25
2-7-2015,77,72,68,64,63,61,83,72,49,30.14,30.08,30.04,6,6,6,22,13,,0,4,,55
2-8-2015,84,73,62,64,61,59,88,66,35,30.13,30.05,29.98,6,6,6,18,10,,0,3,,76
2-9-2015,84,76,67,68,65,61,88,75,35,30.03,29.93,29.86,6,5,1,12,5,,0.67,5,Rain,69
2-10-2015,81,73,64,68,60,55,94,64,32,29.9,29.84,29.77,6,6,6,21,10,,0,3,,282
2-11-2015,78,67,56,55,51,46,88,58,25,30.04,29.96,29.86,6,6,6,12,6,,0,1,,344
2-12-2015,80,67,54,57,52,48,77,59,28,30.1,30.04,29.98,6,6,6,14,4,,0,1,,17
2-13-2015,79,68,56,57,53,47,94,60,23,30.17,30.1,30.04,6,6,6,13,5,,0,2,,28
2-14-2015,81,70,58,59,54,50,88,62,26,30.18,30.13,30.09,6,6,6,14,6,,0,1,,360
2-15-2015,76,66,56,54,52,48,88,62,27,30.21,30.14,30.09,6,6,6,15,8,,0,1,,34
2-16-2015,82,70,57,59,54,50,83,61,24,30.17,30.09,30.01,6,6,6,18,9,,0,1,,76
2-17-2015,88,74,60,64,61,59,94,65,31,30.1,30.02,29.98,6,6,3,15,6,,0,2,,151
2-18-2015,87,75,63,68,64,60,94,72,32,30.09,30.03,29.98,6,6,4,16,5,,0.02,3,Fog-Rain,261
2-19-2015,68,62,57,66,53,46,88,66,48,30.22,30.15,30.06,6,6,6,18,12,,0,5,,3
2-20-2015,71,64,56,55,49,45,68,57,36,30.25,30.19,30.15,6,6,6,21,12,,0,3,,35
2-21-2015,78,71,64,66,61,55,94,68,45,30.21,30.16,30.12,6,6,6,23,14,,0,4,,72

[R-es] PCA output

2016-06-15 Thread Dr. José A. Betancourt Bethencourt
Estimados

�C�mo se puede exportar  los resultados de summary(res.pca)  a un archivo
que deseemos sea en pdf o htlm? Como hemos visto en otros an�lisis

 

saludos

 

library(FactoMineR)

res.pca<-PCA(X, scale.unit = TRUE, dimensiones ,  graph = TRUE, axes =
c(1,2))

summary(res.pca) 

 

 

X


Max,TemperatureF

Mean,TemperatureF

Min,TemperatureF

Max,Dew,PointF

MeanDew,PointF

Min,DewpointF

Max,Humidity

Mean,Humidity

Min,Humidity


1

82

75

68

62

58

54

72

56


2

82

74

67

61

58

56

74

58


3

84

75

66

62

58

52

77

58


4

82

73

65

61

57

54

76

59


5

83

75

66

60

58

54

73

57


6

84

73

62

64

59

56

83

64


7

85

75

65

61

58

54

76

58


8

85

76

66

60

56

49

73

53


9

84

74

64

60

55

49

72

53


10

84

75

65

60

57

52

72

54


11

84

75

65

61

56

49

74

53


12

87

76

64

63

55

50

72

52


13

86

76

64

61

57

50

75

54


14

86

77

67

59

56

51

70

51


15

85

76

68

62

60

57

76

58


16

87

77

67

66

61

56

81

59


17

89

77

66

68

62

56

84

63


18

88

77

67

68

61

55

85

62


19

88

79

70

65

63

60

78

61


20

85

76

68

66

62

58

85

66


21

84

76

68

67

64

60

88

69


22

85

76

66

68

63

57

92

70


23

81

73

66

70

66

64

98

83


24

85

76

67

68

64

61

88

70


25

85

77

70

67

64

62

82

65


26

84

77

70

68

64

62

85

66


27

84

76

68

67

64

62

84

68


28

82

74

66

68

64

60

90

73


29

83

74

65

69

63

57

96

73


30

83

76

68

66

63

60

82

67


31

85

77

69

65

62

58

82

62


32

85

76

68

67

65

61

90

72


33

84

75

67

68

63

59

90

67


34

84

75

66

68

63

57

91

69


35

85

76

66

69

64

59

95

71


36

84

75

66

68

63

58

88

68


37

82

74

66

71

65

62

100

78


38

84

75

66

70

66

62

99

81


39

85

75

64

70

66

61

93

79


40

84

75

66

70

65

60

93

76


41

84

74

65

69

65

61

93

77


42

82

74

67

70

67

64

95

82


43

81

74

66

70

66

62

94

81


44

83

73

64

70

65

60

94

78


45

83

74

65

70

65

61

93

80


46

82

75

66

71

66

63

94

80


47

81

73

65

67

64

60

88

75


48

84

74

64

64

60

56

79

62


49

84

76

67

66

62

59

83

65


50

84

75

65

69

64

61

89

72


51

85

76

67

65

62

58

80

63

 



--
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/




[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] RV: aov

2016-05-24 Thread Dr. José A. Betancourt Bethencourt
 

 

 

Estimados

 

¿Cómo se podría modificar el script o su source para lograr que en la salida
se vean de manera clara los números? En este ejemplo quedan parcialmente
afuera. Adjunto datos y script

Saludos

José

 

##



rm(list = ls())  

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



data<-read.csv("dengueutiuti10.csv",header=TRUE, sep=";", dec=",")

head(data)

depVar=data$plt

grpVar=data$bleed



library(sjPlot)



sjp.aov1(depVar, grpVar, meansums = TRUE, title = NULL,

 axisLabels.y = NULL, reverseOrder = FALSE,

 stringIntercept = "(Intercept)", showAxisLabels.y = TRUE,

 axisTitle.x = "", axisLimits = NULL, geom.colors = c("#3366a0",

 "#aa"), geom.size = 3, breakTitleAt = 25, breakLabelsAt = 25,

 gridBreaksAt = NULL, expand.grid = TRUE, showValueLabels = TRUE,

 labelDigits = 2, y.offset = 0.1, showPValueLabels = TRUE,

 showModelSummary = TRUE, printPlot = TRUE)





Mi correo anterior parece que no sirve ya

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

Enviado el: martes, 24 de mayo de 2016 06:57
Para: 'r-help-es@r-project.org' <r-help-es@r-project.org
<mailto:r-help-es@r-project.org> >
Asunto: aov





--
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/





script.r
Description: Binary data


dengueutiuti10.csv
Description: MS-Excel spreadsheet
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] RV: heplots

2016-03-07 Thread Dr. José A. Betancourt Bethencourt
 

 

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

Enviado el: lunes, 7 de marzo de 2016 05:40
Para: 'r-help-es@r-project.org' <r-help-es@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/





duda.csv
Description: MS-Excel spreadsheet
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] RV: fechas

2016-01-16 Thread Dr. José A. Betancourt Bethencourt
Disculpen, A veces me falla el otro correo y lo reenvío por este

 

De: R-help-es [mailto:r-help-es-boun...@r-project.org] En nombre de
jbetancourt
Enviado el: sábado, 16 de enero de 2016 05:49
Para: 'Lista R' 
CC: luisfo89 
Asunto: [R-es] fechas

 

 

Estimados

Muchas gracias por su ayuda. No he logrado convertir los datos de climas
diarios a datos por semana, es decir

 que junto a la identificación de las semanas, quisiera que los datos
climaticos de cada semana se reflejaran 

en cada semana .Adjunto datos

Ejemplo

semana48  Max.TemperatureF Mean.TemperatureF Min.TemperatureF etc

Perdonen mis limitados conocimientos

 

LOS SCRIPT ENVIADOS FUNCIONARON PARA IDENTIFICAR LA SEMANA

#

rm(list = ls()) 

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

datIn <- read.table("clima.csv", header=T, as.is=TRUE, sep=",")

datIn

library(lubridate)

#datIn$fecha <- mdy(datIn$CST)

datIn$semana <- week(datIn$fecha)

datIn$semana

##

library("lubridate")

data <- read.csv("clima.csv", stringsAsFactors=F)

# mdy() para formato m/d/y. existe dmy() para formato d/m/y

woy <- week(mdy(data$CST)) #week of year

woy

SEMANAS

1  1  1  1  1  1  1  2  2  2  2  2  2  2  3  3  3  3  3  3  3  4  4  4  4  4
4

 [28]  4  5  5  5  5  5  5  5  6  6  6  6  6  6  6  7  7  7  7  7  7  7  8
8  8  8  8

 [55]  8  8  9  9  9  9  9  9  9 10 10 10 10 10 10 10 11 11 11 11 11 11 11
12 12 12 12

 [82] 12 12 12 13 13 13 13 13 13 13 14 14 14 14 14 14 14 15 15 15 15 15 15
15 16 16 16

[109] 16 16 16 16 17 17 17 17 17 17 17 18 18 18 18 18 18 18 19 19 19 19 19
19 19 20 20

[136] 20 20 20 20 20 21 21 21 21 21 21 21 22 22 22 22 22 22 22 23 23 23 23
23 23 23 24

[163] 24 24 24 24 24 24 25 25 25 25 25 25 25 26 26 26 26 26 26 26 27 27 27
27 27 27 27

[190] 28 28 28 28 28 28 28 29 29 29 29 29 29 29 30 30 30 30 30 30 30 31 31
31 31 31 31

[217] 31 32 32 32 32 32 32 32 33 33 33 33 33 33 33 34 34 34 34 34 34 34 35
35 35 35 35

[244] 35 35 36 36 36 36 36 36 36 37 37 37 37 37 37 37 38 38 38 38 38 38 38
39 39 39 39

[271] 39 39 39 40 40 40 40 40 40 40 41 41 41 41 41 41 41 42 42 42 42 42 42
42 43 43 43

[298] 43 43 43 43 44 44 44 44 44 44 44 45 45 45 45 45 45 45 46 46 46 46 46
46 46 47 47

[325] 47 47 47 47 47 48 48 48 48 48 48 48 49 49 49 49 49 49 49 50 50 50 50
50 50 50 51

[352] 51 51 51 51 51 51 52 52 52 52 52 52 52 53



--
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/





clima.csv
Description: Binary data
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] RV: least square for time series

2016-01-09 Thread Dr. José A. Betancourt Bethencourt
 

 

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

Enviado el: sábado, 9 de enero de 2016 07:40
Para: 'r-help-es@r-project.org' <r-help-es@r-project.org>
Asunto: least square for time series

 

Estimados

Logré mejorar la predicción con un r2 del 73%. 

 

A este script quisiera añadirle un comando que me permita ver los datos  que
se pronostican al menos en las próximas 4 semanas y una evaluación del
modelo, a través del error . Adjunto serie de casos semanales desde enero a
noviembre. Apreciaría su ayuda.

saludos cordiales

Betancourt

##



rm(list=ls())

datos<-scan("data1.csv")



fit1<- lm (log10(datos)~time(datos)) # regress gtemp on time



x<-(log10(datos))

x

plot(x, type="l", ylab="dengue cases CR")

abline(fit1)

summary(fit1)





--
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/





data1.csv
Description: MS-Excel spreadsheet
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] RV: ANOVA R

2015-12-28 Thread Dr. José A. Betancourt Bethencourt
 

Estimados

He buscado un script para ANOVA en medidas repetidas en los paquetes aov y
easyanova y no los encuentro

�Podr�an enviarme un script para para ANOVA para medidas repetidas?

Saludos y feliz 2016

Jos� (mi correo usual no funciona ahora)



--
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/




[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] pregunta

2015-06-13 Thread Dr. José A. Betancourt Bethencourt
Gracias a todos!

 

De: Freddy Omar López Quintero [mailto:freddy.vat...@gmail.com] 
Enviado el: viernes, 12 de junio de 2015 11:06
Para: jbetancourt
CC: Lista R
Asunto: Re: [R-es] pregunta

 

​Holap.​

 

x - c 
(1,3,1,1,1,1,2,1,1,2,1,1,2,1,1,1,2,1,2,2,1,1,1,3,1,1,2,2,2,2,1,1,1,2,1,1,1,1,1,2,2,5,1,1,1,1,2,1,1,1,2,4,1,2,1,3,1,1,1,1,3,1,2,1,1,3,1,3,3,1,2,1,2,2,2,3,1,2,1,1,1,1,2,1,1,1,1,1,1)

 

names(x) - format(index(zooreg(1:89, start =as.Date(2015-01-01))), 
'%Y-%m-%d')

 

​Puedes usar:

 

x-scan('ruta/a/tu/archivo/a.txt')​

 

y luego el código que ya utilizas.

 

​¡Salud!​

 

-- 

«No soy aquellas sombras tutelares
que honré con versos que no olvida el tiempo.»

JL Borges



--
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/




[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es


[R-es] pregunta

2014-08-27 Thread Dr. José A Betancourt Bethencourt
 

Estimados

�Existe en R alguna manera de que las tablas de estad�stica descriptiva como
las de codebook de epicalc,  salgan en una tabla como las de SPSS  ?

 

Saludos

Jos�



--
Nunca digas nunca, di mejor: gracias, permiso, disculpe.

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/




[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es


[R-es] pregunta

2014-08-20 Thread Dr. José A Betancourt Bethencourt
Estimados

Estoy entrenando hacer funciones que respondan a comandos, 

 

en esta caso en la salida gráfica se observa que dice :  Exposure=var3  y
outcome=var 1 

 

quisiéramos que se reflejan los nombres de la base de datos : var1=estado,
var2=cake, var3=chocolate

 

Espero haberme explicado adecuadamente

Adjunto tabla con datos

 





#Comando que llama  a una función

rm(list=ls())

#setwd(D:/DEMO_new/demo_scripts/OR/) 

#setwd(D:/Public/Documents/R/EPICALC/funciones/OR/)



data= mydata-read.csv(OR.csv,header=TRUE, sep=,, dec=.)

use(data)

attach(data)



var1=estado

var2=cake

var3=chocolate

library(epicalc)

source(function_or.r)

odratios(data,var1,var2,var3)





#función

odratios - function (data,var1,var2,var3){

  or1 -cc(var1, var2)   

  or2 - cc(var1, var3)

}



--
Nunca digas nunca, di mejor: gracias, permiso, disculpe.

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/





OR.csv
Description: MS-Excel spreadsheet
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es