path="/Users/mauriciogutierrez/Library/Mobile Documents/com~apple~CloudDocs/estadisticas";
pathfunctions=paste(path,"/general_functions_new.R",sep="")
pathdatabase=paste(path,"/Gabriela_Mintegui/fiv.xlsx",sep="")
library(gplots)
##
## Attaching package: 'gplots'
## The following object is masked from 'package:stats':
##
## lowess
library(ROCR)
source(pathfunctions)
## Loading required package: XLConnectJars
## XLConnect 0.2-12 by Mirai Solutions GmbH [aut],
## Martin Studer [cre],
## The Apache Software Foundation [ctb, cph] (Apache POI, Apache Commons
## Codec),
## Stephen Colebourne [ctb, cph] (Joda-Time Java library),
## Graph Builder [ctb, cph] (Curvesapi Java library)
## http://www.mirai-solutions.com ,
## http://miraisolutions.wordpress.com
data=readexcel(pathdatabase,"Hoja1")
means(data,data$edad)
##
## edad
## n 76
## Media 36.25
## Mediana 37
## DS 4.69
## Min 26
## Max 45
## Rango 19
## Q1 33
## Q3 40
## RIQ 7
## Normal 0.031
## ICI 35.18
## ICS 37.32
## Resumen Mediana[Q1-Q3] -->37(33-40)
means(data,data$FSH)
##
## FSH
## n 75
## Media 7.72
## Mediana 7.12
## DS 2.65
## Min 2.53
## Max 16
## Rango 13.47
## Q1 5.76
## Q3 9.2
## RIQ 3.44
## Normal 0.059
## ICI 7.11
## ICS 8.33
## Resumen Media[ICI-ICS] --> 7.72(7.11-8.33)
means(data,data$HAM)
##
## HAM
## n 67
## Media 2.8
## Mediana 1.41
## DS 3.22
## Min 0.19
## Max 15.3
## Rango 15.11
## Q1 0.63
## Q3 3.68
## RIQ 3.05
## Normal 0
## ICI 2.01
## ICS 3.58
## Resumen Mediana[Q1-Q3] -->1.41(0.63-3.68)
means(data,data$RCA)
##
## RCA
## n 34
## Media 9.74
## Mediana 10
## DS 5.87
## Min 0
## Max 30
## Rango 30
## Q1 5
## Q3 11
## RIQ 6
## Normal 0.01
## ICI 7.69
## ICS 11.78
## Resumen Mediana[Q1-Q3] -->10(5-11)
means(data,data$ovocit.met2)
##
## ovocit.met2
## n 76
## Media 6.62
## Mediana 5
## DS 4.33
## Min 0
## Max 22
## Rango 22
## Q1 4
## Q3 10
## RIQ 6
## Normal 0
## ICI 5.63
## ICS 7.61
## Resumen Mediana[Q1-Q3] -->5(4-10)
lineal(data$ovocit.met2,data$edad)
##
## Call:
## lm(formula = y ~ x)
##
## Coefficients:
## (Intercept) x
## 22.6873 -0.4433
##
## [1] "r2=0.23"
## [1] "p=0"
lineal(data$ovocit.met2,data$FSH)
##
## Call:
## lm(formula = y ~ x)
##
## Coefficients:
## (Intercept) x
## 11.3567 -0.6091
##
## [1] "r2=0.138"
## [1] "p=0.001"
lineal(data$ovocit.met2,data$HAM)
##
## Call:
## lm(formula = y ~ x)
##
## Coefficients:
## (Intercept) x
## 3.9953 0.9832
##
## [1] "r2=0.509"
## [1] "p=0"
lineal(data$ovocit.met2,data$RCA)
##
## Call:
## lm(formula = y ~ x)
##
## Coefficients:
## (Intercept) x
## 3.7062 0.2749
##
## [1] "r2=0.144"
## [1] "p=0.027"
reg=(lm(data$ovocit.met2 ~ data$edad + data$FSH + data$HAM + data$RCA))
summary(reg)
##
## Call:
## lm(formula = data$ovocit.met2 ~ data$edad + data$FSH + data$HAM +
## data$RCA)
##
## Residuals:
## Min 1Q Median 3Q Max
## -5.7643 -1.0221 0.2044 1.1492 7.5121
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 14.56625 5.48316 2.657 0.01309 *
## data$edad -0.19198 0.12058 -1.592 0.12301
## data$FSH -0.38306 0.20964 -1.827 0.07873 .
## data$HAM 0.70103 0.23486 2.985 0.00596 **
## data$RCA -0.01334 0.10335 -0.129 0.89826
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 2.804 on 27 degrees of freedom
## (44 observations deleted due to missingness)
## Multiple R-squared: 0.6306, Adjusted R-squared: 0.5758
## F-statistic: 11.52 on 4 and 27 DF, p-value: 1.38e-05
plot(data[,c(4:6)])
data(data)