22 jun 2011
Estadística en Java
Como para ejecutarla desde el celular, esta chevere. http://www.socr.ucla.edu/htmls/SOCR_Analyses.html
15 jun 2011
Cardioide en R.
# Figura 2: Generacion de una cardioide.
# Construccion tomada de: http://mathworld.wolfram.com/Cardioid.html
library(grid)
ori <- c(7.5,8.5)
r <- 2.5
grid.circle(x=ori[1], y=ori[2], r=r, default.units="cm",
gp=gpar(col="red", lwd=2) )
a <- ori+c(0,r)
tet <- seq(pi/2,5*pi/2,length=40)
x <- r*cos(tet)+ori[1]
y <- r*sin(tet)+ori[2]
rad <- sqrt( (x-a[1])^2+(y-a[2])^2 )
grid.circle(x=x, y=y, r=rad, default.units="cm", gp=gpar(col="blue") )
# Construccion tomada de: http://mathworld.wolfram.com/Cardioid.html
library(grid)
ori <- c(7.5,8.5)
r <- 2.5
grid.circle(x=ori[1], y=ori[2], r=r, default.units="cm",
gp=gpar(col="red", lwd=2) )
a <- ori+c(0,r)
tet <- seq(pi/2,5*pi/2,length=40)
x <- r*cos(tet)+ori[1]
y <- r*sin(tet)+ori[2]
rad <- sqrt( (x-a[1])^2+(y-a[2])^2 )
grid.circle(x=x, y=y, r=rad, default.units="cm", gp=gpar(col="blue") )
7 jun 2011
5 jun 2011
Función Para crear un Gráfico Cuantil-Cuantil con ggplot2
install.packages("ggplot2")
library(ggplot2)
ggQQ <- function(LM) # argument: a linear model
{
y <- quantile(LM$resid[!is.na(LM$resid)], c(0.25, 0.75))
x <- qnorm(c(0.25, 0.75))
slope <- diff(y)/diff(x)
int <- y[1L] - slope * x[1L]
p <- ggplot(LM, aes(sample=.resid)) +
stat_qq(alpha = 0.5) +
geom_abline(slope = slope, intercept = int, color="blue")+
opts(title = "Gráfico Cuantil-Cuantil")+ylab("Muestrales")+
xlab("Teoricos")
return(p)
}
ggQQ(LM)
library(ggplot2)
ggQQ <- function(LM) # argument: a linear model
{
y <- quantile(LM$resid[!is.na(LM$resid)], c(0.25, 0.75))
x <- qnorm(c(0.25, 0.75))
slope <- diff(y)/diff(x)
int <- y[1L] - slope * x[1L]
p <- ggplot(LM, aes(sample=.resid)) +
stat_qq(alpha = 0.5) +
geom_abline(slope = slope, intercept = int, color="blue")+
opts(title = "Gráfico Cuantil-Cuantil")+ylab("Muestrales")+
xlab("Teoricos")
return(p)
}
ggQQ(LM)
Suscribirse a:
Entradas
(
Atom
)