search for: estima

Hi all, I want to make the following: I want to run a linear regression on each column of a matrix "estima" on the correspondent column on the matrix "estima2". You see I want to regress estima[,1] on estima2[,1] this way to all columns.... At the same time I want to make a regression adding each time a new observation. You see, the first regression will regress only one observation wi...

...ror but the histogram plot doesn´t appear. Should I install some special library to run sapply? pru<-function(){ randz<-matrix(rnorm(200000),100,2000) H<-matrix(0,100,2000) for (j in 2:2000){ for (i in 2:100){ H[1,]<-randz[1,] H[i,j]<-H[i-1,j]+randz[i,j] }} hy<-nrow(H)-1 estima<-H[2:nrow(H),] estima2<-H[1:hy,] a<-estima b<-estima2 mycoef <- function(x, y) coefficients( lm(y ~ x-1) ) rest <- sapply(2:2000, function(i){ y <- a[,i] x <- b[,i] mycoef(x,y) } ) print(summary(...

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

...ne for each panel - the dots with Urb=1 filled in black; and the dots with Urb=0 filled in white I've been back and forwards, trying to use panel function, but I didn't manage to do it right... Thank you for your help in advance. Regards, C. ------------------------------------ Cristina Estima Ramalho PhD Candidate MSc in GIS Ecosystem Restoration Lab School of Plant Biology (M090) The University of Western Australia 35 Stirling Highway, Crawley WA 6009 Perth, Australia Phone: +61 (08) 6488 4655; Fax: +61 (08) 6488 7461 www.plants.uwa.edu.au/research/ecosystem_restoration --------------...

I have been asked to analyse the results of (what is to me) a very complicated experiment. The dependent measure is the estimated distance, which is measured as a function of the actual distance. There are also several other IVs. The plot of log estimated distance as a function of log distance is linear. So in the rest of the analysis I will use logestimate and logdistance. My plan is to see how the other IVs affect the...

Dispongo de un experimento en el que cinco tratamientos ha sido aplicados a cinco grupos de voluntarios. En cada grupo había tres personas y a cada persona se le tomaron 3 medidas, En total dispongo de 45 medidas, pero evidentemente no son independientes entre sí. Si no tomo en cuenta que las medidas de la misma persona son más parecidas entre sí (bloques anidados) estaría incurriendo en

Dear List Has anyone used R to distnguish between alternative forecasting models? In particular is the Diebold Mariano test available for use within R. Any assistance would be greatly appreciated. Graham Leask Lecturer in Strategy Economics & Strategy Group Aston University Aston Triangle Birmingham B4 7ET Tel: 0121 204 3150 E Mail: g.leask at aston.ac.uk

Dear All, I'd like to fit a "kind" of logistic model to small data-set using nonlinear least-squares regression. A transcript of R-script are reproduced below. Estimated B and T (the model's coeff, herein B=-8,50 and T=5,46) seem appropriate (at least visually) but are quite diff from those obtained w/ SPSS (Levenberg-Marquardt): B=-19,56 and T=2,37. Am I doing something wrong in R (or at least non-comparable "methodologies")? Please, feel free to...

Hola buenas a todxs, comento mi duda: Calculo la correlación entre dos variables (calculo por años que: mayor tasa de desempleo implica menos venta de coches) y me sale positiva la correlación por lo tanto tienen una relación directa. Mi duda: ¿para que esta prediccíón sea cierta, además de una correlación positiva tengo que calcular el nivel de significancia entre las dos variables? Saludos y

Buenas. Abajo pongo la salida de un modelo de cox , dónde he estratificado por una variable de país (Countryb) y por otra (Q6). Además hay interacción entre la variable mobilityPDurG2 (es una variable 0,1, y 0 es la categoría de referencia) país. La categoría de referencia para país es "united kingdom". Mi duda surge si quiero calcular el hazard ratio para los que tienen un 1

Buen dia!! Pues me encuentro trabajando con un conjunto de datos simulados que se ajustan a un modelo Ar(1) y pues queria saber si existe un comando para realizarle una prueba de homocedasticidad, pues la prueba que hay en el R es la bartlett.test pero pues no estoy muy seguro para usarla. Gracias [[alternative HTML version deleted]]

...elo = paste(modelo, " + x1[,", i, "] + x2[,", i, "]", sep="") } res = lm (as.formula(modelo))   Este modelo deriva en resultados correctos de su análisis correspondiente (summary(res)). Sin embargo, si añado regresión polinómica, el summary no me devuelve los estimadores correctos para las componentes lineal, cuadrática y cúbica.   La línea del modelo la modifico del siguiente modo:   modelo = "y ~ I(t) + I(t^2) + I(t^3) + x1[, 1] + x2[, 1]"   He probado también así:   modelo = "y ~  t +  t^2 + t^3 + x1[, 1] + x2[, 1]"   ¡Pero no hay manera...

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 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

Hola, ¿qué tal? Tengo una pregunta de esta