search for: oscorespl

All, I am working with NIR spectral data and it was great to find that the example in ?plsr also used spectral data. Unfortunately, I am having difficulty figuring out how the "yarn" dataset is structured to allow for the plsr model to read: library(pls) data(yard) yarn.oscorespls <- mvr(density ~ NIR, 6, data = yarn, validation = "CV", method = "oscorespls") dim(yarn) yarn Specifically, it is VERY convenient to be able to have the ~256 variables to be used in the model defined by NIR term. When I use dim(yarn) it claims [1] 28 3. When I call &...

..., I took inspiration from lm() and implemented it like this: plsr <- function(..., method = pls.options()$plsralg) { cl <- match.call() cl$method <- match.arg(method, c("kernelpls", "widekernelpls", "simpls", "oscorespls", "model.frame")) cl[[1]] <- as.name("mvr") res <- eval(cl, parent.frame()) ... Recently, Prof. Brian Ripley kindly pointed out that this doesn't work properly when the 'pls' package in not attached: > data(yarn, package='pls') &...

...quot;, preProc=c("scale")) data(iris) training1=iris[inTrain1,] datvars=training1[,1:4] dat.sc=scale(datvars) n=nrow(dat.sc) dat.indices=seq(1,n) timematrix=with(training1, classvec2classmat(Species[dat.indices])) pls.dat=plsr(timematrix ~ dat.sc, ncomp=3, method="oscorespls", data=training1) x=crossval(pls.dat, segments=10) summary(x) summary(plsFit2) I see two different R2 values and I cannot figure out how to get the Q2 value. Any insight as to what my errors may be would be appreciated. Regards, -- Charles [[alternative HTML version deleted]]

...ve been using a data set to build pls models for three different soil properties. Two of the three models run fine; however I receive the following error for the final model. > libs.IC.cal <- mvr(libs.IC.fmla, data = libsdata.cond.cal, ncomp=20,validation = "LOO", method = "oscorespls") Error in colMeans(x, n, prod(dn), na.rm) : 'x' must be numeric There are many "0" for this soil property. Could this cause the error? Best, Ross ******************************************************************* Ross Bricklemyer Dept. of Crop and Soil Scien...

Dear list, I want to calculate SPE (squared prediction error) in x-space, can someone help? Here are my codes: fit.pls<- plsr(Y~X,data=DAT,ncomp=3,scale=T,method='oscorespls',validation="CV",x= T) actual<-fit.pls$model$X pred<-fit.pls$scores %*% t(fit.pls$loadings) SPE.x<-rowSums((actual-pred)^2) Am I missing something here? Thanks in advance. Stella Sim DISCLAIMER:\ This email contains confidential informatio...{{dropped:11}}

...ng on the algorithm used, the size of the matrices, and the number of components used, one can expect from 5% to 65% reduction in computation time. - Scaling of scores and loadings of kernel PLS and svd PCR algorithm has changed. They are now scaled using the `classic' scaling found in oscorespls. - The arguments `ncomp' now always means "number of components", and `comps' always means "component number". The argument `cumulative' has been removed. - A new data set 'gasoline' has been included. - The 'NIR' and 'sensory' data sets...

...ng on the algorithm used, the size of the matrices, and the number of components used, one can expect from 5% to 65% reduction in computation time. - Scaling of scores and loadings of kernel PLS and svd PCR algorithm has changed. They are now scaled using the `classic' scaling found in oscorespls. - The arguments `ncomp' now always means "number of components", and `comps' always means "component number". The argument `cumulative' has been removed. - A new data set 'gasoline' has been included. - The 'NIR' and 'sensory' data sets...

...ay x1, x2, ..., x30).  Aside from studying the relationship of Y and x1, ..., x30, I am also interested in studying the effect of quadratic terms (x1^2, ..., x30^2) and two-way interactions (x1*x2, x1*x3, ..., x29*x30) on Y. My R code is dat.pls<-plsr(Y ~ X^2 + I(X^2), ncomp = 2, method="oscorespls") This code does not work for me.  I got this from Chapter 11 of "An Introduction to R." Any suggestion of what the appropritae code should be is highly appreciated. Thanks, Jaclyn M. [[alternative HTML version deleted]]