similar to: Does R have function/package works similar to SAS's 'PROC REG'?

Hi dear R-users I try to reproduce the steps included in a LDA. Concerning the eigenvectors there is a difference to SPSS. In my textbook (Bortz) it says, that the matrix with the eigenvectors V usually are not normalized to the length of 1, but in the way that the following holds (SPSS does the same thing): t(Vstar)%*%Derror%*%Vstar = I where Vstar are the normalized eigenvectors. Derror

Dear list, how can I compute the inverse of the X'WX matrix ("inverse of the weighted sum of squares and crossproducts matrix") from an object of class "lm" from a weigthed linear regression? Thanks, Sven

Can R do a repeated measures MANOVA and tell what dimensionality the statistical variance occupies? I have been using MATLAB and SPSS to do my statistics. MATLAB can do ANOVAs and MANOVAs. When it performs a MANOVA, it returns a parameter d that estimates the dimensionality in which the means lie. It also returns a vector of p-values, where each p_n tests the null hypothesis that the mean

Hi all, I'm trying to convers some Cisco SSCP phones to the SIP formware. The phone boots, I see it tries to fetch a bunch of files on my TFTP: Jun 17 09:37:45 firewall dnsmasq-dhcp[21202]: DHCPACK(eth2) 192.168.10.103 6c:50:4d:da:f0:67 SEP6C504DDAF067 Jun 17 09:38:10 firewall in.tftpd[22666]: RRQ from 192.168.10.103 filename CTLSEP6C504DDAF067.tlv Jun 17 09:38:10 firewall in.tftpd[22666]:

Dear Skinny/SCCP lovers :-) I've just completed & uploaded to the cvs the newest version with fixed redial key AND implementation of speed dials. please test extensively and report any bugs. i know that the display is not yet set correctly but the buttons are working as expected. Enjoy testing... --jan (*1) http://chan-sscp.sf.net (*2) yes, bugtracker is down at the moment, will fix

Version 0.9-99 of package "subselect" is now on CRAN. The subselect package has functions that search for k-variable subsets of a p-variable (p>k) data set that are optimal under some criterion. Search algorithms include a full search algorithm "leaps", based on Furnival and Wilson's leaps and bounds algorithm, and three random search algorithms: a genetic algorithm

Version 0.9-99 of package "subselect" is now on CRAN. The subselect package has functions that search for k-variable subsets of a p-variable (p>k) data set that are optimal under some criterion. Search algorithms include a full search algorithm "leaps", based on Furnival and Wilson's leaps and bounds algorithm, and three random search algorithms: a genetic algorithm

A new version (0.7.1) of package 'subselect' has been uploaded to CRAN. Package 'subselect' provides functions which assess the quality of variable subsets as surrogates for a full data set, in an exploratory data analysis, and search for subsets which are optimal under various criteria. As of version 0.7 a new function 'leaps' has been added. 'Leaps' performs a

A new version (0.7.1) of package 'subselect' has been uploaded to CRAN. Package 'subselect' provides functions which assess the quality of variable subsets as surrogates for a full data set, in an exploratory data analysis, and search for subsets which are optimal under various criteria. As of version 0.7 a new function 'leaps' has been added. 'Leaps' performs a

Hi, I have a huge number of categorical variables, say at least 10000, and I put them into a matrix, each column is one variable. The question is: how can I make all of the pairwise cross tabulation tables efficiently? The straightforward solution is to use for-loops by looping two indexes on the table() function, but it was just too slow. Is there a more efficient way to do that? Any guidance

dear all, i'm trying to reproduce an spss-anova in R. It is an 2x3x3 repeated measures desingn with repeated contrasts. In R i've coded a contrast matrix for all factors and made a split in the aov summary - but I can't get the repeated interaction contrasts. The output from SPSS looks like this: TaskSw * CongNow * CongBefore: SS df Mean Square F Sig. 1 vs. 2 1 vs. 2 1 vs. 2

Hi, I am not sure if this is a bug and I apologize if it is something I didn't read carefully in the R extension manual. My initial search on the R help and R devel list archive didn't find useful information. I am using .Call (as written in the R extension manual) for the C code and have found that the .Call didn't release the memory claimed by allocVector. Even after applying

Hi, I am not sure if this is a bug and I apologize if it is something I didn't read carefully in the R extension manual. My initial search on the R help and R devel list archive didn't find useful information. I am using .Call (as written in the R extension manual) for the C code and have found that the .Call didn't release the memory claimed by allocVector. Even after applying

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 I'm trying to understand some C code in an R package I'm using. I'm address this question here as it's matrix algebra...and I'm no pro at that! the C command reads: double alpha = 1.0, beta = 0.0; dsyrk_("L", "N", nGenes, nGenes, & alpha, mat1, nGenes, & beta, mat2, nGenes); - From google,

Dear list: I have tried MASS's mca function and SAS's PROC corresp on the farms data (included in MASS, also used as mca's example), the results are different: R: mca(farms)$rs: 1 2 1 0.059296637 0.0455871427 2 0.043077902 -0.0354728795 3 0.059834286 0.0730485572 4 0.059834286 0.0730485572 5 0.012900181 -0.0503121890 6

Hi all, I'm just wondering if there is a equivalent of SAS's FIRST. and LAST. variables in R? For example, suppose this is a snapshot of the data: ClientCode CaseCode open close Important 1 37 28 2003-07-08 2003-09-02 1 2 37 310 2003-11-01 2004-09-10 1 3 37 1562 2007-04-03 2007-07-27 1 4

Hi try this function ive written it should be self explantory but let me know if you have any problems. I've only been using R for a few eeeks so apologies if its not the most efficient! rankit2<-function(rankvar,cuts,data,factor) { ranker<-rankvar ranker<-0 range<-c(1:cuts) range2<-range/cuts range3<-quantile(factor,range2) over<-length(factor) for (i in 1:over){

Hello- I wonder if there's a way to replicate SAS rank procedure where it ranks a variable by a certain number of groups. For example, it's very easy to calculate quintile rank in SAS, but I couldn't find the similar function in R. Does anyone know how to do this? Thanks. -- View this message in context:

Hi Ruser I'm trying to replicate some SAS code. I have to add a spline to my longitudinal spaghetti plot. I have the plot, but I can't add the spline, a overall trend line. In the SAS code they use the command 'I=SM50S' and I would prefer something similar. I?m using R 2.10.1 on windows XP? I have made this working example. tolerance.pp <-

If anyone could tell me what this error is all about, I would be very grateful. Oct 12 01:42:53 WARNING[2724]: app.c:1109 ast_lock_path: Failed to lock path '/var/spool/asterisk/voicemail/default/5933/INBOX': Operation not permitted Oct 12 01:42:53 WARNING[2724]: app.c:1109 ast_lock_path: Failed to lock path '/var/spool/asterisk/voicemail/default/5933/Old': Operation not