search for: dsq

I'm looking for something in R to fit a multivariate linear model robustly, using an M-estimator or any of the myriad of other robust methods for linear models implemented in robustbase or methods based on MCD or MVE covariance estimation (package rrcov). E.g., one can fit an mlm for the iris data as: iris.mod <- lm(cbind(Sepal.Length, Sepal.Width, Petal.Length, Petal.Width) ~ Species,

...he Spearman rho between two ranked lists. I am getting results with cor.test that differ in comparison to my own spearman function: > my.spearman function(l1, l2) { if(length(l1) != length(l2)) stop("lists must have same length") r1 <- rank(l1) r2 <- rank(l2) dsq <- sapply(r1-r2,function(x) x^2) 1 - ((6 * sum(dsq)) / (length(l1) * (length(l1)^2 - 1))) } Perhaps I'm doing something wrong in that code, but it's a pretty straightforward calculation, so it's hard to see what, especially with rank() handling the ties correctly. One example...

I'm using rseptbeta for w95, I tried to do: library(mva) data(quakes) quakes<-quakes[sample(seq(1:1000),100),1:5] cmd.dsq<-cmdscale(dsq) Error in .C("dblcen", x, as.integer(n)) : C/Fortran function not in load table how can I solve this problem? Andrea Rossetti _______________________________________________________ Statistica & Informatica per la Gestione delle Imprese | Universit? degli Studi di P...

I'm using the WINDOWS 95 version of R (with midified menu "exe.zip"): I tried to do library(mva) data(quakes) quakes<-quakes[sample(seq(1:1000),100),1:5] cmd.dsq<-cmdscale(dsq) Error in .C("dblcen", x, as.integer(n)) : C/Fortran function not in load table how can I solve this problem? -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send...

...t data. Two species of trees were planted in different mixtures in 10 (nearly-adjacent) plots and measured on four occasions over 10 years. The model is constructed in terms of the diameter increments (per year; DI) in the 3 intervals, in which DI is modelled as a function of mid-interval D and DSQ. The details of the fixed part of the model are not so important here, but four pertinent variables are distance-dependent competition indices: C_spp1fromspp1, C_12, C_21, C_22. The random structure is: random = ~ 1 | PLOT/TREE. More complex random structures aren't required, there's n...

Hi, Thanks for the help. I tried this out, but on x86_64 and with llvm 2.6 and llvm-gcc 4.2.1, I dont seem to get an aligned variable. This is the function definition that I have define internal fastcc void @walksub(%struct.hgstruct* noalias nocapture sret %agg.result, %struct.node* %p, double %dsq, %struct.hgstruct* byval align 64 %hg, i32 %level) nounwind { And these are the call sites, call fastcc void @walksub(%struct.hgstruct* noalias sret %2, %struct.node* %tmp399, double %tmp398, %struct.hgstruct* byval align 64 %1, i32 0) nounwind call fastcc void @walksub(%struct.hgstruct* notcc...

Hi Arushi, > If I want a particular alignment on a ByVal argument, should I specify > it as a function attribute, or should I also add the attribute to > every callsite that calls that function. in theory you should do both. For example, if the call is an indirect call (i.e. to a function pointer where the function is not known) then it should be clear that you need to put the

Hi, If I want a particular alignment on a ByVal argument, should I specify it as a function attribute, or should I also add the attribute to every callsite that calls that function. Thanks, Arushi