similar to: html command index

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

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]

Greetings RListers, I have a data set containing two types of outcomes; success and failure. Associated with each outcome are 12 different measurements. I'm trying to find out, for example, if some of the 12 measures are associated more with success or failure or, if there's any relationship at all between the measures and the outcomes (success or failure). I don't have (as yet)

I'm attaching an small text file on "Getting your stuff organized in R". (Sorry if sending an attachment is not considered a correct etiquette in r-help, but this is only 7911 bytes, plain ascii text and I cannot post it in a web page at the moment). Probably all the information in this document is scattered in one or more R introduction guides, but I think that it is useful to have

Dear r-users, I am looking for an implementation of the Naive Bayes classifier for a multi-class classification problem. I can not even find the Naive Bayes classifier for two classes, though I can not believe it is not available. Can anyone help me? Uschi -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

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

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

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

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, 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 want to do something like this: par(mfrow=c(7,7)) for (i in 1:49) hist(RATDACOM[SUBJNO==i],breaks=0.5+(0:6), main="",xlab="",ylab="",xaxt="n",yaxt="n") (Don't think about what RATDACOM and SUBJNO are.) I get an error Error in plot.new() : Figure margins too large. 36 histograms with mfrow=c(6,6) work. But the 36 histograms

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

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

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