Parcial
Sebastián Pulido Almeida
2 de noviembre de 2015
Parcial: Punto 1:
set.seed(1234)
f<-function(x,a,b,c){
(exp(-((abs(x-a)/b)^(2/c))*(1/2)))/(b*(2^((c/2)+1))*gamma(1+(c/2)))
}
gri<-seq(-1000000,1000000,0.1)
v<-f(gri,0,47500,2)
w<-v/sum(v)
X<-sample(gri,1000,prob=w, replace=TRUE)
hist(X,freq=FALSE)
curve(f(x,0,47500,2),add=TRUE)
Punto 2:
set.seed(1234)
f<-function(x,a,b,c){
(exp(-((abs(x-a)/b)^(2/c))*(1/2)))/(b*(2^((c/2)+1))*gamma(1+(c/2)))
}
gri<-seq(-1000000,1000000,0.1)
v<-f(gri,0,47500,2)
w<-v/sum(v)
X<-sample(gri,1000,prob=w, replace=TRUE)
hist(X,freq=FALSE)
curve((exp(x/95000))/(2*95000),add=TRUE)
curve((exp(-x/95000))/(2*95000),add=TRUE)