similar to: from .csv file to a pca plot

Does anyone have nice quality controlls for affy arrays,.... Can't find any tools as are being used for 2 dye arrays..... cheers, marinus This e-mail and its contents are subject to the DISCLAIMER at http://www.tno.nl/disclaimer/email.html [[alternative HTML version deleted]]

I am attempting to develop a multiple regression model using selected model variables that should all be treated as numeric (mostly real) values. However, R considers one specific variable "mass" automatically to be of class "factor", probably because "mass" consists of integer values that are repeated. I now want to force R to treat "mass" as a numeric

I am hoping to get some advise on the following: I am looking for an automatic variable selection procedure to reduce the number of potential predictor variables (~ 50) in a multiple regression model. I would be interested to use the forward stepwise regression using the partial F test. I have looked into possible R-functions but could not find this particular approach. There is a function

Hi to all, very strange! I tried to debug why imap process takes 100% CPu and it is very slow when I move an email from an imap folder to another imap folder and then I typed: strace -tt -o /log-imap.txt -p <process number> well, during the trace the imap process doesn't crash and it is very fast. If I stop the tracing, I have again the problem and sometimes the crash of the

[was sent a wrong R-help address; manually resent by MM] Hello I'am using the 'biplot' and 'biplot.pincomp' functions of the 'mva' package for my studies. The biplot represents both the observations and the variables of a matrix of multivariate data on the same plot. The observations are represented by their numbers (the line of the data matrix), but I would need to

Dear R Helper, I am writing because I seek to perform a varimax rotation on my Principal Components Analysis (PCA) solution. (I have been performing PCA's using the eigen command in R.) If you can tell me how to perform this rotation when I use the eigen command (or the princomp command) I would be thrilled. Thanks so much! Wendy Treynor Ann Arbor, MI USA

Hi all, I have relatively small dataset on which I would like to perform a PCA. I am interested about a package that would also combine a method for determining the number of components (I know there are plenty of approaches to this problem). Any suggestions about a package/function? thanks, Nick -- View this message in context:

Dear expe-R-ts, I try some test data for a factorAnalysis (resp. pca) in the sense of Prof. Ripley's MASS ? 11.1, p. 330 ff., just to prepare myself for an analysis of my own empirical data using R (instead of SPSS). 1. the data. ## The test data is (from the book of Backhaus et al.: Multivariate ## Analysemethoden. Springer 2000 [9th ed.], p. 300 ff):

Dear All, I want to know if there is some easy and reliable way to estimate the number of dominant eigenvalues when applying PCA on sample covariance matrix. Assume x-axis is the number of eigenvalues (1, 2, ....,n), and y-axis is the corresponding eigenvalues (a1,a2,..., an) arranged in desceding order. So this x-y plot will be a decreasing curve. Someone mentioned using the elbow (knee)

Hello I am able to extract partial regression coefficients from a fitted model object "model", i.e. model <- lm(var.sel.gkm, weights = count.gkm, data = DATA) summary(model) write.table(model$coef, file = "C:/coef_CO_gkm.txt", row.names = TRUE, col.names = TRUE) I was wondering if anyone could advise me how to extract other object items such as std. error, t-values

Hello all, I have performed PCA on stacked wind vector data for 20 different spatial locations. In this case I therefore am in effect working out a PCA for 40 different stations (2 paired pieces of information [ie u and v vector information] for each station), so have 40 PCs. In my previous PCA (for none vector data) i have plotted the PC loadings onto a map of the area of concern, and would

Dear All, I am analysing a dataset on levels of herbivory in seedlings in an experimental setup in a rainforest. I have seven classes/categories of seedling damage/herbivory that I want to analyse, modelling each separately. There are twenty maternal trees, with eight groups of seedlings around each. Each tree has a TreeID, which I use as the random effect (blocking factor). There are two

Hello, The estimate of glm dispersion can be based on the deviance or on the Pearson statistic. I have compared output from R glm() to another statastical package and it appears that R uses the Pearson statistic. I was wondering if it is possible to make use R the deviance instead by modifying the glm(...) function? Thanks for your attention. Kind regards, Robin Smit This e-mail and its

On Wed, Feb 11, 2015 at 06:41:20PM +0800, Chen Hanxiao wrote: > + /* Output of `btrfs scrub -R status' is like: > + * > + * scrub status for 346121d1-1847-40f8-9b7b-2bf3d539c68f > + * scrub started at Mon Feb 2 17:39:38 2015, running for 93 seconds > + * data_extents_scrubbed: 136670 > + * tree_extents_scrubbed: 30023 > + *

Hello I am trying to use all subsets regression on a test dataset consisting of 11 trails and 46 potential predictor variables. I would like to use Mallow's Cp as a selection criterion. The leaps function would provide the required output but does not work with this many variables (see below). The alternative function regsubsets should be used, but I am not able to define the function in

While setting my first couple IP phones, I set their voicemail buttons to an extension that runs VoicemailMain. exten => 8500,1,Wait(1) ; voicemail exten => 8500,2,VoicemailMain ; exten => 8500,3,Hangup ; I would like to be able to pass the mailbox number allowing each phone to go in directly but I'd rather tno have

Dear list members, I apologize for putting this (probably) very basic question on the mailing list. I have scanned through the R website (using search) but did not found an answer. (code included below) A factor matrix is simply extracted (which can then subsequently be exported using write.table) by FACT$loadings[1:6,]. I would also like to specifically extract and export

Peoplefone AG jest firm? oferuj?c? us?ugi telekomunikacyjne oparte na technologii Voice over IP (VoIP) w wyj?tkowo niskich stawkach. Peoplefone jest certyfikowanym partnerem firm Siemens<http://www.siemens.ch/index.jsp?sdc_p=c175fi1012637lmno1012637psuz1&sdc_sid=1113876080&>i AVM/FRITZ!Box <http://www.fritz-shop.ch/>. W zwi?zku z naszym dynamicznym rozwojem, do nowo tworzonej

Hello all, I am somewhat new to S-plus and would really appreciate any help you might have on the following problem: I installed S-plus (0.65.1) on a Sun box (Solaris) and it runs fine as far as I can tell. I also installed the cluster and multiv contributed packages. The problem I have is that I get a core dump everytime I try to run a function in the contributed packages. This seems

Dear R users, I'm trying to cluster 30 gene chips using principal component analysis in package mva.prcomp. Each chip is a point with 1,000 dimensions. PCA is probably just one of several methods to cluster the 30 chips. However, I don't know how to run prcomp, and I don't know how to interpret it's output. If there are 30 data points in 1,000 dimensions each, do I have to