search for: comp

Hello all, any help with the following would be appreciated. I am attempting to maximize 28 parameters similtaneously in Optim(). The follwoing is my code: Y<-comp[,9] n<-length(Y) e.1<-.5 e.2<-.5 e.3<-.5 e.4<-.5 e.5<-.5 e.6<-.5 e.7<-.5 e.8<-.5 e.9<-.5 e.10<-.5 e.11<-.5 e.12<-.5 e.13<-.5 e.14<-.5 e.15<-.5 e.16<-.5 e.17<-.5 e.18<-.5 e.19<-.5 e.20<-.5 e.21<-.5 e.21<-.5 e.22<-.5 e.23<-....

...:2] 0.99083 0.00852 0.93983 0.23186 0.85931 ... ..- attr(*, "dimnames")=List of 2 .. ..$ : chr [1:160] "1" "2" "3" "4" ... .. ..$ : chr [1:2] "1" "2" $ terms :Classes 'terms', 'formula' length 3 groups ~ Comp.1 + Comp.2 + Comp.3 + Comp.4 + Comp.5 + Comp.6 + Comp.7 + Comp.8 + Comp.9 + Comp.10 + Comp.11 + Comp.12 + ... .. ..- attr(*, "variables")= language list(groups, Comp.1, Comp.2, Comp.3, Comp.4, Comp.5, Comp.6, Comp.7, Comp.8, Comp.9, Comp.10, Comp.11, Comp.12, Comp.13, ......

...results to a file using a write.table() I get the following error message "Error in as.data.frame.default(x[[i]], optional = TRUE) : can't coerce loadings into a data.frame" Any idea why write.table is doing that and any possible solutions?? thanks hector > loadings(M2.princomp.corr) Loadings: Comp.1 Comp.2 Comp.3 Comp.4 Comp.5 Comp.6 Comp.7 Comp.8 Comp.9 Comp.10 f61 0.158 -0.138 -0.341 f63 0.234 -0.359 f67 0.178 -0.113 -0.329 f75 0.187 -0.199...

Hi, I just wonder if this is a rounding error by the princomp command in R. Although this does not make much sense, using a hypothetical dataset, a, a<-matrix(runif(1000),100,10) I did PCA with the princomp, and compared it with the results estimated with the eigen and the prcomp commands. And I found some differences in the results: opposite signs in...

I start a clean session of R 1.9.1 on Windows and I run the following code: > library(MASS) > data(painters) > pca.painters <- princomp(painters[ ,1:4]) > loadings(pca.painters) Loadings: Comp.1 Comp.2 Comp.3 Comp.4 Composition 0.484 -0.376 0.784 -0.101 Drawing 0.424 0.187 -0.280 -0.841 Colour -0.381 -0.845 -0.211 -0.310 Expression 0.664 -0.330 -0.513 0.432 Comp.1 Comp.2 Comp.3 Comp.4...

Hi!! I'm having trouble selecting 10 out of 41 attributes of the KDD data set. In order to identify the components with the higher variance I'm using princomp. the result i get for summary(pca1) is: Comp.1 Comp.2 Comp.3 Comp.4 Comp.5 Comp.6 Comp.7 Comp.8 Comp.9 Comp.10 Standard...

I have a question about princomp(mva) that I hope isn't too newbie. I used the sample data from Table 1.1 in "Manly (1986/1994) Multivariate Statistical Methods: a primer. Chapman and Hall" on sparrow body measurements. I rescaled the data to mean 0 and SD 1, and the covariance matrix is: V1 V2...

Full_Name: Sahotra Sarkar Version: 2.2.0 OS: Windows XP Professional Submission from: (NULL) (146.6.130.180) The following two commands should give the same results for the eigenvectors but do not (there is a sign reversal for the first one): > summary(princomp(bumpus),loading = TRUE) Importance of components: Comp.1 Comp.2 Comp.3 Comp.4 Comp.5 Standard deviation 6.2801653 2.4285636 1.13995042 0.560733747 0.03458915 Proportion of Variance 0.8399662 0.1256084 0.02767525 0.006696272 0.00002548 Cumulative Proport...

I have the manual for S+ 6 and I'm trying to use R for the Principle Component Analysis example and I'm getting a few interesting answers... The log is as follows: R : Copyright 2001, The R Development Core Team Version 1.3.0 (2001-06-22) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type `l...

Hi I''m having a issue with one object when try to update the attributes in this object previously saved in the database I have one object comp and one object reg they have this relationship comp has_many regs reg belongs_to comp when run use the method valid? like this current_contribuyente.comps.build(params[:comp]).valid? it returns: IndexError (string not matched) And when I try to update directly the comp directly it return...

...ogic EP9315 @ 200 MHz (ARM920T) o.s.: Linux 2.6.9-rc2-ep93xx GCC: gcc version 3.3.3 (DENX ELDK 3.1 3.3.3-10) without MaverickCrunch support ogg library: 1.0 speex command line configuration: ./configure --enable-fixed-point --enable-arm4-asm The code is natively built on the target (i.e. no cross compilation) I performed some tests with this simple script (the test-gab.pcm file is about 10.1 seconds): echo -------------------------------------------------------------------------- echo comp 0, bitrate 6600 time ../src/speexenc -w --bitrate 6600 --comp 0 test-gab.pcm test-gab_b6600_c1.arm.spx...

...2.1/mkl/lib/intel64/libmkl_gf_lp64.a /usr/caen/intel-12.1/mkl/lib/intel64/libmkl_sequential.a /usr/caen/intel-12.1/mkl/lib/intel64/libmkl_core.a -Wl,--end-group -lpthread" ./configure --prefix=/home/software/rhel6/R/2.15.1 \ -with-x=no --with-blas --with-lapack --enable-R-shlib Everything compiles OK, but 'make check' shows some differences. I tried this first with gcc 4.7.0 and then again with gcc 4.4.6 (RHEL 6 installed). Same results. There are differences for the package 'tools', and 'stats'. The first, perhaps, I could ignore because of a change in format...

Hi, Newbie question about the pls package. Setup: Mac OS 10.3.9 R: Aqua GUI 1.01, v 2.0.1 I want to get R^2 and Q^2 (LOO and Leave-10-Out) values for each component for my model. I was running into a few problems so I played with the example a little and the results do not match up with the comments in the help pages. $ library(pls) $ data(NIR) $ testing.plsNOCV <- plsr(y ~ X, 6, data = NIR, method="kernelpls", validation="none"...

...ginal bug report PR# 6005 on Mon, 22 Dec 2003) ------------------------ Michael Drescher Ontario Forest Research Institute Ontario Ministry of Natural Resources 1235 Queen St East Sault Ste Marie, ON, P6A 2E3 Tel: (705) 946-7406 Fax: (705) 946-2030 ------------------------ Code: > sfalls.plot.comp <- matrix(nrow=plot.count, ncol=spec.count, byrow=T) > colnames(sfalls.plot.comp) <- levels(SPECIES) ### SPECIES, SPP_VOL, & PLOT are columns/variables in a previously read data file > sfalls.plot.comp <- data.frame(sfalls.plot.comp) > attach(sfalls.plot.comp) > sfalls.plot...

...ycle(tso), tso), file=paste(file,".dat",sep=""), sep=" ", quote=F, row.names=F, col.names=F) } readx12out <- function(file, adtype = "M", calendar = F, tblnames = NULL) { notbls <- 6 + (calendar != F) + length(tblnames) comp <- vector("list", notbls) if(calendar == F) { names(comp) <- c("original", "seasonal", "adjusted", "trend", "irregular", "weights", tblnames) tblnames <- c("b1", "d10...

Hello, I tried to do a 'sem' analysis for data of how blueberry consumption by birds is influenced by a pollution gradient, using distance and vegetation structural and composition variables, but I got the following error message: Error in sem.default(ram = ram, S = S, N = N, param.names = pars, var.names = vars, : S must be a square triangular or symmetric matrix This may be very obvious for R specialist, but I cant find the problem! #Symbolic ram model mod.BB...

Dear all, I have a question about the 'sdev' value returned by the princomp function (which does principal components analysis). On the help page for princomp it says 'sdev' is 'the standard deviations of the principal components'. However, when I calculate the principal components for the USArrests data set, I don't find this to be the case: Here is...

...t the values produced are nearly of the same magnitude but of opposite sign. When I use the example data (sensory) I find this result reproduced. I am prepared to work this through but I have a feeling that there could be a possible error in the code. (?!) > sens.pcr$loadings Comp 1 Comp 2 Comp 3 Comp 4 yellow 75.186621 -0.4780473 3.212149 1.750123 green -90.490256 8.5880530 1.634961 1.042239 brown -2.861241 -11.3600509 -15.920789 -1.105799 glossy 13.347090 19.3103902 -3.121693 2.781282 transp 20.126987 24.0653312 -6.656764 -1.842907 syrup -...

...riable names rather than numbers. >str(x) loadings [1:94727, 1:10] -0.00113 -0.03001 -0.00059 -0.00734 -0.02969 ... - attr(*, "dimnames")=List of 2 ..$ : chr [1:94727] "PCIList1" "PCIList2" "PCIList3" "PCIList4" ... ..$ : chr [1:10] "Comp 1" "Comp 2" "Comp 3" "Comp 4" ... - attr(*, "explvar")= Named num [1:10] 14.57 6.62 7.59 5.91 3.26 ... ..- attr(*, "names")= chr [1:10] "Comp 1" "Comp 2" "Comp 3" "Comp 4" ... Here is the structure o...

Hi I am attempting to do a loadings plot from a plsr object. I have managed to do this using the gasoline data that comes with the pls package. However when I conduct this on my dataset i get the following error message. >plot(BHPLS1, "loadings", comps = 1:2, legendpos = "topleft", labels = "numbers", >xlab = "nm") Error in loadingplot.default(x, ...) : Could not convert variable names to numbers. str(BHPLS1_Loadings) loadings [1:8892, 1:60] -0.00717 0.00414 0.02611 0.00468 -0.00676 ... - attr(*, "...