similar to: nonparametric covariance analysis

Dear all, For my research I need to perform a Levene's test on my data.I have not large experience with this test and with R,so I have some problems to use the leveneTest function. I read some previous posts on this topic but I cannot understand what I have to do. I have a five-columns data frame organised in the following way: number of observation, number of subject, score, type of

Hi, i would to know, if someone have ever write the code to estimate the parameter (mixing proportion, mean, a var/cov matrix) of a mixture of two multivariate normal distribution. I wrote it and it works (it could find mean and mixing proportion, if I fix the var/cov matrix), while if I fix anything, it doesn't work. My suspect is that when the algorithm iterates the var/cov matrix, something

Hi, i would to know, if someone have ever write the code to estimate the parameter (mixing proportion, mean, a var/cov matrix) of a mixture of two multivariate normal distribution. I wrote it and it works (it could find mean and mixing proportion, if I fix the var/cov matrix), while if I fix anything, it doesn't work. My suspect is that when the algorithm iterates the var/cov matrix, something

Hi. I have tried to use Package "SEM". As a learning, I try to convert a program running well of EQS which is as follows to SEM: ### EQS ### /SPECIFICATION CAS=100; VAR=5 MAT=COR; ANA=COR; /EQUATIONS V1=*F1+E1; V2=*F1+E2; V3=*F1+*F2+E3; V4=**F1+*F2*E4; V5=*F2+E5; /VAR E1 TO E5=*; F1*1.0; F2=1.0; /COV E1,E2=*; F1,F2=*: /PRINT FIT ALL; /MATRIX ...... /END This is the converted SEM

Dear list, I don't know what's the correct term for this in the OOP context, but is it possible to "clone"/copy an instance of a reference class (say 'a') so that I get an *autonomous* second instance 'b'? Autonomous in the sense that changes to 'a' do not affect 'b'. I know that this is somewhat against the pass-by-reference paradigm, but the

hello, I try to run simulation of levene's test to find the p-value but the error of replacement has length zero occur, could anyone help me to fix this problem? asim <- 1000 pv<-rep(NA,asim) for(i in 1:asim) {print(i) set.seed(i) g1 <- rnorm(20,0,2) g2 <- rnorm(20,0,2) g3 <- rnorm(20,0,2) x <- c(g1,g2,g3) group<-as.factor(c(rep(1,20),rep(2,20),rep(3,20))) library(Rcmdr)

Dear Mariano, See the function leveneTest() in the car package. I hope that this helps, John ----------------------------- John Fox, Professor Emeritus McMaster University Hamilton, Ontario, Canada Web: socialsciences.mcmaster.ca/jfox/ > -----Original Message----- > From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Marcelo > Mariano Silva > Sent: Monday, January

Hola, lo pude corregir de la siguiente forma: Edité el archivo Rprofile.site agregando las siguientes líneas: options(download.file.method = "libcurl") Sys.setenv("http_proxy"="http://mi_proxy::mi_puerto") Sys.setenv("https_proxy"="http:// http://mi_proxy::mi_puerto") Sys.setenv("ftp_proxy"="http://

Hi R-friends, Attached is the SAS XPORT file that I have imported into R using following code library(foreign) mydata<-read.xport("C:\\ctf.xpt") print(mydata) I am trying to maximize logL in order to find Maximum Likelihood Estimate (MLE) of 5 parameters (alpha1, beta1, alpha2, beta2, p) using NLM function in R as follows. # Defining Log likelihood - In the function it is noted as

Sebastián, tengo dificultades para entender si sigues con el ejemplo anterior o esto es algo nuevo. Tú descripción se debe entender que dentro de aggregate(cbind(X1,X2) ~ B + C, t, function(x) mean(x < 2)) utilizas un table? Aparentemente ese no es el problema sino que el resultado que obtienes es una lista, ¿te sirve ?unlist ? De cualquier manera no estoy frente a tu ordenador y solo estoy

Daniel, en vez de mean utilicé table. Entonces me genera una lista y quiero que quede una matriz. Unlist no me sirve porque me pasa todos los valores a un vector y me pierdo de saber a que caso corresponde. Saludos, Sebastián. El 26 de septiembre de 2013 16:50, daniel <daniel319@gmail.com> escribió: > Sebastián, tengo dificultades para entender si sigues con el ejemplo > anterior

Intento cargar el paquete reshapeGUI y me da el siguiente error: > library(reshapeGUI) Loading required package: reshape2 Loading required package: plyr Loading required package: gWidgets Attaching package: ‘gWidgets’ The following object(s) are masked from ‘package:plyr’:     id Loading required package: gWidgetsRGtk2 > reshapeGUI() Error en eval(expr, envir, enclos) : objeto

I'm trying to treat a continuous variable as discrete for plotting multiple regression lines in a scatterplot as a function of the level on the moderating variable. In the example below, there is only one regression line plotted to the whole data. I would like a separate regression line for each discrete level of the moderator. The moderator is continuous, so I'd like to treat it as

Gracias Javier. Voy a probar. Saludos, Sebastián. Enviado desde Correo para Windows 10 De: Javier Marcuzzi Enviado: domingo, 10 de marzo de 2019 14:53 Para: Sebastian Kruk CC: Lista R Asunto: Re: [R-es] Problema codificación archivos Estimado Sebastian Kruk Yo utilizo Microsoft R-open 3.5.1 y la r-project se encuentra en 3.5.3, ¿si intenta con cualquiera de estas dos? Javier Rubén Marcuzzi

Estimados amigos, en el siguiente link dejo una direccion para la descarga de Runner una interfaz gratuita para R (sólo para windows xp o posterior ) que puede ser últil para alguno de Uds. https://sites.google.com/site/runnerstatisticalsoftware/download Posiblemente mirar el manual podria ayudarlos a evaluar si esta interfaz puede ayudarlos en el trabajo diario. Cualquier sugerencia será

Sebastián, Una forma es aggregate(cbind(X1,X2) ~ B + C, t, function(x) mean(x < mean(x))) Saludos, Jorge.- 2013/9/25 Sebastian Kruk <residuo.solow@gmail.com> > Hola Juan. > > ¿Si yo en vez de querer obtener la media quiero sacar el porcentaje de > casos que son menores que la media como lo haría? > > Saludos, > > Sebastián. > > > El 17 de septiembre

Sebastián No se para que vas a usar los datos, pero yo no usaría ceros cuando el dato no esta disponible ( me saltearía el último paso de la rutina que te mando). Seguramente hay mejores y mas eficientes maneras de hacerlo en R pero esta es una quizás no muy elegante: original <- matrix( c(1 ,1 ,1, 1, 223, 2, 1, 224, 3, 3, 447, 1, 78, 335.75, 8, 2901, 186.3, 2, 2901, 886.18, 3, 2901,

Dear Sirs, My name is Marcelo Grimaldi and I work for GE Lighting Brazil in Six Sigma Program. In order to accomplish a statistical project I need to generate exact tables for the Kruskal-Wallis statistic ("Three or More Populations Location Test") and I wonder if R could help me. In the past I have heard of generating tables for Wilcoxon statistic using R. Could you help me ? Thank

The new package "benchden" 1.0.0 implements 28 benchmark densities for nonparametric density estimation that were introduced by A. Berlinet and L. Devroye ("A Comparison of Kernel Density Estimates", Pub. Inst. Stat. Univ. Paris, XXXVIII, fasc. 3, 1994, 3-59, http://cg.scs.carleton.ca/~luc/devs.html ). This collection includes a variety of densities with different degrees of

The new package "benchden" 1.0.0 implements 28 benchmark densities for nonparametric density estimation that were introduced by A. Berlinet and L. Devroye ("A Comparison of Kernel Density Estimates", Pub. Inst. Stat. Univ. Paris, XXXVIII, fasc. 3, 1994, 3-59, http://cg.scs.carleton.ca/~luc/devs.html ). This collection includes a variety of densities with different degrees of