search for: sub_x

...Y should have same size !") return } if (is.null(lag.before)) { lag.before <- l_X / 10 #should do 10*log10(N/m) } if (is.null(lag.after)) { lag.after <- l_X / 10 # idem lag.before. } res <- list() res$cor <- rep(0.0, 1 + lag.before + lag.after) sub_X <- X[ 1 + lag.before : (l_X - lag.after)] res$delay <- seq(-lag.before, lag.after) for (i in res$delay) { sub_Y <- Y[(lag.before+i+1) : (l_X-lag.after+i+1)] OK <- !(is.na(sub_X) | is.na(sub_Y)) x_select <- sub_X[OK] y_select <- sub_Y[OK] res$cor[i+lag.b...

...ata$group==3),3]~data[which(data$group==3),2]) abline(reg_3.1, col="green", lty=1) #Segmented package to estimate breakpoints and determine slopes for group2 and 3: subset_data <- subset(data, subset=data$group!=3) names(subset_data) <- c("sub_group", "sub_x", "sub_y") attach(subset_data) X <- model.matrix(~0+sub_group)*sub_x time.1 <- X[,1] time.2 <- X[,2] time.3 <- X[,3] olm <- lm(sub_y~0 + sub_group + time.1 + time.2) os <- segmented(olm, seg.Z= ~ time.1 + time.2, psi=list(time.1=c(14,30)...