similar to: variance components

How can I perform a Student-Newman-Keuls-Test for multiple comparision of means in R? (I did not manage to find any specific function in the libraries) Jörg Peter Baresel Technische Universität München Institut für Ackerbau und INformatik im Pflanzenbau D-85354 Freising Lange Point 51 [[alternative HTML version deleted]]

Being relatively new to R, I want to use R in connection with a MySQL database. After installing (apparently sucessfully) the package RMySQL (R Version 1.6.2) the result is as follows: > library(RMySQL) > con <- dbConnect(MySQL(), dbname= "test") Error: couldn''t find function "dbConnect" Thank You in advance for help Peter Jörg Peter Baresel Technische

Hi, under the new version R 2.7.1 I can't use WinEdt. It seems that there is a bug in that version. To use older versions of the package RWinEdt doens't bring success. The Use of R 2.7.0 doesn't cause that Problem. Regards, Hans-J?rgen Auinger -- ---------------------------------------------- Hans-J?rgen Auinger Lehrstuhl f?r Pflanzenz?chtung Wissenschaftszentrum Weihenstephan

Dear all, I'm analysing a split-plot experiment, where there are sometimes one or two values missing. I realized that if the data is slightly unbalanced, the effect of the subplot-treatment will also appear and be tested against the mainplot-error term. I replicated this with the Oats dataset from Yates (1935), contained in the nlme package, where Variety is on mainplot, and nitro on

Hi, I'm going to migrate from SAMBA 3.2 to 3.5 (Debian Lenny -> Squeeze). The server is an AD member and uses ACLs. Are there any preliminary steps to make the upgrade as smooth as possible? What kind of problems could I expect? Best Regards, Christian -- __________________________________________ Christian Reischl Fraunhofer Institut f?r Verfahrenstechnik und Verpackung Giggenhauser

Dear List, I´m trying to implement the functionalities from WEKA into my modeling project in R through the RWeka package. In this context I have a slightly special question about the filters implemented in WEKA. I want to convert nominal attributes with k values into k binary attributes through the NominalToBinary filter ("weka.filters.supervised.attribute.NominalToBinary"). But

Hello, I want to draw 3D cylinder objects in R. Given is the length and the diameter of the cylinder. Has anybody an example? Thank you very much! Best regards -- ---------------------- Dr. Hans-Joachim Klemmt Forstoberrat Organisationsprogrammierer IHK Bayerische Landesanstalt f?r Wald und Forstwirtschaft zugewiesen an Lehrstuhl f?r Waldwachstumskunde Technische Universit?t M?nchen

Hi list, I´m working on a predictive modeling task using the caret package. I found the best model parameters using the train() and trainControl() command. Now I want to evaluate my model and make predictions on a test dataset. I tried to follow the instructions in the manual and the vignettes but unfortunately I´m getting an error message I can`t figure out. Here is my code: rfControl <-

Hi, I would like to use the beanplot() function from the beanplot package. Unfortunately, I can't find out how to suppress the dashed horizontal line, that shows the overall mean. In the help I've found the argument "overallline", but it only allows for "mean" or "median" . I have tried overallline = F, overallline="n", and

Hi! In a current project, I am fitting loess models to subsets of data in order to use the loess predicitons for normalization (similar to what is done in many microarray analyses). While working on this I ran into a problem when I tried to predict from the loess models and the data contained NAs or NaNs. I tracked down the problem to the fact that predict.loess will not return a value at all

Hi, I use a mboost model to predict my dependent variable on new data. I get the following warning message: In bs(mf[[i]], knots = args$knots[[i]]$knots, degree = args$degree, : some 'x' values beyond boundary knots may cause ill-conditioned bases The new predicted values are partly negative although the variable in the training data ranges from 3 to 8 on a numeric scale. In order to

Dear ExpeRts, I am trying to plot a bunch of growth curves and would like to get some more control over groups and line colors than I seem to have. Example: # make some data dat <- Orange dat$group <- ifelse(dat$Tree%in%c('1','4','5'), 'A', 'B') # plot xyplot(circumference~age, dat, groups=group) # now use lines to make the growth curve more

Hello List, I`m trying to make prediction using a bagged tree with the package ipred. I tried to follow the manual but I`m getting an error message. Also browsing through the list-archive I didn`t find any hint. Maybe someone can help me? selbag <- bagging(SOIL_UNIT ~., data=traindat.bin, coob=TRUE) Error in dimnames(X) <- list(dn[[1L]], unlist(collabs, use.names = FALSE)) :

Dear ExpeRts,t I am trying to read tab delimted data produced by somewhat brain dead software that seems to think it's a good idea to have an extra tab character after the last column - except for the header line. As explained in the help page, read.delim now assumes that the first column contains the row.names (which is not even wrong) but now and all col.names get shiftet by one column.

Hello Everybody, we have authentication problems with the mentioned configuration. Current situation: We have two Windows 2008 R2 domains (currently on 2003 level) in separate forests. Recently we created a two-way trust with selective authentication between them. A Debian Squeeze LTS machine running SAMBA 3.6.6 (latest Backports version) is member of our primary domain and provides file

Dear List, I´m working on a classification problem. My response has 60 levels. I`m very interested in boosted trees like AdaBoost or gradient boosting machine as implemented in the package "gbm". Unfortunately gbm is only applicable for 2-class problems. Is anybody out there who can help me? Is there a way to use gbm() for multi-class problems? Maybe there is a way to transform my

-------- Original-Nachricht -------- Datum: Wed, 21 May 2008 18:12:59 +0800 Von: "Jon Ervine" <jon.ervine@gmail.com> An: "Jörg Kühne" <joerg.kuehne@gmx.net> Betreff: Re: [Xen-users] opensuse10.3 host and windows XP guest via Xen/Vbx 2008/5/21 "Jörg Kühne" <joerg.kuehne@gmx.net>: > Hello > > did anyone may made some experience with the

Hallo everyone I have two univariate time series (class ts) describing the same variable. They have the same resolution, but span different periods in time with a big gap in between. I need to append one to the other such that they are one object, with the gap filled with NAs. The method ts.union produces a multivariate time series where the time axis is correct, but the individual time series

Full_Name: Jörg Polzehl Version: 1.8.0 OS: Windows XP Submission from: (NULL) (62.141.176.1) I encountered a problem when playing with the mle library and specifying negative starting values for the parameters. The reason seems to be an incorrect behaviour of function formals: glike<-function(a=1,b=1,c=1) a > formals(glike) $a [1] 1 $b [1] 1 $c [1] 1 > unlist(formals(glike)) a b c

Dear R users, I have been doing rotation on Principal components analyse, with varimax function and promax. Following those changes, I cannot find now how to get the spatial evolution and the variance. Indeed, with the function princomp that has been used to get the principal components, arguments such as $scores and $sdev were available to get the spatial evolution and the variance, with