similar to: Calculate SPE in PLS package

Greetings, I have been exploring the use of the caret package to conduct some plsda modeling. Previously, I have come across methods that result in a R2 and Q2 for the model. Using the 'iris' data set, I wanted to see if I could accomplish this with the caret package. I use the following code: library(caret) data(iris) #needed to convert to numeric in order to do regression #I

All, I have 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

Hi Jean-Marc, We needed to use a release version of the SDK because this part of production code for our Sony implementations. I'm not sure how stable that branch is. Is the encoding and decoding already stable? In other words, what is not stable about that branch? Thanks, --------------------------- Saad Nader Middleware Engineer Powered By Gamespy, IGN -----Original Message----- From:

Hello, I am relatively new to R and I am working through the code that is provided in the book Numerical Ecology with R and I have run across an error message that I can't seem to figure out. I am using the vegan, ade4, gclus and cluster packages. The code is as follows: # Ordered community table # Species are ordered by their weighted averages on site scores or <- vegemite(spe,

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

Dear developeRs, I maintain a package 'pls', which has a main fit function mvr(), and functions plsr() and pcr() which are meant to take the same arguments as mvr() and do exactly the same, but have different default values for the 'method' argument. The three functions are all exported from the name space. In the 'pre namespace' era, I took inspiration from lm() and

Hi all, >From a list of strings, I desire to filter out the followings: 1. Digits at the beginning of the strings 2. Character "SPE" following the digits (if it exists) 3. Any characters followed by hyphen The following produces the desired result, but would like to know whether this can be done more efficiently. Any suggestions would be much appreciated. dat <- c("2148

Dear all: I am comparing the PLS outputs of R and SAS for the following data set: Y x1 x2 x3 3 6 2 2 3 1 5 5 4 7 4 1 5 6 5 6 2 4 3 2 8 5 0 9 where Y is the dependent variable and x1, x2, x3 are the independent variables. I found several PLS algorithms in R (NIPALS,SIMPLS,KERNEL PLS). SAS has SIMPLS and NIPALS. The following are the NIPALS calculations of

Hi Jean-Marc, Jim, Saad has been keeping me in the loop on your recent discussions. Since all of our testing has been against 1.0.5, based on that being the last non-beta version, that's the particular scope of the task that Saad is working on right now. I like what I'm reading about as far as encoder/decoder quality improvements e.g. in the 1.2 betas, and am going to push for

Please correct me if I am wrong, Jean-Marc, but I do not think that any patches are getting applied to 1.0.5 anyway. Also, if you expect a patch to be applied, you will need to provide the changes as a patch, not as a modified copy of the source tree. The 1.2 branch includes a mechanism for private memory allocation from a static buffer. You provide a usermisc.h file that replaces the

Joost Schuur wrote: > My primary concern is the use of the word 'unstable' on the current > download page for 1.2b2. One of our major devloper partners in > particular saw that reference and opted to use 1.0.5 on their PS3 title, > which is why we based our work for them on 1.0.5. The kind of commercial > game developers that are our customers aren't going to have the

Dear all, I have run the confint() function according to below and I get the following error message: > confint(stepAIC.glm.spe.var.konn2.abund, level=0.95) Waiting for profiling to be done... Error: missing value where logical needed In addition: Warning message: NaNs produced in: sqrt((fm$deviance - OriginalDeviance)/DispersionParameter) or > confint(stepAIC.glm.spe.var.konn2.abund,

Version 2.0-0 of the pls package is now available on CRAN. The pls package implements partial least squares regression (PLSR) and principal component regression (PCR). Features of the package include - Several plsr algorithms: orthogonal scores, kernel pls and simpls - Flexible cross-validation - A formula interface, with traditional methods like predict, coef, plot and summary - Functions

Version 2.0-0 of the pls package is now available on CRAN. The pls package implements partial least squares regression (PLSR) and principal component regression (PCR). Features of the package include - Several plsr algorithms: orthogonal scores, kernel pls and simpls - Flexible cross-validation - A formula interface, with traditional methods like predict, coef, plot and summary - Functions

Saad Nader wrote: > We needed to use a release version of the SDK because this part of > production code for our Sony implementations. I'm not sure how stable > that branch is. Is the encoding and decoding already stable? In other > words, what is not stable about that branch? It's listed as "unstable" in the sense of "I might change the API" and that is

Skipped content of type multipart/alternative-------------- next part -------------- A non-text attachment was scrubbed... Name: speex-1.0.5.zip Type: application/x-zip-compressed Size: 641376 bytes Desc: speex-1.0.5.zip Url : http://lists.xiph.org/pipermail/speex-dev/attachments/20071022/2447cf50/speex-1.0.5-0001.bin

Dear R-gurus, Just supose I have a dara.frame that looks like myDF<-read.table(stdin(),head=T,sep=",") codID,namesp,k1,k2,k3,k4 1,spA,2,5,6,3 2,spB,4,5,4,6 3,spC,2,1,5,6 4,spC,5,4,3,2 5,spD,1,2,3,4 6,spE,2,4,3,1 I need to update the columns k1-k4 with the namesp, but considering the math between Kx and codID. My desired output must looks like: codID,namesp,k1,k2,k3,k4

Hello, I've got an ISDN PRI circuit terminating in a Cisco AS5350, which in turn is talking to an Asterisk server via SIP for call origination and termination. Seems simple enough, and it works for the most part, but: 1) Caller ID name data comes in on the PRI, but doesn't appear to get handed off to the Asterisk server via SIP, at least not in any format that Asterisk

Joost Schuur wrote: > Saad has been keeping me in the loop on your recent discussions. > Since all of our testing has been against 1.0.5, based on that being the > last non-beta version, that's the particular scope of the task that Saad > is working on right now. I'm just saying it's useless to do any work on 1.0.5. It is inferior to the latest beta in all aspects,

Dear all, Is there any function for calculating confidence limits for coefficients in an lm() object? I know of the confint() function in the MASS library working very well on my binomial GLMs and I have tried it (using glm () , family=gaussian) but it gives NAs according to below. Does the confint() function not accept gaussian GLMs? Could there be convergence problems in the GLM? Note the