Top Banner
Webapp en Shiny para modelos de regresi´ on con polinomios fraccionarios de orden 1 Antonio Torres Ruiz Isabel Mar´ ıa Ortiz Rodr´ ıguez Universidad de Almer´ ıa IX JORNADAS DE USUARIOS DE R Granada, Noviembre de 2017
46

Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Oct 04, 2018

Download

Documents

nguyenkhue
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Webapp en Shiny para modelos de regresioncon polinomios fraccionarios de orden 1

Antonio Torres RuizIsabel Marıa Ortiz Rodrıguez

Universidad de Almerıa

IX JORNADAS DE USUARIOS DE RGranada, Noviembre de 2017

Page 2: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Indice

1 Polinomios fraccionarios

2 Modelos de regresion con polinomios fraccionarios de orden 1

3 Eleccion del mejor modelo fraccionario de orden 1

4 Webapp en Shiny para modelos fraccionarios de orden 1

5 Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacional

6 El paquete mfp de R

7 Conclusiones

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 3: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Caracterısticas del polinomio fraccionario de orden 1Caracterısticas del polinomio fraccionario de orden 2

Indice

1 Polinomios fraccionarios

2 Modelos de regresion con polinomios fraccionarios de orden 1

3 Eleccion del mejor modelo fraccionario de orden 1

4 Webapp en Shiny para modelos fraccionarios de orden 1

5 Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacional

6 El paquete mfp de R

7 Conclusiones

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 4: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Caracterısticas del polinomio fraccionario de orden 1Caracterısticas del polinomio fraccionario de orden 2

Polinomios fraccionarios

Un polinomio fraccionario de orden m se denota por PFm y vienedefinido como:

φm(x ;β, p) = β0 +m∑j=1

βjHj(x),

con Hj(x): H1(x) = x (p1) para j = 1

Hj(x) =

{x (pj ) si pj 6= pj−1

Hj−1(x) ln x si pj = pj−1para j = 2, ...,m

y

x (p) =

{xp p 6= 0ln x p = 0

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 5: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Caracterısticas del polinomio fraccionario de orden 1Caracterısticas del polinomio fraccionario de orden 2

Royston y Altman (1994): tomar la potencia p en el conjunto

P = {−2,−1,−0.5, 0, 0.5, 1, 2, 3}

Transformacion cuando x < 0:

φm(x ;β, p) = β0 +m∑j=1

βjHj(x − c)

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 6: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Caracterısticas del polinomio fraccionario de orden 1Caracterısticas del polinomio fraccionario de orden 2

Caracterısticas del polinomio fraccionario de orden 1

PF1(p) = β0 + β1x(p) =

{β0 + β1x

p p ∈ P, p 6= 0β0 + β1 ln x p = 0

Representacion para β0 = β1 = 1:

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 7: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Caracterısticas del polinomio fraccionario de orden 1Caracterısticas del polinomio fraccionario de orden 2

Caracterısticas del polinomio fraccionario de orden 2

PF2(p1, p2) =

{β0 + β1x

(p1) + β2x(p2) p1, p2 ∈ P, p1 < p2

β0 + β1x(p) + β2x

(p) ln x p1, p2 ∈ P, p1 = p2 = p

Comportamiento grafico para varios tipos de curvas:

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 8: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Hipotesis basicas sobre los errores aleatoriosResiduos del modeloEstimacion de los parametros del modeloI.C., Contrastes y Bondad del ajusteComprobacion de las hipotesis del modelo

Indice

1 Polinomios fraccionarios

2 Modelos de regresion con polinomios fraccionarios de orden 1

3 Eleccion del mejor modelo fraccionario de orden 1

4 Webapp en Shiny para modelos fraccionarios de orden 1

5 Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacional

6 El paquete mfp de R

7 Conclusiones

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 9: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Hipotesis basicas sobre los errores aleatoriosResiduos del modeloEstimacion de los parametros del modeloI.C., Contrastes y Bondad del ajusteComprobacion de las hipotesis del modelo

Modelos de regresion con PF1

Modelo de regresion con polinomios fraccionarios de orden 1:

Y = β0 + β1X(p) + ε, X (p) =

{X p p 6= 0lnX p = 0

,

con p ∈ P = {−2,−1,−0.5, 0, 0.5, 1, 2, 3}.

Para una muestra de tamano n:

Yi = β0 + β1X(p)i + εi , i = 1, ..., n.

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 10: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Hipotesis basicas sobre los errores aleatoriosResiduos del modeloEstimacion de los parametros del modeloI.C., Contrastes y Bondad del ajusteComprobacion de las hipotesis del modelo

Hipotesis basicas sobre los errores aleatorios εi

Yi = β0 + β1X(p)i + εi , i = 1, ..., n.

E (εi ) = 0, i = 1, ..., n.

Var (εi ) = E(ε2i

)= σ2 = cte, i = 1, ..., n.

Cov (εi , εj) = E (εiεj) = 0, para todo i 6= j .

Hipotesis adicional: εi ∼ N (0, σ2), i = 1, ..., n.

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 11: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Hipotesis basicas sobre los errores aleatoriosResiduos del modeloEstimacion de los parametros del modeloI.C., Contrastes y Bondad del ajusteComprobacion de las hipotesis del modelo

Residuos del modelo

Muestra de tamano n: (Xi ,Yi ), i = 1, ..., n.

Valores observados: Yi , i = 1, ..., n.

Modelo estimado: Y = β0 + β1X(p).

Valor estimado para Xi : Yi = β0 + β1X(p)i .

Residuos: ei = Yi − Yi , i = 1, ..., n.

A partir de ahora,X (p) = Z

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 12: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Hipotesis basicas sobre los errores aleatoriosResiduos del modeloEstimacion de los parametros del modeloI.C., Contrastes y Bondad del ajusteComprobacion de las hipotesis del modelo

Estimacion de los parametros del modelo

Estimacion de los parametros β0 y β1:{β0 = Y − Syz

S2zzZ

β1 =Syz

S2zz

con

Syz =

∑ni=1 YiZi

n− Y Z

S2zz =

∑ni=1 Z

2i

n− Z 2

Estimador insesgado para σ2:

σ2 = S2R =

∑ni=1 e

2i

n−2

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 13: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Hipotesis basicas sobre los errores aleatoriosResiduos del modeloEstimacion de los parametros del modeloI.C., Contrastes y Bondad del ajusteComprobacion de las hipotesis del modelo

Intervalos de confianza 100(1− α) %

Para β0:β0 − tn−2,1−α/2

√S2R

n

(1 +

Z 2

S2zz

), β0 + tn−2,1−α/2

√S2R

n

(1 +

Z 2

S2zz

)Para β1: β1 − tn−2,1−α/2

√S2R

nS2zz

, β1 + tn−2,1−α/2

√S2R

nS2zz

Para σ2: [

(n − 2)S2R

χ2n−2,1−α/2

,(n − 2)S2

R

χ2n−2,α/2

]Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 14: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Hipotesis basicas sobre los errores aleatoriosResiduos del modeloEstimacion de los parametros del modeloI.C., Contrastes y Bondad del ajusteComprobacion de las hipotesis del modelo

Contrastes de hipotesis

Para β0 y β1:

Parametro β0 β1

Contraste

{H0 : β0 = β00

H1 : β0 6= β00

{H0 : β1 = β10

H1 : β1 6= β10

Estadıstico T = β0−β00√S2Rn

(1+ Z2

S2zz

) ∼ tn−2 T = β1−β10√S2R

nS2zz

∼ tn−2

Region de rechazo |T | > tn−2,1−α/2 |T | > tn−2,1−α/2

Contraste H0 : β1 = 0. Tabla ANOVA:

Fuente de variacion SC gl MC Estadıstico

Regresion SCR 1 MCR FError SCE n − 2 MCETotal SCT n − 1

Rechazar H0 si F > F1,n−2,1−α.

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 15: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Hipotesis basicas sobre los errores aleatoriosResiduos del modeloEstimacion de los parametros del modeloI.C., Contrastes y Bondad del ajusteComprobacion de las hipotesis del modelo

Intervalos de prediccion

Para la respuesta media en X = X0 (Z0 = X(p)0 ):Y0 − tn−2,1−α/2

√S2R

(1

n+

(Z0 − Z)2

nS2zz

), Y0 + tn−2,1−α/2

√S2R

(1

n+

(Z0 − Z)2

nS2zz

)

Para la prediccion en X = X0 (Z0 = X(p)0 ):

Y0 − tn−2,1−α/2

√√√√S2R

(1 +

1

n+

(Z0 − Z)2

nS2zz

), Y0 + tn−2,1−α/2

√√√√S2R

(1 +

1

n+

(Z0 − Z)2

nS2zz

)

Bandas de confianza y de prediccion: unir los extremos de losintervalos anteriores calculados para cada valor de X .

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 16: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Hipotesis basicas sobre los errores aleatoriosResiduos del modeloEstimacion de los parametros del modeloI.C., Contrastes y Bondad del ajusteComprobacion de las hipotesis del modelo

Coeficiente de determinacion R2

R2 =SCR

SCT= 1− SCE

SCT

0 ≤ R2 ≤ 1.

Cuanto mayor sea R2 mejor sera el modelo estimado.

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 17: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Hipotesis basicas sobre los errores aleatoriosResiduos del modeloEstimacion de los parametros del modeloI.C., Contrastes y Bondad del ajusteComprobacion de las hipotesis del modelo

Comprobacion de las hipotesis del modelo

Homocedasticidad y linealidad Grafico qq-plot

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 18: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Indice

1 Polinomios fraccionarios

2 Modelos de regresion con polinomios fraccionarios de orden 1

3 Eleccion del mejor modelo fraccionario de orden 1

4 Webapp en Shiny para modelos fraccionarios de orden 1

5 Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacional

6 El paquete mfp de R

7 Conclusiones

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 19: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Eleccion del mejor modelo PF1

Medidas para comparar los modelos:

Medida Expresion Criterio

SCE SCE ↘MCE SCE

n−2 ↘R2 1− SCE

SCT ↗

R2 1−(1− R2

) n − 1

n − 2↗

AIC n ln(

SCEn

)+ n ln(2π) + n + 6 ↘

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 20: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

La aplicacion desarrollada

Indice

1 Polinomios fraccionarios

2 Modelos de regresion con polinomios fraccionarios de orden 1

3 Eleccion del mejor modelo fraccionario de orden 1

4 Webapp en Shiny para modelos fraccionarios de orden 1

5 Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacional

6 El paquete mfp de R

7 Conclusiones

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 21: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

La aplicacion desarrollada

La aplicacion desarrollada

Enlace de la aplicacion:

https://antoniotorres.shinyapps.io/shiny/

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 22: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Indice

1 Polinomios fraccionarios

2 Modelos de regresion con polinomios fraccionarios de orden 1

3 Eleccion del mejor modelo fraccionario de orden 1

4 Webapp en Shiny para modelos fraccionarios de orden 1

5 Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacional

6 El paquete mfp de R

7 Conclusiones

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 23: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Relacion entre edad gestacional (X , semanas) y longitud de la mandıbula (Y , mm)

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 24: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Introduccion de los datos (Chitty et al. 1994):

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 25: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Eleccion del mejor modelo:

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 26: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Comprobacion de las hipotesis basicas:

El modelo considerado no es el adecuado.

Y = ln(Longitud)

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 27: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Relacion entre edad gestacional (X ) y Y = ln(Longitud)

Introduccion de los datos:

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 28: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Eleccion del mejor modelo:

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 29: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Comprobacion de las hipotesis basicas:

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 30: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Representacion grafica de las observaciones y el modelo estimado:

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 31: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Estimacion de parametros, I.C. y contrastes:

Modelo estimado: ln(Longitud) = 4.69− 30.41

SemanasAntonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 32: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Bandas de confianza y prediccion:

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 33: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Modelo con datos originales:

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 34: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Indice

1 Polinomios fraccionarios

2 Modelos de regresion con polinomios fraccionarios de orden 1

3 Eleccion del mejor modelo fraccionario de orden 1

4 Webapp en Shiny para modelos fraccionarios de orden 1

5 Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacional

6 El paquete mfp de R

7 Conclusiones

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 35: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

El paquete mfp (Sauerbrei et al. 2006)

Contraste: {H0 : Modelo restringidoH1 : Modelo completo

mediante test de razon de verosimilitud:

RV = −2 lnL0

L1

∼ χ2r

Ejemplo:{H0 : Modelo lineal (p = 1)H1 : Modelo no lineal (p 6= 1)

⇔{

H0 : E(Y ) = β0 + β1XH1 : E(Y ) = β0 + β1X

(p), p 6= 1

⇒ RV ∼ χ23−2 = χ2

1

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 36: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Esquema del proceso iterativo:

Relacion entre X e Y :{H0 : Modelo nuloH1 : Mejor PF2

Linealidad:{H0 : Modelo linealH1 : Mejor PF2

Simplificacion:

{H0 : Mejor PF1

H1 : Mejor PF2

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 37: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Aplicacion del algoritmo mfp al ejemplo

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 38: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Indice

1 Polinomios fraccionarios

2 Modelos de regresion con polinomios fraccionarios de orden 1

3 Eleccion del mejor modelo fraccionario de orden 1

4 Webapp en Shiny para modelos fraccionarios de orden 1

5 Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacional

6 El paquete mfp de R

7 Conclusiones

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 39: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 40: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 41: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 42: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 43: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 44: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Bibliografıa I

L.C. AcostaAjuste de polinomios fraccionarios a curvas de crecimientoDisponible en: Enlace al texto completo. [Consulta 13 nov. 2017]

J. Aparicio, M.A. Martınez, J. MoralesModelos lineales aplicados en RDisponible en: Enlace al texto completo. [Consulta 13 nov. 2017]

J.E. CallejasMetodos de aproximacion polinomica en el diseno experimentalDisponible en: Enlace al texto completo. [Consulta 13 nov. 2017]

L.S. Chitty, D.G. Altman, A. Henderson, S. CampbellCharts of fetal size: 3, Abdominal measurementsBr. J. Obstet. Gynaecol., 101 (1994), 125–131

D. C. Montgomery, E. A. Peck, G. G. ViningIntroduction to linear regression analysisWiley, 1992

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 45: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

Bibliografıa II

P. Royston, D.G. AltmanRegression using fractional polynomials of continuous covariatesJ. Royal Stat. Soc., C 43 (1994), 429–467

W. Sauerbrei, C. Meier-Hirmer, A. Benner, P. RoystonMultivariable regression model building by using fractional polynomials.Comput. Stat. Data An., 50 (2006), 3464–3485

Shiny: A web application framework for RDisponible en: Enlace web. [Consulta 13 nov. 2017]

J. ZamoraUna metodologıa para manejar variables continuas en los modelos depronostico: Polinomios fraccionales

Disponible en: Enlace al texto completo. [Consulta 13 nov. 2017]

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion

Page 46: Webapp en Shiny para modelos de regresi on - r-es.orgr-es.org/9jornadasR/pdf/9JUR_paper_6.pdf · I.C., Contrastes y Bondad del ajuste Comprobacion de las hip otesis del modelo Modelos

Polinomios fraccionariosModelos de regresion con polinomios fraccionarios de orden 1

Eleccion del mejor modelo fraccionario de orden 1Webapp en Shiny para modelos fraccionarios de orden 1

Ejemplo: Relacion de Longitud de la mandıbula con Edad gestacionalEl paquete mfp de R

Conclusiones

¡MUCHAS GRACIAS!

Antonio Torres, Isabel M. Ortiz Webapp en Shiny para modelos de regresion