similar to: 49 histograms on one page

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 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)

Hello list, I want to apply svm from library e1071, and I want to supply class weights. I do not really understand the help entry (and there is no example) class.weights: a named vector of weights for the different classes, used for asymetric class sizes. Not all factor levels have to be supplied (default weight: 1). All components have to be named. I have two

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 group, I would like to document more than one function on one help page using the Rd language. I tried something like \value{ \code{foo1} returns a list with components \item{arg1}{Argument 1} \item{arg2}{Argument 2} \code{foo2} returns a list with components \item{arg3}{Argument 3} \item{arg4}{Argument 4} } Unfortunately this makes the text "\code{foo2} returns a list

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, 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,]),

Dear group, I want to use the mclust package on large data, and therefore I want to use a subset in the initial clustering phase. From help(mclust): k: If `k' is specified, the hierarchical clustering phase will use a sample of size `k' of the data in the initial hierarchical clustering phase. The default is to use the entire data set. m2 is a

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: >

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

Hello list, I try to optimize some R code and it turns out that the function as.matrix takes a lot of time in my code. There is only one reason why I need as.matrix: My code should work for matrices with ncol=p and p should be allowed to be 1 or larger. Now I have a matrix x with dim(x)=(n,p), and I need to work with the submatrix y <-x[gv,], (gv a vector of n logicals) and to calculate

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

Dear list, I use R 1.3.0 on Solaris. Today I realized that the html help search engine does not work at the moment (this means that I always get the search engine page again if I try to specify a search word, and click "search" or one of the prespecified topics), presumably since the latest netscape update. Do you have an idea what can be wrong? Perhaps this must be fixed not by the

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/

Dear list, on p. 11 of "Writing R-extensions" I read "See the ``Guidelines for Rd-files'' for guidelines ... which should be useful..." Where can I find these guidelines? Best, Christian -- *********************************************************************** Christian Hennig Seminar fuer Statistik, ETH-Zentrum (LEO), CH-8092 Zuerich (current) and Fachbereich

Dear list, what is the fastest way to compute a multivariate mean and cov-matrix? I presume that the mean is computed in cov, so it may be a waste of time to compute the mean first and then a second time inside of cov. Is it faster to use cov.wt, which gives cov-matrix and center? And: If mean and cov should be computed on a part of the data, is it faster to use cov.wt with some weights zero, or

Hi, I want to fit a kernel density estimator by bkde of library KernSmooth. I need only the density value at the point 0. I do not understand the following behaviour: > q <- rnorm(100) > bkq <- bkde(q, bandwidth=0.11, gridsize=1, range.x=c(0,0)) Error in 0:L : NA/NaN argument > bkq <- bkde(q, bandwidth=0.11, gridsize=1, range.x=c(-1,1)) > bkq $x [1] -1 $y [1] NA > bkq

Hi, exist a special cluster-analysis algorithms which can work with NA's. a further "problem" is that i want cluster variables not cases to identify special variable-set's. Is it a common way turn the data.frame and use kmeans,because this works with NA's, or have anybody another method for finding "variable-sets" , with exception of factor analysis. thanks for

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, 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