
Disculpen por demorar en responder, aun tengo el problema de las randomizaciones, sin embargo para al final realizar lo que necesito me he dado cuenta queno es necesario hacer con todas las posibilidades, sin embargo tengo una duda con un grafico de contorno de probabilidad.
Usando el package car, he podido hacer lo siguiente:
Datos 1:
#####################################x=rnorm(100)y=1+0.3*x+0.3*rnorm(100)dataEllipse(x ,y, levels=0.80, center.pch=F, plot.points=F, add=F)###################################################
Obtengo un grafico sin mostrar los valores de x e y pero si la region de la probabilidad indicada,sin embargo me gustaria introducir en ese grafico otros dos pares de datos para ver si es que se ubican dentro de la region de probabilidad
########################################################data.frame(x=c(1.5,0.2),y=c(-0.2,2.3),row.names=c("A","B"))###########################################################
Aca tengo dos pares de datos, uno de nombre de fila "A" y otro de nombre de fila "B"Me gustaria que estos dos pares de datos aparescan dentro del grafico con una etiqueta para cada uno que diga "A" y "B" para reconocer cual es cual
Mi segunda duda es poder superponer dos lineas de contorno de probabilidad, con el mismo levels de 0.80, sin embargo que estas dos provengan de un origen de datos distintos
Por ejemplo mi linea de contorno que proviene de Datos 1 y otra de 
Datos 2:
################################################x=rnorm(100)y=1+0.4*x+0.4*rnorm(100)dataEllipse(x, y, levels=0.80, center.pch=F, add=F)################################################
y ambas ubicarlas en un mismo grafico y poner una etiqueta a cada linea de probabilidad, una que diga "Datos1" y otra "Datos2" correspondiente a cada origen de datos que le pertence a cada una
Agradesco de antemano por su atencion



Gabriel Antonio Trujillo Paucar
Estudiante de Biologia
     997603768



Date: Mon, 14 Jan 2013 18:54:56 -0300
Subject: Re: [R-es] Generacion de numeros al azar en data.frame
From: daniel319@gmail.com
To: gabriel_enter@hotmail.com
CC: javier.ruben.marcuzzi@gmail.com; r-help-es@r-project.org

Permíteme repasar tu mensaje, soy algo lento y eso me ayuda.
# Tu tienes una matrix  diagonal de 300 x 300 con obviamente 300 unos. Seguramente solo a titulo ilustrativo.data <- diag(1, 300, 300)
data[1:10, 1:10]class(data)str(data)
# Con esa matriz generas 1000 listas de una lista de 300 números que tienen en algún lugar unos, uno o no los tienen, y que en total son 1048 (para el seed 123)
set.seed(123)a=replicate(1000,sample(data,size=300),simplify=F) class(a)class(a[1])class(a[1][[1]])sum(a[1][[1]])sum(a[7][[1]])sum(a[33][[1]])
sum(a[375][[1]])sum(a[999][[1]])abis=replicate(1000,sample(data,size=300),simplify=TRUE) colSums(abis)sum(colSums(abis))str(a)
# No entiendo para que haces esta segunda parte.
df<- lapply(lapply(lapply(a,t),colSums),data.frame)class(df)class(df[1])class(df[1][[1]])sum(df[1][[1]])sum(df[7][[1]])sum(df[33][[1]])sum(df[375][[1]])
sum(df[999][[1]])
# Hasta aquí llegué a la primer parte de tu planteo y me parece que no entiendo lo que dices o no entiendo tu objetivo ya que no me parece que estes obteniendo lo que declaras. 

# Para la segunda parte, tomando prestada la función de Jorge Velez, me parece que quieres hacer para la cantidad de unos por réplica i de 1 a 300, lo siguiente:
jv <- f(i,300,1000)                                           # no creo que por hacer un bucle 1:300 se pueda colgar, aunque no se si eso buscas.
colSums(jv)
Lamento pero puedo ayudarte con el motivo por el cual se te "cuelga" o cierra R con tu función, no la conozco, pero apostaría a que no es un problema de R.

Daniel Merino


El 14 de enero de 2013 17:33, Gabriel Trujillo Paucar <gabriel_enter@hotmail.com> escribió:

lapply(lapply(lapply(a,t),colSums),data.frame)



-- 
Daniel 		 	   		  
	[[alternative HTML version deleted]]

