similar to: necesito ayuda para realizar contrastes

Hola a todos! perdón por molestarlos nuevamente con los contrastes, pero estoy trantando de entender que es lo que está haciendo R y de donde vienen los valores que informa pero  no lo logro. Creí haberlo entendido pero a la hora de usar mis datos los resultados no dan como deberían. Tengo dos variables explicativas que son factores con 3 niveles cada uno. Esta es la tabla de medias de la

Estimados usuaRios: deseo hacer una consulta, debo realizar una prueba no paramétrica a un conjunto de datos que posee una variable de respuesta con tres factores que poseen interacción, o al menos debo probar ello. Es decir, realizar una prueba Scheirer-Ray-Hare o de Contrastes de rangos alineados. Estuve trabajando el paquete ARTool aunque no me queda claro

Buenas tardes. Tengo un ensayo de dos factores (Variedad y dosis) Quisiera hacer anovas para cada una de las dos variedades que estoy evaluando. Lo puedo hacer cogiendo los datos primero de una variedad y luego de la otra. Pero creo que hay un comando "by variedad" que lo hace si necesidad de andar modificando el archivo de datos. Pero no sé dónde hay que colocarlo. Si alguno me puede

Estimados usuarios: Estoy intentando reproducir el ejemplo 6.4 de Thomas Lumley. Complex Survey. Editorial Wiley. 2010 (ver la página en google:

estimados, ajusté una superficie de respuesta con el paquete "rsm" de R y hago un grafico 2D de la superficie de respuesta con contour(), pero como el argumento de esta función es la superficie de respuesta ajustada y ésta se ha ajustado codificando las variables, de acuerdo a la recomendación de la documentación del paquete rsm, obtengo un gráfico con los ejes en unidades codificadas

Hola compañeros de la lista, qué tal. Los molesto con la siguiente duda: Tengo un experimento con dos factores A y B, cada uno de los cuales tiene los siguientes niveles (que son concentraciones de dos hormonas vegetales aplicadas a plantas): niveles del factor A: 0, 0.2, 0.5, 1 niveles del factor B: 0, 0.1, 0.2, 0.5, 1 y mi variable de respuesta es continua, todo dentro del set de datos

Hola, tenemos como equivalencia no paramétrica al Anova de 1 vía el test Kruskal Wallis, pero ¿existe un test no paramétrico equivalente al Anova de 2 vías? Gracias. Saludos. -- Sec. Apoyo Estadístico. Servicio de Apoyo a la Investigación (SAI) Vicerrectorado de Investigación e Internacionalización. Universidad de Murcia Edif. SACE. Campus de Espinardo. 30100 Murcia @. mariaelvira.ferre en

I don't understand the coefficients returned from the lm function. I expected these to be the mean values for each factor in the model. Given this data and model: data<-c(rnorm(10,mean=0,sd=1),rnorm(10,mean=1,sd=1),rnorm(10,mean=-.5,sd=1)) ftr<-as.factor(rep(1:3,each=10)) fit<-lm(data ~ ftr) the mean values of the three facotrs from the data: c(mean(data[1:10]),

Hola Ruben, Sí precisamente es lo que comentas, en matemáticas no se suele llamar bucketización (este término se emplea más en informática) sino datos agrupados. Pero la idea es la que tu mismo dices. Respecto a las gráficas que has puesto, me han aclarado mucho sobre el tema, gracias. Si realizo lo mismo, por ejemplo con nbucket=1000 sigo obteniendo un p-valor de 1. Es decir, que casi le

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

Hi, I am a very occasional user of R, and will be grateful for some help in constructing a regression across groups. Here is an example: library(MASS) attach(cats) Sex[120:144]<-factor(TG) #Renaming some males to transgender, to create 3 groups, male, female and transgender out<-lm(Bwt~Sex/Hwt) #Gives me 3 separate linear regressions for groups M, F and TG What I now want to do

Hola. Algún paquete que calcule la potencia de un U-test? -- Antonio Maurandi [[alternative HTML version deleted]]

Hello, Does anyone know of a way I can speed this up? Basically I'm attempting to get the data item on the same row as the report date for each report date available. In reality, I have over 11k of columns, not just A, B, C, D and I have to do that over 100 times. My solution is slow, but it works. The loop is slow because of merge. # create sample data z.dates =

hi Michael, am sorry for the incomplete reply . csv file data havinmg like this , >getSymbols("IBM") >weekly_data = to.weekly(week_name) >dataframe=data.frame(weekly_data) >outputfile_name="F:\\R-programs\\Outputfile.csv" >write.table(dataframe, file =outputfile_name,sep = ",",col.names =TRUE,row.names = T) > datafrom_table <-

Hola a tod en s, queremos hacer una comparación entre dos lugares muy alejados entre sí en relación a la temperatura de cada sitio usando medias horarias de un período determinado. Sólo hay medidas de un sensor en cada sitio y queremos saber si las diferencias son significativas o no entre sitios/curvas. Hemos usado un test de Mann?Whitney U con la función wilcox.test (paired=F) ya que los

Hi All, I've been having a little trouble with creating a loop that will run a a series of t.tests for inspection, Below is the code i've tried, and some checks i've looked at. I've used the get(paste()) idea as i was told previously that the use of the eval should try and be avoided. I've run a single syntax to check that my systax is correct and works without any problems