search for: myresp

Displaying 5 results from an estimated 5 matches for "myresp".

Did you mean: myres
2008 Jul 01
1
Help in using PCR
Hi, Currently I have a dataset of 2400*408. And I would like to apply PCR method to study the any correlation between the tests. My current data is in data.frame and I have formed horizontal(1-407) to be the exact data, and (408) to be my results data(Yes and No) I have also binarized these Yes and No to 1 and -1s. However, when I refer to PCR manual on R, the example of yarn.pcr <-
2017 Dec 07
2
parallel computing with foreach()
...", "psme") covset <- c("PEN", "Thorn") foreach(i = 1:length(splist)) %:% foreach(j = 1:length(covset)) %dopar% { spname <- splist[i]; spname myTorP <- covset[j]; myTorP DataSpecies = data.frame(prsabs = rep(1, 10), lon = rep(30, 10), lat = rep(80, 10)) myResp = as.numeric(DataSpecies[,1]) myRespXY = DataSpecies[, c("lon", "lat")] # directory of a bunch of raster files specific to each R environment rastdir <- paste0(rootdir, "Current/", myTorP); rastdir rasterc = list.files(rastdir, pattern="\\.tif$", full.name...
2013 Mar 15
0
Poisson and negbin gamm in mgcv - overdispersion and theta
Dear R users, I am trying to use "gamm" from package "mgcv" to model results from a mesocosm experiment. My model is of type M1 <- gamm(Resp ~ s(Day, k=8) + s(Day, by=C, k=8) + Flow + offset(LogVol), data=MyResp, correlation = corAR1(form= ~ Day|Mesocosm), family=poisson(link=log)) where the response variable is counts, offset by the log of sample volume. Unfortunately, the residuals from the model show heteroscedasticity. While trying to follow up on this, I have run...
2017 Dec 07
0
parallel computing with foreach()
...PEN", "Thorn") > > foreach(i = 1:length(splist)) %:% > foreach(j = 1:length(covset)) %dopar% { > > spname <- splist[i]; spname > myTorP <- covset[j]; myTorP > > DataSpecies = data.frame(prsabs = rep(1, 10), lon = rep(30, 10), lat = > rep(80, 10)) > myResp = as.numeric(DataSpecies[,1]) > myRespXY = DataSpecies[, c("lon", "lat")] > # directory of a bunch of raster files specific to each R environment > rastdir <- paste0(rootdir, "Current/", myTorP); rastdir > rasterc = list.files(rastdir, pattern="\\.t...
2002 Jan 16
1
getting the response (dependent) from formula?
I'm trying to build a vector (fitting systems of equations) and I'm having a little trouble getting the response variables from a list of formulas... bdgmodel <- lbdg ~ d1sqr + ld1 + lhg hgmodel <- lhg ~ ht2 + lht + whc inst <- ~ d1sqr + ld1 + ht2 + lht + whc systemeq <- list( bdgmodel, hgmodel ) # manually generate the y matrix... y <- cbind( matrix( lbdg ),