search for: matris

Hi I'm trying to utilize the break command for breaking the loop when the p-value is less than 10 per cent using the urca package. But it does not break the loop, anyone that can help me? library(urca) set.seed(1) a1 <- runif(100) lag.max <- function(object, n = 12){ matris <- matrix(NA, nrow = n) for(i in 1:n) { matris[i] <- ur.df(object, lags = i, type = "trend")@testreg$coefficients[i+3,4] if (i < 0.1) {break(i)} } list(matris = round(matris, 3)) } a2 <- lag.max(a1) /With regards Serda...

If you submit these lines, you end up with variable "vare.dis". I want to export vare.dis to csv. Stuck I am. library(vegan,logical.return = TRUE) #return=true verifies package is available library(MASS,logical.return=TRUE) #return=true verifies package is available data(varespec) #varespec is an example data file in the vegan package vare.dis <- vegdist(varespec)

> version _ platform sparc-sun-solaris2.8 arch sparc os solaris2.8 system sparc, solaris2.8 status major 1 minor 6.0 year 2002 month 10 day 01 language R

Dear all, I have 2 matrices, one is a 8x3 matrix, called X; the other matrix is a 3x3 indicator matrix with the diagonal element as 0 or 1. when a variable is included in the model, the corresponding diagonal element is 1, otherwise, it is 0. Let A be a set of matrices that contain the possible indicator matrix. e.g., A= [A1, A2, A3], where A1= [1,0,0;0,0,0,0,0,0], A2 =[1,0,0;0,1,0,0,0,0], A3

...ue teclear esas operaciones a cada Mapa y por ende prolongar la vida de mi precioso teclado con valor de 50 pesos mexicanos #Operación general Mapa1 <- readRAST6 ('Mapa1') ###Para importar un raster de GRASS-GIS a R Mapa1.v <- as.vector(Mapa1.m <- as.matrix(Mapa1)) ###Convierto a matris y después a vector Mapa1.vn <- Mapa.v[!is.na(Mapa.v)] ### Elimino datos nulos del vector SP1 <- sample(Mapa.vn, 50) ### Extraigo una muestra de 50 datos. El asunto ahí no termina, después de realizar esas operaciones necesito calcular las "medidas de separabilidad" mediante el pa...

I am relatively new to R so I am not confident enough in what I am doing to be certain this is a bug. I am running R 2.1.1 on a Windows XP machine and the lme4 package version 0.98-1. The following code fits the model I want using the nlme package version 3.1-60. mltloc$loc <- factor(mltloc$loc) mltloc$block <- factor(mltloc$block) mltloc$trt <- factor(mltloc$trt) Mltloc <-