similar to: Simulacion matrices de varianza-covarianza

Hola, ¿qué tal? Tu problema es que tienes agrupaciones de países (p.e., ASIA_MERIDIOL) con un solo respresentante. Las matrices con las que construyes las elipses son degeneradas (más bien, no existen: sospecho que tienen que ver con la covarianza... ¿de un único caso?). Elimina las agrupaciones triviales. Un saludo, Carlos J. Gil Bellosta http://www.datanalytics.com El día 4 de diciembre de

Buenas tardes, Estoy interesado en generar observaciones de una distribucion binomial bivariada en la que hay _cierto_ grado de correlacion (denotemoslo rho). Podria por favor alguien indicarme como hacerlo en R? Este es el contexto. Supongamos que se tienen dos experimentos en los que la variable respuesta _sigue_ una distribucion binomial, i.e., X_i ~Binomial(n_i, p_i), i=1,2 y que, por ahora,

Cordial saludo adjunto la base de datos y el script: eu60<- read.csv(file.choose(), header=T, sep=";", dec=".", row.names=1) eu60.pca <- PCA(eu60, quali.sup=19) eu60.data <- cbind.data.frame(eu60[,19], eu60.pca$ind$coord) eu60.ellipse <- coord.ellipse(eu60.data, bary=TRUE) plot.PCA(ellipse=eu60.ellipse, cex=0.8) El 4 de diciembre de 2013 17:35, Camilo Calle

Estimados Escribo para consultar sobre el uso de modelos mixtos anidados. Los datos que estoy analizando provienen de censos de malezas en cuatro tipos de paisajes de la región pampeana, en los que seleccioné al azar igual número de lotes agrícolas cultivados con tres cultivos (maíz, soja y trigo-soja). En cada lote censé el número de especies de malezas en tres posiciones: el alambrado, el borde

Cordial saludo, soy principiante en R y estoy haciendo un análisis multivariado de datos, y necesito hacer las elipses de confianza donde la variable 19 es cualitativa suplementaria y el resto son variables activas; tengo el siguiente código: eu60 <- read.csv(file.choose(), header=T, sep=";", dec=",", row.names=1) ### capturar datos eu60.pca <- PCA(eu60, quali.sup=19) ##

Alguien tiene experiencia en escribir una pdMat. Para aquellos que no lo recuerden son las matrices de covarianzas de los efectos aleatorios que ajusta la función lme de la librería nlme Estas matrices tiene especial importancia en aplicaciones de genética de poblaciones y en particular en mapeo de asociación. Pinheiro y Bates dicen que el usuario puede crear sus propias pdMat y sugiere como

Hola, LLevo buscando desde hace tiempo como hacer el Goodness of fit test en R. Es decir, me explico, intento hacer una cosa parecida que se hace en Minitab, por ejemplo, yo tengo un conjunto de datos, y lo que quiero es sabes que tipo de distibución es, en minitab se hace un histograma para ver si se ajusta bien o no a la campana de Gauss, luego vemos si aproximar la distribución de la muestra

Hola! Estoy teniendo un problema para testar la normalidad y la homogeneidad de varianzas para mi modelo de manova Es este mi modelo: m1<-manova(cbind(pred,fit,oniv)~tratamento/parcela+bloco,data.vegetacao) Aparte de tener tres variables respuestas (predador, fitófago, omnívoro), también es un diseño anidado. Alguien sabría responder esto? Gracias! Carol [[alternative HTML version

Hi all, im triyng to implement a bayesian model with R and c++. I have a strange problem. I can't reproduce the error with a small script and then i post the original one. The problem is after the line for(MCMC_iter2=0;MCMC_iter2<thin;MCMC_iter2++) For the first 34 iterations all work fine, after, all the simulations of mu_acc_P return an "nan". If i delete the line

Muchas gracias por la rapidez y por la ayuda tan efectiva, me quedó el gráfico de maravilla, al parecer el problema era que habían individuos en una sola categoría, así quedó la gráfica: [image: Imágenes integradas 1] El 4 de diciembre de 2013 11:04, Carlos J. Gil Bellosta < cgb en datanalytics.com> escribió: > Hola, ¿qué tal? > > Tu problema es que tienes agrupaciones de

Estimados Al leer este interesante script lo corrí, pero al hacer la elipse me devuelve esa información, no encuentro el porque saludos library(FactoMineR) eu60<- read.csv(file.choose(), header=T, sep=";", dec=".", row.names=1) eu60.pca <- PCA(eu60, quali.sup=19) eu60.data <- cbind.data.frame(eu60[,19], eu60.pca$ind$coord) eu60.ellipse <-

¿Viene implementado en algún paquete el clásico test para la varianza de una normal (basado en la distribución chi-cuadrado)? Sí, ya se que es muy sencillo hacerlo uno mismo... Gracias. Saludos a tod@s y felices navidades. H. Gómez [[alternative HTML version deleted]]

Saludos. Bueno quiero realizar el analisis multivariado especificamente el analisis de correspondencia en R , cual es el paquete implementado? [[alternative HTML version deleted]]

??Gracias Javier,? 2017-04-01 12:07 GMT-03:00 <javier.ruben.marcuzzi en gmail.com>: > Mi duda es la siguiente, la varianza residual en su modelo, es homogénea o > heterogénea, ?La varianza es homogénea, común a todas las observaciones. Digamos que el modelo es el siguiente, y=f(x, betas)+e con f alguna función no lineal cuyos parámetros son betas y con e~N(0, sigma2). Uno suele

Estimados, tengo que analizar un experimento que fue ejecutado bajo el diseño strip strip plot, lo que es bastante infrecuente al menos en mi area, y no quisiera hacerlo a mano. Me pregunto si en R existira algun paquete que ejecute el analisis de varianza de este tipo de experimentos. Recuerdo que un par de veces se han dado aqui paginas para consultar acerca de paquetes, funciones y esas

Hola Alguien me podría decir como hacer una grafica del tipo persp() de la densidad de una distribucion normal bivariante estandarzada con correlacion 0.5?? gracias [[alternative HTML version deleted]]

Hola Javier, gracias por contestar tan rapido (tarde en la noche) ... olvide decir que, obviamente, he revisado la internet buscando y no encuentro ni siquiera la tabla de analisis de varianza del diseño ... habia revisado el paquete agricolae porque este trae diseños parecidos pero no exactamente lo que necesito. Sigo buscando, muchas gracias, Eric. On 29/04/15 22:57,

Saludos a todos como analizo las curvas o funciones (ANALISIS DE COMPONENTES PRINCIPALES O CORRELACION CANONICA) que estan adjuntadas en el archivo JPEG, Necesito determinar el patron o una funcion equivalente, gracias a todos CARLOS ANDRES ------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL:

Muy buenas tardes, les escribo con el fin de saber si alguno de ustedes, me puede dar una guia de como generar una simulacion para ilustrar el insesgamiento de los estimadores B0 y B1 de regresion lineal simple lo he logrado hacer para un tamaño de muestra fijo pero no he podido hacer las iterciones. En particular lo que requiero es hacer la simulacion que me grafique para varios diferentes

Hola amigos, Supongamos que se quiere ejecutar un modelo no lineal con nls. Pensemos en el ejemplo de la ayuda: DNase1 <- subset(DNase, Run == 1) fm1DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal), DNase1) summary(fm1DNase1) Aquí se está modelando la densidad óptica de un ensayo relacionada de forma no lineal (logística) con (el logaritmo) de la concentración de una proteína.