search for: fm3a004

Sorry if this bug is fixed in the meantime, but what is this?? > m [,1] [,2] [,3] [1,] 1 2 1 [2,] 0 2 0 [3,] 1 2 1 > eigen(m) $values [1] 2 2 0 $vectors [,1] [,2] [,3] [1,] -4.194304e+06 0.7071068 -0.7071068 [2,] -4.656613e-10 0.0000000 0.0000000 [3,] -4.194304e+06 0.7071068 0.7071068 > eigen(m, symmetric=T) $values [1]

Dear list, usually, I can break a running R process with Ctrl-C. This does not seem to work if R was started from within ESS. Is there anything analogous? Thanks, Christian *********************************************************************** Christian Hennig University of Hamburg, Fachbereich Mathematik - SPST (Schwerpunkt Mathematische Statistik und Stochastische Prozesse, Zentrum fuer

Dear list, I am faced with the following model: y=E+P+M+H+PxE+Error y is a response, E and P are factors with fixed effects. M is a random effect nested in P and H is a random effect nested in M. PxE is interaction of P and E. It seems that I should fit such a model with the lme function of library nlme, but I was not able to figure out from the help page how to specify the formula. In the

Dear list, I tried to install R-1.1.0 under SuSe Linux 5.2. ./configure --prefix=/usr/applic/packages/R make make install worked without error. make check led to make[1]: Entering directory `/export/old/usr/applic/packages/R/R-1.1.0/tests' make[2]: Entering directory `/export/old/usr/applic/packages/R/R-1.1.0/tests' make[3]: Entering directory

Hi, this is rather a (presumed) bug report than a question because I can solve my personal statistical problem by working with hclust instead of agnes. I have done a complete linkage clustering on a dist object dm with 30 objects with agnes (R 1.8.0 on RedHat) and I want to obtain the partition that results from a cut at height=0.4. I run > cl1a <- agnes(dm, method="complete")

Hi, this is rather a (presumed) bug report than a question because I can solve my personal statistical problem by working with hclust instead of agnes. I have done a complete linkage clustering on a dist object dm with 30 objects with agnes (R 1.8.0 on RedHat) and I want to obtain the partition that results from a cut at height=0.4. I run > cl1a <- agnes(dm, method="complete")

Hi people, Does anybody know some Density-Based Method for clustering implemented in R? Thanks, Fernando Prass _______________________________________________________

Hi, I want R to be installed on a UNIX network (Solaris). I am not the system administrator and so I cannot do it myself. The system administrator wants to know which packages I want, and it may be a lot. Is there an easy way to download and install all packages at once? Is it a good idea? (There is a good chance that I do not need the some few packages that do not install well during such a

Hi, this is to announce a new version (1.1-2) of my package fpc. Apart from the stuff already present in the older version (methods for fixed point clustering and clusterwise regression, somewhat bug-cleaned and with faster examples) there is now a function cluster.stats, which computes some distance-based statistics often used for cluster validation, description and decision about the number of

Hi, I have run rpart to construct a regression tree. Is there any simple method to draw a nice picture of it, as it is usually done in books and paper to visualize the tree? Thank you, Christian *********************************************************************** Christian Hennig Fachbereich Mathematik-SPST/ZMS, Universitaet Hamburg hennig at math.uni-hamburg.de,

Hi, this is to announce a new version (1.1-2) of my package fpc. Apart from the stuff already present in the older version (methods for fixed point clustering and clusterwise regression, somewhat bug-cleaned and with faster examples) there is now a function cluster.stats, which computes some distance-based statistics often used for cluster validation, description and decision about the number of

Hi, is there any method for goodness of fit testing of an (as general as possible) univariate distribution with parameters estimated, for normal, exponential, gamma distributions, say (e.g. the corrected p-values for the Kolmogorov-Smirnov or Chi-squared with corresponding ML estimation method)? It seems that neither ks.test nor chisq.test handle estimated parameters. I am aware of function

Dear list, dear Frank, I try to fit a Weibull survival regression model with package Design: sclear <- psm(sobj~V1+V2,dist="weibull") sobj is a one-dimensional survival object (no event indicators), V1 and V2 are factors. I get the following result: Parametric Survival Model: Weibull Distribution psm(formula = sobj ~ V1 + V2, dist = "weibull") Obs Events

Hi, is there a test for the H0 skewness=0 (or with skewness as test statistic and normality as H0) implemented in R? Thank you, Christian *********************************************************************** Christian Hennig Fachbereich Mathematik-SPST/ZMS, Universitaet Hamburg hennig at math.uni-hamburg.de, http://www.math.uni-hamburg.de/home/hennig/

Hi, here are three results of random forest (version 4.0-1). The results seem to be more or less the same which is strange because I changed the classwt. I hoped that for example classwt=c(0.45,0.1,0.45) would result in fewer cases classified as class 2. Did I understand something wrong? Christian x1rf <- randomForest(x=as.data.frame(mfilters[cvtrain,]),

Hi, I would like to use the format function to get numbers all with three digits to the right of the decimal point, even in cases where there is no significant digit left. For example, I would like to get c(0.3456789,0.0000053) as "0.346" "0.000". It seems that it is not possible to force format to print a "0.000", i.e. without any significant decimal places. Is it

Hi, the following strange error appears when I use qda: > qda1 <- qda(as.data.frame(mfilters[cvtrain,]),as.factor(traingroups)) Error: function is not a closure That's also strange: > qda1 <- qda(mfilters[cvtrain,],as.factor(traingroups)) Error in qda.default(mfilters[cvtrain, ], as.factor(traingroups)) : length of dimnames must match that of dims Some backgroud: >

Dear list, the following should fit the model log(PEPC.Wert)=fKTemp+fHerk+interaction(fKTemp,fHerk)+fMub+error, where fKTemp, fHerk are fixed effect factors and fMub is a random effects factor nested in fHerk (values are different fopr different values of fHerk). > logpepcr1 <- lme(log(PEPC.Wert) ~ fKTemp*fHerk, random= ~ 1 | fMub, na.action=na.omit) The following should be without the

Dear list, (also sent to Roger Bivand, but perhaps somebody of you can help me also) I am trying to use package spdep for fitting an SAR model with errorsarlm. However, I am not sure how to make a valid nb object out of my neighborhood. As far as I have seen, there is no documentation for nb.object. I have done the following: class(pschmid$nb) <- "nb" # pschmid is a prab object as

Hi list, not really a problem to solve but a strange unwanted behavior (at least I do not like it): I'm loading library(multiv) to use sammon. Later I load library(MASS) which tells me that sammon is masked now, because there is also a MASS version of sammon. Since I believe that the MASS version is at least as good as the multiv version, I would also use the MASS version, but help(sammon)