5.Simulacion estadistica del Ejercicio 4.8 (Recta de regresion de una variable bidimensional continua)

Programa en R-Project:

> library(Ryacas);
>

> yacas("mncentrado(k1,k2,lox,upx,loy,upy,f):=Integrate(x,lox,upx)( Integrate(y
,loy,upy) ((x^k1)*(y^k2)*f)))");

> 
> alpha10 <- Eval(yacas("mncentrado(1,0,0,1,0,1,x+y)"));
> 
> alpha01 <- Eval(yacas("mncentrado(0,1,0,1,0,1,x+y)"));

> 
> alpha11 <- Eval(yacas("mncentrado(1,1,0,1,0,1,x+y)"));
> 
> mu11 <- alpha11- alpha01*alpha10;
> 
> alpha20 <- Eval(yacas("mncentrado(2,0,0,1,0,1,x+y)"));
> 
> mu20 <- alpha20 - alpha10^2;
> 

> #Recta de regresion de la funcion de densidad x+y
> 
> linearegresion <- function(x){ alpha01+(mu11/mu20)*(x-alpha10)}
> 
> plot(linearegresion)