similar to: R version of SAS Proc Varclus

Dear List What is the R package equivalent of Proc Varclus or Information Value. ANy assistance in determining R equivalents of f Oblique Component Analysis (PROC VARCLUS), Information Value (IV) and Weight Of Evidence (WOE) analysis, and business intelligence http://www.nesug.org/proceedings/nesug06/an/da23.pdf Regards, Ajay Websites- http://decisionstats.com [[alternative HTML version

Hi, I'll apreciate your guidance on how can I re-create the output from SAS PROC VARCLUS in R. I've found the varclus function in Hmisc. However, is it possible to use that function to compute for each variable the 1-R**2 ratio (this is the ratio of 1 minus the R-squared with Own Cluster to one minus the R-squared in the Next Closest cluster)? Thanks in advance for any help, Lars.

Hi All, I'd like to build a package for the community that replicates the output produced by SAS "proc varclus". According to the SAS documentation, the first few steps are: 1. Find the first two principal components. 2. Perform an orthoblique rotation (quartimax rotation) on eigenvectors. 3. Assign each variable to the rotated component with which it has the higher squared

Dear R users: is there a package, similar to varclus in SAS or varclus in S, ported or written for R? Also, is there any other package in R that was designed for grouping the variables under different measures of distance (in cases where data is non-Gaussian, autocorrelated, and so on). Janusz. -- ** Janusz Kawczak ** ** UNC at Charlotte,

Dear R users, Two questions: 1) Is there a way to simplify the mtext() line below ? beta=c(1,-1) m=5 plot(1) mtext( bquote(paste( beta == .(paste( "(", paste(beta, collapse=", "), ")" )) )), outer=TRUE,line=-3) 2) How do I get the embedded carriage return "\n" below to work, i.e for the text that follows it to appear on the next line? beta=c(1,-1) m=5

The workshop organizing committee of the Southern California Chapter of the American Statistical Association announces the 28th Annual Workshop in Applied Statistics. Professor Colin Cameron from the department of Economics at UC Davis will give a one-day workshop titled "Advances in Count Data Regression." The event will take place on Saturday, March 28, 2009 at UCLA. Details on

Hi there, I'm sure this must have been discussed before, but I can't find anything that addresses exactly my problem in the archives. The problem I ran into is, that a samba-share relying on POSIX file-permissions stops working as soon as the shared folder is a NFS-mount itself. In that case, it is no longer possible for any user to access files or folders. So in my example, if

I've searched the archives and read the xyplot help but can't figure out the 2 lattice questions below? Consider: library(lattice) DF <- data.frame(x=rnorm(20), y=rnorm(20), g1=rep(letters[1:2], 10), g2=rep(LETTERS[1:2], each=10), g3=rep(rep(letters[3:4],each=5),2)) xyplot(y ~ x | g1 + g2, groups=g3, data=DF) 1) Is there a way to get one strip per row and column

I'm getting repeated email (KMail) error messages about one apparently over-large post: -------------------------- Received: from helen.gayleard.com (localhost.localdomain [127.0.0.1]) by helen.gayleard.com (Postfix) with ESMTP id E4500294A0 for <tim at helen.gayleard.com>; Sat, 24 Jan 2015 14:00:05 +0000 (GMT) Received: (from tim at localhost) by helen.gayleard.com

I don't understand Windows permissions under Samba. I am running an smb server under CentOS-5.4, and have a Windows XP Pro machine ("harriet") as a client. I have set up shares C and D on harriet corresponding to the 2 drives C: and D: on the machine, and I have asked that files should be shared. Now when I mount harriet:C I can see the top level folders, but cannot see Documents

I don't understand Windows permissions under Samba. I am running an smb server under CentOS-5.4, and have a Windows XP Pro machine ("harriet") as a client. I have set up shares C and D on harriet corresponding to the 2 drives C: and D: on the machine, and I have asked that files should be shared. Now when I mount harriet:C I can see the top level folders, but cannot see Documents

Dear R list, I've read the function writing sections on both "An introduction to R" and "R language Definition" manuals but still don't understand why the following gives an error message: fun <- function(x, ...) x + y fun(1, y=2) I get: Error in fun(1, y = 2) : object "y" not found I'd appreciate any help in understanding this. R version 2.4.1

Dear R-Help, I am trying to compute a new variable, let's call it "target cannon orientation (tco)" based conditionally on old variables, "TargetColor," "CannonOriB," and "CannonOriR." For every case in the data set, if TargetColor is "B" then I want tco to equal the value for that case of CannonOirB, else CannonOriR. I've tried writing

I need to make a data frame out of the data that I currently have in a list. This works, but is ugly: ineffData<-rbind(ineffFilesList[[1]], ineffFilesList[[2]], ineffFilesList[[3]], ineffFilesList[[4]], ineffFilesList[[5]], ineffFilesList[[6]], ineffFilesList[[7]], ineffFilesList[[8]], ineffFilesList[[9]], ineffFilesList[[10]], ineffFilesList[[11]], ineffFilesList[[12]], ineffFilesList[[13]],

I'd be interested in any corrections or comments on the following instructions (basically for myself): We assume that BackupPC has been installed: sudo yum install BackupPC 1. BackupPC must be run by the user backuppc. Accordingly the lines User apache Group apache in /etc/httpd/conf/httpd.conf should be changed to User backuppc Group backuppc 2. The user backuppc must be able to

Dear R helpers, I have a data frame where missing values for numeric variables are coded as 999. I want to recode those as NAs. The following only partially succeeds because numeric variables are converted to character in the process: df <- data.frame(a=c(999,1,999,2), b=LETTERS[1:4]) is.na(df[2,1]) <- TRUE df a b 1 999 A 2 NA B 3 999 C 4 2 D is.numeric(df$a) [1] TRUE

Hi, I've been using R for many years and have always tried to keep my R version up to date, and when I switched from 2.7.1 to 2.7.2 I'm suddenly having trouble printing from the graphics device. I have Windows XP, and installed both versions of R from the binaries. If I start 2.7.2 and simply do: > plot(1:10) the default R Graphics Device window opens as usual. However, when I

influence.measures gives several measures of influence for each observation (Cook's Distance, etc) and actually flags observations that it determines are influential by any of the measures. Looks good! But how does it discriminate between the influential and non- influential observations by each of the measures? Like does it do a Bonferroni-corrected t on the residuals identified by

Hi. I have repeated measures data of the form where each observation is a trial, and trials are grouped by subject, and variables encode whatever level of a factor was present during that trial, and the dependent variable is response time (RT). I want to transpose the data to a form suitable for MANOVA such that there is one observation per subject and RT is recoded across many

I have a strange problem on a CentOS-5.8 machine. I can only login as root. If I try to login with one of the user's names, it hangs for a long time. I thought it hung forever, but I just found that I do login after "su tim" after 5 minutes. It seems that the problem lies in repeated messages in /var/log/messages --------------------------- May 3 12:14:13 helen su: nss_ldap: failed