search for: aleatoria

Buenas tardes amigos, De nuevo por aqui con un incoveniente, tengo el siguiente arreglo: > MuestraS [1] 1 0 0 0 1 0 1 1 1 1 1 0 Deseo realizar todas las permutaciones posibles para luego tomar una muestra aleatoria pequeña, esto lo debo hacer varias veces incrementando el largo del arreglo "MuestraS". El inconveniente esta en que al hacer las permutaciones con este arreglo de 12 espacios ya tengo problemas al parecer de memoria y posiblemente el arreglo deba ir hasta 100. Al ejecutar el siguiente c...

Hola a todos He estado trabajando con el paquete osmar y funciona muy bien.Puede alguien decirme donde encontrar informacion o algun script para generar una lista con N pares de puntos (longitud/latitud) aleatorios o siguiendo una funcion determinada dentro del marco de coordenadas definido por corner_bbox(left, bottom, right, top) para despues plotearlos. saludos eloy ortiz -- Este

Cordial saludo. Estoy trabajando con una densidad conjunta algo elaborada (suma y cociente de dos variables aleatorias positivas usando la cópula Frank) y requiero encontrar la distribución marginal de cada una de las variables. A manera de ejemplo (no es el caso específico, que tiene mucho más código), si partimos de la distribución exponencial bivariada #----- theta<-0.5 f<-function(x,y){exp(x+y+theta*x*y)...

...ratando de correr un survival analysis usando un Cox regression model. Tengo una duda respecto a la organizacion del script. Tengo una variable que es -tamano del individuo- y quiero ver si hay diferencia en sobrevivencia respecto a tamano. Como diseno de campo los tamanos fueron ubicados de forma aleatoria en bloques al azar. Cuado planteo el script tengo algo como: modelo<-coxp(surv(deaths,day)~size+strata(block)) El factor bloque es random pero no se si esta es la nomenclatura correcta para testar mi hipotesis, o simplemente debo omitir el termino -strata- Gracias,saludos Eric Flo...

Buena tarde, Estoy llevando a cabo un trabajo y no encuentro la forma de llenar una matriz con el p_value y un estadístico calculado. Un poco mas detallado, tengo muestras aleatorias, calculo por ejemplo la prueba t, wilcoxon, etc y requiero llevar por una parte los p_value de k muestras a una matriz creada anteriormente y por otro lado el valor del estadístico de las mismas k muestras a otra matriz creada anteriormente. Las k muetras estarían determinadas por un ciclo. Agrad...

...el. O eso me parece ¿Por qué? pues porque la probabilidad de sacar 8 "ceros" u 8 "unos" seguidas es 0.003 y a mi me está pareciendo que lo hace mucho más a menudo de lo debido. Y ora cosa, ¿sabéis si cuando apago el ordenador me aparecen al iniciarlo de nuevo las mismas series aleatorias, en el mismo orden que en el trabajo previo? y de ser así ¿existe alguna función en R que las haga iniciar por otro lugar de la serie interna predefinida? ¿Sabéis si hay test de aleatoriedad para series predefinidos en R? Es que el buscador de funciones ?? me busca en los paquetes incluidos pero...

hola! podrian ayudarme con esto: Genere 1000 muestras aleatorias de tamaño n=10 de una población normal con media =2,5 y desviación 0,16. Con cada una de las muestras construya un intervalo de confianza al 95% para la media y para la varianza. Cuente que porcentaje de los 1000 intervalos atrapan la verdadera media y la verdadera varianza poblacional respecti...

...inomial M(n; N/n), es decir: rmultinom(1, n, N/n) o si deseas generar ?b? remuestras bootstrap: rmultinom(b, n, N/n) Fíjate que el proceso anterior es equivalente aunque más compacto y posiblemente más rápido al típico proceso que se suele designar como ?bootstrap no paramétrico?, generar muestras aleatorias y con reemplazamiento de la muestra original. Las frecuencias anteriores serían asimilables a una gran muestra, un vector de ?n? valores 1, 2, ?, 6, en la que hubiese 12 valores 1, ningún valor 2, 8 valores 3, etc. Tomando una muestra aleatoria (sample, etc.) de tamaño n del vector anterior y lueg...

Buenos dias para todos, Estoy trabajando en una aplicación que involucra análisis de clusters. Básicamente el objetivo es determinar a qué cluster pertenece cada observación de una matriz de datos "mydata" y luego generar muestras aleatorias de los mismos datos para determinar la proporción de veces que cada observación es clasificada en el cluster k. Este tipo de análisis es muy común en aplicaciones de Phylogeny < http://en.wikipedia.org/wiki/Phylogeny> Suponiendo que el método de cluster a utilizar es fijo y que de alguna m...

...l: 1. Tengo en el vector Muestras: rm(list = ls()) set.seed(2085089)# Fija Datos data<-matrix(-1,30,ncol=10) > MuestraS [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 1 0 0 0 1 1 0 1 1 1 Ahora, dentro de un bucle for voy a generar muestras aleatorias, pero requiero que en ese mismo proceso, apenas se tenga una muestra duplicada esta se sobrescriba, no he logrado que el bucle se quede en una posición y de esta manera generar una nueva muestra que no este duplicada. Para hace lo descrito en el parrafo anterior uso el siguiente código: a<-1...

...ertinentemente, el codigo que corro es el siguiente, la idea es que a medida que crece el tamaña muestral pues la potencia de la prueba aumente, pero no se como hacer para que, repita la prueba variando el tamaño de la muestra pues no logro acomodar el elemento del vector para que genere la muestra aleatoria de datos normales. Mx=170 M0=168 alpha=0.05 sigma=8 T=qt((1-alpha/2),n-1) T2=qt(alpha/2,n-1) Nr=1000 Nrechazo=0 n=seq(10,900, 70) dF <- data.frame(PromX = rep(0,length(n)), VarX = rep(0,length(n))) for(i in 1:Nr){ x <- rnorm(n[i], Mx, sigma) dF$PromX[i] <- mean(x) dF$VarX[i] <- var...

...babilidad puntual y <- c(1, 1.3, 1.5, 2.1, 2.8) k <- length(y) # Frecuencia o número de veces de cada observación y_j s <- c(1, 1, 2, 3, 1) n <- sum(s) # Número de observaciones # Conjunto de riesgo. Número de observaciones mayores o iguales a cada valor # que toma la variable aleatoria (a cada y_j) r <- sort(cumsum(s), decreasing = TRUE) # Kernel gamma # Probabilidades discretas a suavizar p <- s/n # Parámetro de forma común de las distribuciones gamma mi2 <- sum(y**2*s) mi4 <- sum(y**4*s) alfa <- sqrt(n/(mi4/mi2 - 1)) # Función de distribución su...

...ct),data = efi_sem, method = "REML") lme.efi<-lme(efi ~year* plot,random=~1|year/week | plot/transect,data = efi_sem, method = "REML") Mi variable respuesta será el número de grillos, como variables fijas he puesto parcela y año (solo hay dos años de muestreo), como variable aleatoria he puesto el transecto (cada transecto se corresponde con una única parcela). El problema es que no sé cómo definir el muestreo semanal como medida repetida dentro de cada año, ¿podríais ayudarme? Muchas Gracias de antemano por vuestra atención y colaboración, Lucía [[alternative HTML version d...

Hello I have problem with wrr and more than 1300 classes when i try to add qdisc like: tc qdisc add dev imq1 parent 1:21 handle 1000: wrr dest ip 1512 0 then i have a message like: RTNETLINK answers: Cannot allocate memory And some other problem: When I do more than 2048 filters I have message like: Classes: tc class add dev imq1 parent 1:1 classid 1:9 htb rate 1024kbit ceil 2048kbit .. tc

Queridos R-users: Tengo una duda que hace mucho tiempo que estoy intentando resolver, os explico a modo de ejemplo: Tengo estos efectos: Año(5 niveles),Localidad (10 niveles) y genotipo (3 niveles), año y localidad son aleatorios y genotipo es fijo (los he escogido yo). Me gustaría hacer obtener una tabla parecida a la Tabla Anova donde aparezca cada factor y sus interacciones y

Buena tarde amigos, En días pasados hice algunas consultas y ya pude salir de las dudas que tenia en ese momento, ahora requiero de su colaboración con lo siguiente: Tengo un vector dicotomico (Binario) con la siguiente información que me surgio de algunos procesos anteriores: > MuestraS [1] 1 1 1 1 0 1 1 0 1 1 Ahora necesito hallar todas las posibles combinaciones que se puedan hacer con

Hola a todos. Estoy utilizando la función lmer del paquete lme4 para ajustar un modelo mixto. Tengo varias variables en mi data.frame, unas son a nivel individual y otras a nivel de comarcas. Listo algunas. ingre_6 : Ingresos (nivel individual) iscs_a : un indicador sintético resumen de otras variables, calculado mediante componentes principales. sau_com : superficie agraria útil de cada una

...supervivencia a 3 a?os de pacientes con c?ncer colorrectal en funci?n de una serie de factores. De la muestra total de unos 1000 casos, cre? dos grupos: uno de entrenamiento y otro de test. El modelo lo hice con el grupo de entranamiento. La asignaci?n de los pacientes a cada grupo lo hice de forma aleatoria a trav?s de la librer?a CaTools, en una proporci?n 70/30. Bueno, he conseguido hacer el nomograma a partir del modelo basado en la regresi?n de Cox con la librer?a rms, pero ahora no s? c?mo calibrarlo y validarlo. He intentado utilizar "calibrate" del paquete "rms", pero me p...

...supervivencia a 3 a?os de pacientes con c?ncer colorrectal en funci?n de una serie de factores. De la muestra total de unos 1000 casos, cre? dos grupos: uno de entrenamiento y otro de test. El modelo lo hice con el grupo de entranamiento. La asignaci?n de los pacientes a cada grupo lo hice de forma aleatoria a trav?s de la librer?a CaTools, en una proporci?n 70/30. Bueno, he conseguido hacer el nomograma a partir del modelo basado en la regresi?n de Cox con la librer?a rms, pero ahora no s? c?mo calibrarlo y validarlo. He intentado utilizar "calibrate" del paquete "rms", pero me p...

...estos tuneles. Explico mi situaci?n, tengo un servidor centos 5.5, con sistema ssh que viene por defecto la versi?n 4.3. A este servidor se conectan alrededor de 100 terminales, cada uno de estos terminales levantan 3 tuneles contra el servidor, puerto 80, 3306, 1752, pero constantemente y en forma aleatoria estos tuneles se cierran, por lo que debo volver a crear los tuneles. He aplicado pol?ticas tanto en el servidor como en los clientes de keeplive, serveraliveinterval, serveralivecountmax, clientaliveinterval, clientalivecountmax. He modificado par?metros a nivel de capa 4, protocolo TCP, para mant...