similar to: newbie question on dist

The following problem occurs in R 1.4.0 and 1.3.1 for Windows95, but not in R 1.2.0 for Windows95. The problem does not occur in R 1.4.0 for Linux PC, Linux Alpha and HP-UX. Sometimes, the type of 'Size' of an object of type 'dist' changes from integer into double. Running cmdscale on such a 'dist' object gives invalid results. I don't know what should be considered

Hi, i need a sample-code, for which the llvm alias-analysis finds a *must-aliases*. I have tried codes like followings. In all cases, i see just *may-aliases* when i use "opt -aa-eval -print-all-alias-modref-info foo.bc": Regards Raad 1 ========================================== void foo() { int i = 2; int& r = i; } 2 =========================================== void

I have a list of vectors $A "AB" "BC" "CD" $B "GF" "HG" "FH" "FJ" and I want to convert it into a dataframe of form A AB A BC A CD B GF B HG B FH B FJ Just can't quite come up with a nice "R" solution for it. Thanks, Sean

In applying multidimensional scaling, it seems to me that sometimes the underlying dimensionality of the matrix is 1. However I found a case where cmdscale failed when I tried k=1. Here it is: m<-matrix( c(.5,.81,.23,.47,.61, .19,.5,.06,.17,.28, .77,.94,.5,.74,.85, .53,.83,.26,.5,.64, .39,.72,.15,.36,.5), nrow=5) # BTW I think cmdscale uses only the lower triangle--how to enter only # that

Hi We have a few questions regarding the use of the "isoMDS" function. When we run "isoMDS" function using 60,000 x 60,000 data matrix, we get the following error message: ------------------------------------ cmdscale(d, k) : invalid value of 'n' Calls: isoMDS -> cmdscale ------------------------------------ We checked the source code of "cmdscale" and

> > On 9 Feb 2017, at 16:00, G?ran Brostr?m <goran.brostrom at umu.se> wrote: > > > > In my package 'glmmML' I'm using old C code and linpack in the optimizing procedure. Specifically, one part of the code looks like this: > > > > F77_CALL(dpoco)(*hessian, &bdim, &bdim, &rcond, work, info); > > if (*info == 0){ > >

I'm looking for someone to explain the difference between these procedures. The function prcomp() does principal components anaylsis, and the function cmdscale() does classical multi-dimensional scaling (also called principal coordinates analysis). My confusion stems from the fact that they give very similar results: my.d <- matrix(rnorm(50), ncol=5) rownames(my.d) <-

Sorry to bother you kind folks again with my questions. I am trying to learn as much as I can about all this, and I will admit that I don't have the proper background, but I hope that someone can at least point me in the correct direction. I have created a test matrix for what I want to do: s1 s2 s3 s4 s5 s1 10 5 0 8 7 s2 5 10 0 0 5 s3 0 0 10 0 0 s4 8 0 0 10 0 s5 7

> On 10 Feb 2017, at 14:53, G?ran Brostr?m <goran.brostrom at umu.se> wrote: > > Thanks to all who answered my third question. I learned something, but: > > On 2017-02-09 17:44, Martin Maechler wrote: >> >>>> On 9 Feb 2017, at 16:00, G?ran Brostr?m <goran.brostrom at umu.se> wrote: >>>> >>>> In my package 'glmmML'

Hello, I have a dilemma that I'm hoping the R gurus will be able to help resolve. For background: My data is in the form of a (dis)similarity matrix created from taking the inverse of normalized reaction times. That is, each cell of the matrix represents how long it took to distinguish two stimuli from one another-- a square matrix of 45X45 where the diagonal values are all zero (since this

Hello, We are trying to find a function to compute "stress" in our multidimensional scaling analysis of a dissimilarity matrix. We've used "dist()" to create the matrix and "cmdscale()" for the scaling. In order to determine the number of dimensions we would like to plot stress vs. dimensions. However, we cannot find a pre-made command. It seems that other

Hi. I know matrices that use distances between places works fine when using cmdscale. However, what about matricies such as: A B C D E A 0 1 23 12 9 B 1 0 10 12 3 C 23 10 0 23 4 D 12 12 23 0 21 E 9 3 4 21 0 i.e. matrices which do not represent physical distances between places (as they would not make sense for real distances such as the one above)

Full_Name: Laurent Gautier Version: 1.3.0-patched OS: IRIX 6.5 Submission from: (NULL) (130.225.67.199) Hello, The function La.eigen, called by cmdscale in the package mva behaves an unexplicable way (for me). The following lines show what happened. I tried the very same on linux, and it worked fine. >a <- matrix(c(1,2,3,2),3,3) >a [,1] [,2] [,3] [1,] 1 2 3 [2,]

Hello, I want to implement LLVM backend for a specific VLIW hardware. I am working on defining its instruction set, and assembly language. The hardware has two pipelines, int and float. Each pipeline can do 3 operations/cycle, 3 operations forms an instruction. One of the Integer Instruction looks like this: add Ri, Rj, Rk; add Rl, Rm, Rn; add Ro, Rp, Rq An int instruction and a float

Hi, I am looking for the best multidimensional configuration for my data (47*47 distance matrix). I ve tried classical metric (cmdscale) and non metric MDS (isoMDS, nmds) but it is now difficult to choose the best solution because of the uncertainties in the definitions of the "stress" function. So, same problem, several questions : 1. Statistical consideration : With

If tbl is an object of class 'dist', you can do this: a <- sammon(tbl, k=3) But you can't do this: b <- isoMDS(tbl, k=3) Wouldn't it be sensible to have identical interfaces to sammon() and isoMDS() ? I think all that would be needed is to change this: isoMDS <- function(d, y=cmdscale(d, 2), maxit=50, trace=TRUE) { ...into this: isoMDS <-

I have a dataframe like this (toy example): x y z "a" "a" 0 "a" "b" 1 "a" "c" 2 "b" "a" .9 "b" "b" 0 "b" "c" 1.3 "c" "a" 2.2 "c" "b" 1.1 "c" "c" 0 The observations are from a matrix like this: c 2.2 1.1 0.0 b 0.9 0.0

If do: > library("e1071") > example(svm) I get: svm> data(iris) svm> attach(iris) svm> ## classification mode svm> # default with factor response: svm> model <- svm(Species ~ ., data = iris) svm> # alternatively the traditional interface: svm> x <- subset(iris, select = -Species) svm> y <- Species svm> model <- svm(x, y) svm>

Dear All, I am in the somewhat unfortunate position of having to reproduce the results previously obtained from (non-metric?) MDS on a "kinship" matrix using Statistica. A kinship matrix measures affinity between groups, and has its maximum values on the diagonal. Apparently, starting with a nxn kinship matrix, all it was needed to do was to feed it to Statistica flagging that the

Hello all- I'm having some problems with memory consumption under R. I've tried increasing the appropriate memory values, but it keeps asking for more; I've even upped the heap size to 600M, significantly eating into swap (256M real, 500+M swap). So, performance slows to a crawl. What I'm trying to do is run isoMDS on a 4000x4000 matrix. My first question is, how much memory