search for: loo

Dear R user, I am working with LOO. Can any one who is working with leave one out cross validation (LOO) could send me the code? Thanks in advance Alex [[alternative HTML version deleted]]

I have two ideas about it. 1- i) Entering variables in quadratic form is done with the command I (variable ^ 2) - plsr (octane ~ NIR + I (nir ^ 2), ncomp = 10, data = gasTrain, validation = "LOO" You could also use a new variable NIR_sq <- (NIR) ^ 2 ii) To insert a square variable, use syntax I (x ^ 2) - it is very important to insert I before the parentheses. iii) If you want to make the interaction between x and x ^ 2 use the command ":" -> x: I(x ^ 2) iv) For mul...

...at 3:07 AM, Luigi Biagini <luigi.biagini at gmail.com> wrote: > I have two ideas about it. > > 1- > i) Entering variables in quadratic form is done with the command I > (variable ^ 2) - > plsr (octane ~ NIR + I (nir ^ 2), ncomp = 10, data = gasTrain, validation = > "LOO" > You could also use a new variable NIR_sq <- (NIR) ^ 2 > > ii) To insert a square variable, use syntax I (x ^ 2) - it is very > important to insert I before the parentheses. True, but better I believe: see ?poly. e.g. poly(cbind(x1,x2,x3), degree = 2, raw = TRUE) is a full qu...

...t;luigi.biagini at gmail.com> wrote: >> I have two ideas about it. >> >> 1- >> i) Entering variables in quadratic form is done with the command I >> (variable ^ 2) - >> plsr (octane ~ NIR + I (nir ^ 2), ncomp = 10, data = gasTrain, validation = >> "LOO" >> You could also use a new variable NIR_sq <- (NIR) ^ 2 >> >> ii) To insert a square variable, use syntax I (x ^ 2) - it is very >> important to insert I before the parentheses. > > True, but better I believe: see ?poly. > e.g. poly(cbind(x1,x2,x3), degr...

I'm reposting this problem (perhaps a bug) now I've got more information on it. This is another point of view of the situation and I hope someone could have run into the same trouble before (and solved it :-)) This is it: * with ntbackup 2000 I create a 22Gb .bkf file in the windows machine. * I can copy that file over a samba share and get correct info form the file in windows

...f) and therefore cannot be estimated. Have I totally misunderstood what you meant or committed some other blunder? Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Sun, Jul 16, 2017 at 7:36 AM, David Winsemius <dwinsemius at comcast.net> wrote: > >> On Jul 13, 2017, at 7:43 AM, Bert Gunter <bgunter.4567 at gmail.com> wrote: >> >> Below. >> >> -- Bert >> Bert Gunter >> &g...

Hi Folks, I am wanting advise on a good soft-phone on Linux. I have looked at Gnophone but cannot seem to get it to compile under debian sarge. I am now looing at sipXphone seem to be picking up that it is not that stable, but perhaps someone here can advise on what softphone I can use on Linux. Thanks in advance, Hamish ---------------------------------------------...

Sorry about the one-man thread but if someone would slap me... floor0_forward is declared "int". It returns "val". "val" is declared "long". "val" is assigned the return from "rint" which is "double". Am I misreading this? ber --- >8 ---- List archives: http://www.xiph.org/archi...

Hello, I am trying to compile a file .Rd into .html using R CMD Rdconv -t=html /file.rd/>/file.html./ It seems that the process works ok but instead of having for example loo.cv(NMRTools) at the top of the html file I have loo.cv{unknown}. Any ideas on how to solve this? Best, Dani // -- Daniel Valverde Saub? Grup de Biologia Molecular de Llevats Facultat de Veterin?ria de la Universitat Aut?noma de Barcelona Edifici V, Campus UAB 08193 Cerdanyola del Vall?s- SPA...

Hi All, I have been trying to get this LOO-Cross Validation method to work on R for the past 3 weeks but have had no luck. I am hoping someone would kindly help me. Essentially I want to code the LOO-Cross Validation for the 'Local Constant' and 'Local Linear' constant estimators. I want to find optimal h, bandwidth. Than...

While doing pls I found the following problem > BHPLS1 <- plsr(GroupingList ~ PCIList, ncomp = 10, data = PLSdata, jackknife = >FALSE, validation = "LOO") when not enabling jackknife the command works fine, but when trying to enable jackknife i get the following error. >BHPLS1 <- plsr(GroupingList ~ PCIList, ncomp = 10, data = PLSdata, jackknife = >TRUE, validation = "LOO") Error: cannot allocate vector of size 289.1...

I have been using samba 2.2.2 for quit awhile running with security= share in my company. I looked into the /var/lock/samba and found a growing file unexpected.tdb. Is there a way to trim it down while samba is running or do I have to start and stop samba? Is this an expected behavior? TIA Joseph Loo

...in 1:100 ) { pp <- round(x[i] * pRatio) if( pp == x[i] ) { pp <- pp-1 } if( pp >= Nmin ) { x[i+1] <- pp } else { break } } x } ## R-SVM core code ## input: ## x: row matrix of data ## y: class label: 1 / -1 for 2 classes ## CVtype: ## integer: N fold CV ## "LOO": leave-one-out CV ## "bootstrape": bootstrape CV ## CVnum: number of CVs ## LOO: defined as sample size ## Nfold and bootstrape: user defined, default as sample size ## output: a named list ## Error: a vector of CV error on each level ## SelFreq: a matrix for the frequency of...

...nt apparently *must* be explicitly named if NIR is not a numeric vector. AFAICS, this is unclear or unstated in ?poly. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Thu, Jul 13, 2017 at 10:15 AM, David Winsemius <dwinsemius at comcast.net> wrote: > >> On Jul 12, 2017, at 6:58 PM, Ng, Kelvin Sai-cheong <kscng at connect.hku.hk> wrote: >> >> Dear all, >> >> I am using the pls package o...

...unction with interaction terms. But I am not sure how. I will use an example to illustrate my problem: Following the example in the PLS manual: ## Read data data(gasoline) gasTrain <- gasoline[1:50,] ## Perform PLS gas1 <- plsr(octane ~ NIR, ncomp = 10, data = gasTrain, validation = "LOO") where octane ~ NIR is the model that this example is fitting with. NIR is a collective of variables, i.e. NIR spectra consists of 401 diffuse reflectance measurements from 900 to 1700 nm. Instead of fitting with predict.octane[i] = a[0] * NIR[0,i] + a[1] * NIR[1,i] + ... I want to fit the...

...ust* be explicitly named if NIR is > not a numeric vector. AFAICS, this is unclear or unstated in ?poly. I still get the same error with: library(pld) data(gasoline) gasTrain <- gasoline[1:50,] gas1 <- plsr(octane ~ poly(as.matrix(NIR), 2), ncomp = 10, data = gasTrain, validation = "LOO") Error in rep.int(rep.int(seq_len(nx), rep.int(rep.fac, nx)), orep) : invalid 'times' value > gas1 <- plsr(octane ~ poly(as.matrix(gasTrain$NIR), degree=2), ncomp = 10, data = gasTrain, validation = "CV") Error in rep.int(rep.int(seq_len(nx), rep.int(rep.fac, nx)...

...... I am now stuck while conducting the plsr. I have tried various methods of creating structured lists etc and have got nowhere. I have also tried many incarnations of BHPLS1 <- plsr(GroupingList ~ PCIList, ncomp = FeaturePresenceExpected[1], data = FullDataListTrans, validation = "LOO") Where am I going wrong. Also what is the easiest method to identify which of the 94,000 peptides are most important to the variance between groups. Thanks in advance for any help Amit Patel [[alternative HTML version deleted]]

...am not sure > how. I will use an example to illustrate my problem: > > Following the example in the PLS manual: > ## Read data > data(gasoline) > gasTrain <- gasoline[1:50,] > ## Perform PLS > gas1 <- plsr(octane ~ NIR, ncomp = 10, data = gasTrain, validation = "LOO") > > where octane ~ NIR is the model that this example is fitting with. > > NIR is a collective of variables, i.e. NIR spectra consists of 401 diffuse > reflectance measurements from 900 to 1700 nm. > > Instead of fitting with predict.octane[i] = a[0] * NIR[0,i] + a[1]...

I just want to set up a folder on my Linux 7.2 running samba 2.2.1a-4. I have looked at the smb.conf over and over. Cannot find out why my share is not showing up. I just want users on my domain to browse the network and read the folder. I am on an NT domain and the furthest i got was seeing the workgroup testgrp on my network neighborhood. But when i click on the workgroup...

I'm running RH7.1 and Samba 2.2.2 . I would like to offer a VFAT partition on the server hard drive for universal read-write access. I mount the partition in fstab with /dev.hda5 /sys vfat defaults,rw,uid=1001,gid=1001 0 0 It mounts fine and shows up on the Windows 98 machine and can be opened. I can't write to it however. User and group 1001 are smbuser account, and they