similar to: prcomp compared to SPAD

Hello again, I was hoping one of you could help me with this problem. Consider the sample data from R: > summary(prcomp(USArrests)) Importance of components: PC1 PC2 PC3 PC4 Standard deviation 83.732 14.2124 6.4894 2.48279 Proportion of Variance 0.966 0.0278 0.0058 0.00085 Cumulative Proportion 0.966 0.9933 0.9991 1.00000 How do I access the

hello everybody, I'm new here. so, I'm running on SimplyMEPIS 7.0 64-bit, and using a 64-bit Wine 1.0 RC1. I'm having a problem that I had in the past too. certain programs give an error about not being able to access a certain file. I copied the software onto the Linux partition in its whole, data files and all, even ran Wine as root, but to no avail. I have the impression this is

Hello. I have a data base with 50 qualitative variables and a lot of individuals. I try to estimate the links between one of these variables (landcover) and the 49 others (geomorphology, hydrography...). I want to use a "discriminant analysis on qualitative variables" (as DISQUAL in SPAD) or a " log-linear model ". Which R-Package(s) or other methods can you advise me.

Hello all, I was wondering if anyone knew how to get R to only spit out a certain portion of PRcomp data; namely, when I enter the following code, I get: > Summary(prcomp(USArrests)) Importance of components: PC1 PC2 PC3 PC4 Standard Deviation 83.732 14.212 6.489 2.483 Proportion of Variance 0.966 0.0278 0.0058 0.00085 Cumulative Proportion 0.966

Dear R-Helpers, I have the following code: data = read.csv("test.csv",header=T,na.strings="-",row.names=1) myData=as.matrix(data) ##myData is a 102x19 matrix myList = alply(myData, 2,function(x){matrix(x,ncol=6,byrow=T)}) ##myList is a list containing 19 matrices with the following names > names(myList) [1] "*ES" "*ET" "*ER"

Hello I would like to transfer some Splus functions in R. But I have a problem first about this assignation in Splus : xnom <- deparse(substitute(x)) I am a bad programmer : I don't understand the R help How to modify these functions ? Thank you very much for your help Here are the four functions and a data test

Hi all, I've become quite enamored of R lately, and have decided to try to teach some of its basics (reading in data, manipulation and classical stats analyses) to my fellow grad students at the University of Toronto. I sent out a mass email and have already received some positive responses. One student, however, wanted to know what differentiates the routines that R uses, from those

Yep, as others have mentioned - using a linker-generated gdb-index is super helpful (10s for my machine to start debugging an unoptimized clang, set a breakpoint at llvm::StringRef::StringRef and run until that breakpoint is hit, as opposed to 1m26s without an index) That's also with/without split DWARF too, but I suspect most of the benefit is in the index there. Split DWARF might save you

Hello. Here is my R code. I used the functional data . Now I need to use the functional data by applying the kernels instead of the xi, yi functions. Bonjour. Voici mon code en R . J'ai utiliser les donn?es fonctionnelles . Maintenant j'ai besoin d'utiliser les donn?es fonctionnelles en appliquant les noyaux ? la place des fontions xi, yi library(MASS)

First of all, I thanks you all for replying me . I've not asked my question in good terms, so I've made a little drawing of what I want to do in my boxplot : - ) | - - | | | | ---- ---- | | | | | | | 2 |

Hi, I would like to know if it is possible to make global variables in R. My problem is that I use the same data frame in many different functions that I have written, but it implies a very big waste of time because it slow down the program. So I would use my data frame without need to copy it each time in another function to have a faster program execution. How can I do ?? Thanks for Help

Hello, I would like to do a partial least square discriminant analysis (PLSDA) in R using the package "ropls" Which is in R available via the R command : source("https://bioconductor.org/biocLite.R") I try to do a PLSDA to illustrate the impact of two genders (AP,C) on 5 compounds measured in persons (samples) should be illustrated. When I try to do a PLSDA I get the warning

Dear List I am trying to fit a multinomial model using the mlogit package. Attempting to load the data into mlogit presents the following error. MLOG<-mlogit.data(Mult3,shape="long",choice="CHOICE",alt.var="mode.ids",indivs = "set3",chid.var = "obs") Error in `row.names<-.data.frame`(`*tmp*`, value = c("1.1", "1.2",

Dear all, I've used the 'prcomp' command to calculate the eigenvalues and eigenvectors of a matrix(gg). Using the command 'principal' from the 'psych' package  I've performed the same exercise. I got the same eigenvalues but different eigenvectors. Is there any reason for that difference? Below are the steps I've followed: 1. PRCOMP #defining the matrix

Hello, I would like to perform a Partial least square discriminate analysis (PLSDA) in R. To do this I use the package mixOmics. I could perform the PLSDA in R. however I would also like to perform a leave-one-out cross validation in order to assess the performance of my model. My supervisor told me that I should focus on the R2/Q2 ratios. However when I read the instruction for running the

Hello, I would like to do a partial least square discriminant analysis (PLSDA) in R using the package "ropls" Which is in R available via the R command : source("https://bioconductor.org/biocLite.R") When I try to do a PLSDA using my own data. The impact of two genders (AP,C) on 5 compounds measured in persons (samples) should be illustrated. When I try to do a PLSDA I get

Hi all, I can't find the error in the binomial GLM I have done. I want to use that because there are more than one explanatory variables (all categorical) and a binary response variable. This is how my data set looks like: > str(data) 'data.frame': 1004 obs. of 5 variables: $ site : int 0 0 0 0 0 0 0 0 0 0 ... $ sex : Factor w/ 2 levels "0","1": NA NA NA

Hi, Can anyone tel me how to plot two boxplot at the same abscisse value ? Because I want to compare the discrimination of two different biological methods of the same data. Thanks Christine -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

Dear all, Last autumn there was some discussion on the list of the warning Warning message: fitted probabilities numerically 0 or 1 occurred in: (if (is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y, when fitting binomial GLMs with many 0 and few 1. Parts of replies: "You should be able to tell which coefficients are infinite -- the coefficients and their standard errors will

Hi, I know that this question has been asked many times but i can't find the correct solution from anywhere. I have 4 computer on the LAN connected as follow. comp1(WinXP)---------------------------------> hub1---->router comp2(winXP)---------------------------------> hub1---->router comp3(WinXP)wireless ---> access point----> hub1--->router comp4(CentOS