similar to: Winbind and AD issues

Dear R developers, there appears to be a small problem with function cmdscale: for non-Euclidean distance matrices, using option add=FALSE (the default), cmdscale misses the smallest eigenvalue. This affects GOF statistic g.1 (See Mardia, Kent + Bibby (1979): Multivariate Analysis, eq. (14.4.7). The corresponding formula in Cox + Cox (2001): Multidimensional Scaling, 2nd ed., p 38, would

Hello I am using the following command but not able to text the values on the graph can someone please make suggestions for improvement #here is the command loc_mds <- cmdscale(dist.r, k = 7, eig = TRUE) loc_mds$eig sum(abs(loc_mds$eig[1:2]))/sum(abs(loc_mds$eig)) sum((loc_mds$eig[1:2])^2)/sum((loc_mds$eig)^2) x <-loc_mds$points[,1] y <-loc_mds$points[,2] plot(x, y,

GENERAL REFERENCE ON NONLINEAR OPTIMIZATION? What are your favorite references on nonlinear optimization? I like Bates and Watts (1988) Nonlinear Regression Analysis and Its Applications (Wiley), especially for its key insights regarding parameter effects vs. intrinsic curvature. Before I spent time and money on several of the refences cited on the help pages for "optim",

Hi, Maybe I have not been looking in the right spot, but, I have not been able to fine a command to automatically calculate the running cumulative sum of a vector. Is there such a command? Example of current code: > eig$values [1] 678.365651 6.769697 2.853783 > prop<-eig$values/sum(eig$values) > prop [1] 0.986012163 0.009839832 0.004148005 >

Hello can anyone help, I have been running the following script to obtain a PCA plot but the end result is rather disappointing as the points are very very small and there are no titles etc geochemdata<-read.csv(file.choose(),header=TRUE) names(geochemdata) library(vegan) bstick<-function(n, tot.var=1) rev(cumsum(tot.var/n:1)/n) geopca<-rda(geochemdata, scale=TRUE) geopca

The following code in clusplot.default (package cluster) is in error: x1 <- cmdscale(x, k = 2, eig = TRUE) var.dec <- sum(x1$eig)/sum(diag(x1$x)) if (var.dec < 0) var.dec <- 0 if (var.dec > 1) var.dec <- 1 x1 <- x1$points x1 has components with names "points" and "eig", not "x", so

Hi all, I've been trying to figure out how to superimpose a histogram on a biplot that shows the relative contribution of each axis. I have been using the NIPALS function (http://biomserv.univ-lyon1.fr/~dray/ files/softwares/nipals.R) to run principal component analyses. Here is a toy example. source("http://biomserv.univ-lyon1.fr/~dray/files/softwares/nipals.R")

Hi folks, I've been trying to get shadow copy / previous versions working under sernet-samba-ad-4.1.6-7.el6.x86_64. I cannot see where I am going wrong. I have: [global] workgroup = COMPANY realm = COMPANY.LOCAL netbios name = DC server role = active directory domain controller server services = s3fs, rpc, nbt, wrepl, ldap, cldap, kdc, drepl, winbind,

Hi I need to extract the data returned by Eigen to plot the eigenvectors. However, when I try and eigv = eigen(covariance); it returns an object with the matrices containing eigenvalues and vectors.. how can I extract the eigenvector matrix from this?? When I try mat = eig["vectors"] it returns a matrix with the "$vectors" string on top , how can I remove this? code: > eig

Hello, I not have attached in this e-mail the zipped list of matrices I am using because it has 1 meg once zipped and anyway we cannot send attached files on r-help mailling list. First, after running the code that is written bellow, I realized that the printout of dudi.pca gives me for both of my matrices $cw length that are unequal to either respective $rank or $eig lengths. Is that

Dear all, I cannot seem to get the R functions step or stepAIC to perform forward or stepwise regression as I expect. I have enclosed the example data in a dataframe at the end of this mail. Note rubbish is and rnorm(17) variable which I have deliberately added to the data to test the stepwise procedure. I have used wateruse.lm<-lm(waterusage~.,data=wateruse) # Fit full model

Hi We have upgraded a very old Samba server to a new server running 3.0.22 on Ubuntu. The the workstations were Windows 98 on very old hardware and are now Windows XP on new hardware. Everything is now very fast compared to the old server, except for two applications - a payroll (Ace Payroll) and and finance system (Profax). Somewhat contrarily, these two apps are MUCH slower on the new server

Hi, Our Debian 3.01 file server, running Samba 3.0.0 is setup for guest access to all of the shares. Now I want to add one username, that is a private share for payroll. I added the username 'pam' with adduser, and configured smb.conf, see below. [Payroll] path = /home/pam comment = Private valid users = pam writeable = yes username level

Dear List, I am having trouble with some R code I have written to perform Redundancy Analysis (RDA) on a matrix of species abundance data (Y) and a matrix of environmental data (X). RDA is a constrained form of PCA and can be thought of as a PCA of the fitted values of a regression of each variable in Y on all variables in X. For info, the first use of RDA is in: Rao, C.R, 1964. The use and

Hi List, I'm following on the discussions regarding the use of exploratory factor analysis with binary data, from nabble: http://www.nabble.com/PCA--for-Binary-data-to11050371.html#a11092127 and i tried using FactorMineR package but unfortunately i have trouble understanding what are the types of the variables which AFDM is suitable for? I guessed from that example in ?AFDM that

Hello all, I'm using R1.7.1 on Linux, generating sammon-optimized MDS plots from distance matrices. This is a calculation I run routinely, often on sample sets of up to 100 samples. This time, with three samples, the sammon function returned an error (shown below), which I tracked down to the cmdscale function it uses to find a starting configuration. In short, cmdscale is returning NaN

Hi all, I had been running a Samba pdc using the standard tdbsam backend on CentOS 4.5 supporting 20+ WinXPsp2 workstations and a Win2k3 Terminal Server. This set up was running fine for arroung 18 months maybe more. The above system was built on ageing hardware that used to be a Win 2k server that died from OS corruption. Due to the expanding client needs there was good reason to replace

Hello to all, I have recently upgraded to SAMBA 3.4.2 on Solaris 10, and reconfigured it to use domain authentication (security = domain). We slapped guest authentication on most shares, with an explicit "valid users = ...." on a small number of sensitive shares. Due to the number of users we were looking at, we set up two UNIX groups "payroll" and "payoff" and

Dear all, I want to compute the eigenvalues of a complex matrix for some statistics. Comparing it to its matlab/octave sibling, I don't get the same eigenvalues in R computing it from the exact same matrix. In R, I used eigen() and arpack() that give different eigenvalues. In matlab/octave I used eig() and eigs() that give out the same eigenvalues but different to the R ones. For real

Hi Could anyone please suggest a good Payroll cum HR management open source rails application? Thanks Tom -- Posted via http://www.ruby-forum.com/. -- You received this message because you are subscribed to the Google Groups "Ruby on Rails: Talk" group. To post to this group, send email to rubyonrails-talk-/JYPxA39Uh5TLH3MbocFF+G/Ez6ZCGd0@public.gmane.org To unsubscribe from this