¿De dónde se tomó la idea?

La idea para esta publicación se tomó de un twitter de Christophe Dervieux @chrisderv, abajo el twitter.

Botón para lanzar el kernel de R

Lo primero que se debe hacer es crear un botón para que el usuario pueda iniciar el kernel de R dentro del documento html.

Por favor dé clic en el botón de abajo para activar el kernel de R, puede tomar unos segundos, sea paciente.

Ejemplo

En este ejemplo el lector va a simular \(n=100\) observaciones de un modelo de regresión lineal simple.

Por favor ingrese abajo dos números para \(\beta_0\) y \(\beta_1\), luego oprima run.

beta0 <- Aquí
beta1 <- Aquí
n <- 100
x <- runif(n)
y <- rnorm(n=n, mean=beta0 + beta1 * x, sd=0.1)

Oprima run para dibujar los datos simulados.

plot(x=x, y=y)

Oprima run para estimar los parámetros \(\beta_0\) y \(\beta_1\) del modelo simulado.

mod <- lm(y ~ x)
coef(mod)

¿Logra ver que las estimaciones son muy similares a los parámetros?

Posibilidades

Hay muchas posibilidades de sacarle provecho a esto, ¿cuál sería una de ellas?

LS0tDQp0aXRsZTogIkVqZW1wbG8gZGUgY8OzZGlnbyBkZSBSIGludGVyYWN0aXZvIGRlbnRybyBkZSB1biBkb2N1bWVudG8gaHRtbCINCmF1dGhvcjogIkZyZWRkeSBIZXJuw6FuZGV6Ig0KZGF0ZTogIjA2LzAzLzIwMjAiDQpvdXRwdXQ6IA0KICBodG1sX2RvY3VtZW50OiANCiAgICBjb2RlX2Rvd25sb2FkOiB5ZXMNCiAgICB0aGVtZTogc2ltcGxleA0KLS0tDQoNCjwhLS0gU29tZSBzdHVmZiBmcm9tIHRoZWJlbGFiIHJlcXVpcmVkIGZvciB0aGUgbWFnaWMgLS0+DQo8IS0tIENvbmZpZ3VyZSBhbmQgbG9hZCBUaGViZSAhLS0+DQo8c2NyaXB0IHR5cGU9InRleHQveC10aGViZS1jb25maWciPg0KICB7DQogICAgcmVxdWVzdEtlcm5lbDogdHJ1ZSwNCiAgICBiaW5kZXJPcHRpb25zOiB7DQogICAgICByZXBvOiAiYmluZGVyLWV4YW1wbGVzL3IiLA0KICAgIH0sDQogICAga2VybmVsT3B0aW9uczogew0KICAgICAgbmFtZTogIlIiLA0KICAgICAga2VybmVsTmFtZTogImlyIiwNCiAgICB9LA0KICB9DQo8L3NjcmlwdD4NCg0KPCEtLSBzY3JpcHQgZm9yIHRoZWJlbGFiIC0tPg0KPHNjcmlwdCBzcmM9Imh0dHBzOi8vdW5wa2cuY29tL3RoZWJlbGFiQGxhdGVzdC9saWIvaW5kZXguanMiPjwvc2NyaXB0Pg0KDQpgYGB7ciBzZXR1cCwgaW5jbHVkZSA9IEZBTFNFfQ0Ka25pdHI6Om9wdHNfY2h1bmskc2V0KGV2YWwgPSBGQUxTRSkNCmBgYA0KDQojIMK/RGUgZMOzbmRlIHNlIHRvbcOzIGxhIGlkZWE/DQpMYSBpZGVhIHBhcmEgZXN0YSBwdWJsaWNhY2nDs24gc2UgdG9tw7MgZGUgdW4gdHdpdHRlciBkZSBDaHJpc3RvcGhlIERlcnZpZXV4IEBjaHJpc2RlcnYsIGFiYWpvIGVsIHR3aXR0ZXIuDQoNCjxibG9ja3F1b3RlIGNsYXNzPSJ0d2l0dGVyLXR3ZWV0Ij48cCBsYW5nPSJlbiIgZGlyPSJsdHIiPjxhIGhyZWY9Imh0dHBzOi8vdHdpdHRlci5jb20vaGFzaHRhZy9yc3RhdHM/c3JjPWhhc2gmYW1wO3JlZl9zcmM9dHdzcmMlNUV0ZnciPiNyc3RhdHM8L2E+IHByb2plY3Qgb2YgdGhlIGRheSA6IG1ha2UgUiBjb2RlIGNodW5rcyBpbnRlcmFjdGl2ZSBpbiBodG1sIG91dHB1dCBmcm9tIGEgPGEgaHJlZj0iaHR0cHM6Ly90d2l0dGVyLmNvbS9oYXNodGFnL3JtYXJrZG93bj9zcmM9aGFzaCZhbXA7cmVmX3NyYz10d3NyYyU1RXRmdyI+I3JtYXJrZG93bjwvYT4gZmlsZS4gPGJyPk5leHQgdHJ5OiBpbnRlZ3JhdGUgd2l0aCA8YSBocmVmPSJodHRwczovL3R3aXR0ZXIuY29tL2hhc2h0YWcvYm9va2Rvd24/c3JjPWhhc2gmYW1wO3JlZl9zcmM9dHdzcmMlNUV0ZnciPiNib29rZG93bjwvYT4gZm9yIGludGVyYWN0aXZlIGV4YW1wbGUgY2h1bmtzLiBRdWljayBEZW1vIGRvYzogPGEgaHJlZj0iaHR0cHM6Ly90LmNvL3RoOTVWeEo5UVkiPmh0dHBzOi8vdC5jby90aDk1VnhKOVFZPC9hPjxicj5XaGF0IGRvIHlvdSB0aGluayBhYm91dCB0aGlzIHNtYWxsIGV4YW1wbGUgPyBJcyB0aGlzIG9mIGFueSBpbnRlcmVzdCA/IDxhIGhyZWY9Imh0dHBzOi8vdC5jby9OeWpsblZlcVUyIj5waWMudHdpdHRlci5jb20vTnlqbG5WZXFVMjwvYT48L3A+Jm1kYXNoOyBDaHJpc3RvcGhlIERlcnZpZXV4IChAY2hyaXNkZXJ2KSA8YSBocmVmPSJodHRwczovL3R3aXR0ZXIuY29tL2NocmlzZGVydi9zdGF0dXMvMTIzNTY2ODE2Mzc0MTgzOTM3Mz9yZWZfc3JjPXR3c3JjJTVFdGZ3Ij5NYXJjaCA1LCAyMDIwPC9hPjwvYmxvY2txdW90ZT4gPHNjcmlwdCBhc3luYyBzcmM9Imh0dHBzOi8vcGxhdGZvcm0udHdpdHRlci5jb20vd2lkZ2V0cy5qcyIgY2hhcnNldD0idXRmLTgiPjwvc2NyaXB0Pg0KDQojIEJvdMOzbiBwYXJhIGxhbnphciBlbCBrZXJuZWwgZGUgUg0KTG8gcHJpbWVybyBxdWUgc2UgZGViZSBoYWNlciBlcyBjcmVhciB1biBib3TDs24gcGFyYSBxdWUgZWwgdXN1YXJpbyBwdWVkYSBpbmljaWFyIGVsIGtlcm5lbCBkZSBSIGRlbnRybyBkZWwgZG9jdW1lbnRvIGh0bWwuDQoNCj4gUG9yIGZhdm9yIGTDqSBjbGljIGVuIGVsIGJvdMOzbiBkZSBhYmFqbyBwYXJhIGFjdGl2YXIgZWwga2VybmVsIGRlIFIsIHB1ZWRlIHRvbWFyIHVub3Mgc2VndW5kb3MsIDxtYXJrPnNlYSBwYWNpZW50ZTwvbWFyaz4uDQoNCmBgYHtyLCBldmFsID0gVFJVRSwgZWNobyA9IEZBTFNFfQ0KIyBhZCBhIGJ1dHRvbg0KaHRtbHRvb2xzOjp0YWdzJGJ1dHRvbigNCiAgaWQgPSAiYWN0aXZhdGVCdXR0b24iLA0KICBzdHlsZT0iY29sb3I6ICNmZmY7IGJhY2tncm91bmQtY29sb3I6ICNlOTU0MjA7IGJvcmRlci1jb2xvcjogI2MzNDExMzsgd2lkdGg6IDI1MHB4OyBoZWlnaHQ6IDUwcHg7IGZvbnQtc2l6ZTogMS41ZW07IiwgDQogICJBY3RpdmFyIGVsIGtlcm5lbCBkZSBSIg0KKQ0KYGBgDQoNCmBgYHtqcywgZXZhbD1UUlVFLCBlY2hvID0gRkFMU0V9DQovLyB0aGViZWxhYiBqcyBzY3JpcHQNCnZhciBib290c3RyYXBUaGViZSA9IGZ1bmN0aW9uKCkgew0KICAgIHRoZWJlbGFiLmJvb3RzdHJhcCgpOw0KfQ0KZG9jdW1lbnQucXVlcnlTZWxlY3RvcigiI2FjdGl2YXRlQnV0dG9uIikuYWRkRXZlbnRMaXN0ZW5lcignY2xpY2snLCBib290c3RyYXBUaGViZSkNCmBgYA0KDQojIEVqZW1wbG8NCkVuIGVzdGUgZWplbXBsbyBlbCBsZWN0b3IgdmEgYSBzaW11bGFyICRuPTEwMCQgb2JzZXJ2YWNpb25lcyBkZSB1biBtb2RlbG8gZGUgcmVncmVzacOzbiBsaW5lYWwgc2ltcGxlLiANCg0KUG9yIGZhdm9yIDxtYXJrPmluZ3Jlc2UgYWJham8gZG9zIG7Dum1lcm9zPC9tYXJrPiBwYXJhICRcYmV0YV8wJCB5ICRcYmV0YV8xJCwgbHVlZ28gb3ByaW1hIHJ1bi4NCg0KYGBge3IsIGF0dHIuc291cmNlID0gYygnZGF0YS1leGVjdXRhYmxlPSJ0cnVlIicsICdkYXRhLWxhbmd1YWdlPSJyIicpfQ0KYmV0YTAgPC0gQXF1w60NCmJldGExIDwtIEFxdcOtDQpuIDwtIDEwMA0KeCA8LSBydW5pZihuKQ0KeSA8LSBybm9ybShuPW4sIG1lYW49YmV0YTAgKyBiZXRhMSAqIHgsIHNkPTAuMSkNCmBgYA0KDQpPcHJpbWEgcnVuIHBhcmEgZGlidWphciBsb3MgZGF0b3Mgc2ltdWxhZG9zLg0KDQpgYGB7ciwgYXR0ci5zb3VyY2UgPSBjKCdkYXRhLWV4ZWN1dGFibGU9InRydWUiJywgJ2RhdGEtbGFuZ3VhZ2U9InIiJyl9DQpwbG90KHg9eCwgeT15KQ0KYGBgDQoNCk9wcmltYSBydW4gcGFyYSBlc3RpbWFyIGxvcyBwYXLDoW1ldHJvcyAkXGJldGFfMCQgeSAkXGJldGFfMSQgZGVsIG1vZGVsbyBzaW11bGFkby4NCg0KYGBge3IsIGF0dHIuc291cmNlID0gYygnZGF0YS1leGVjdXRhYmxlPSJ0cnVlIicsICdkYXRhLWxhbmd1YWdlPSJyIicpfQ0KbW9kIDwtIGxtKHkgfiB4KQ0KY29lZihtb2QpDQpgYGANCg0Kwr9Mb2dyYSB2ZXIgcXVlIGxhcyBlc3RpbWFjaW9uZXMgc29uIG11eSBzaW1pbGFyZXMgYSBsb3MgcGFyw6FtZXRyb3M/DQoNCiMgUG9zaWJpbGlkYWRlcw0KSGF5IG11Y2hhcyBwb3NpYmlsaWRhZGVzIGRlIHNhY2FybGUgcHJvdmVjaG8gYSBlc3RvLCDCv2N1w6FsIHNlcsOtYSB1bmEgZGUgZWxsYXM/DQoNCjxwIGFsaWduPSJjZW50ZXIiPg0KICA8aW1nIHNyYz0idGhpbmsuZ2lmIiB3aWR0aD0iMjAwIj4NCjwvcD4NCg0KDQoNCg==