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

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

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

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

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

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

Dear All: I do see I have weeks of codes in my console when I check with my arrow up keys. I have been clearing them with Control L but it seems to clear it clear the screen temporally. I do see the previous codes again when I open R the next day, after quitting the session! Q: How do I clear this? Thanks; YA -- View this message in context:

This may be asking too much, but I'm wondering if anyone has a solution (even a hack) for creating multiple (overlay) plots in an Sweave file and post-processing the overlays in beamer appropriately. For example, suppose I have a series of figure blocks in my .Rnw file: <<plot1,fig=TRUE>>= [stuff] @ <<plot2,fig=TRUE>>= [stuff] @ <<plot3,fig=TRUE>>=

I have a 64-bit windows box - Intel Xeon CPU E7340 @ 2.4GHz 31.9GB of RAM I have R 2.11.1 (64bit) running on it. My csv data is 3.6 GB (with about 15 million obs, 120 variables.) ------------------------------------------------ I have successfully imported the data above into R. No problem. Now I am trying to run 'rpart' on my data. But I got the following error : Error: cannot

Dear R experts, Suppose I have a data frame of three variables: > foo <- data.frame(row=1:5, col=1:3, val=rnorm(15)) > foo row col val 1 1 1 -1.00631642 2 2 2 0.77715344 3 3 3 0.17358793 4 4 1 -1.67226988 5 5 2 1.08218836 6 1 3 1.32961329 7 2 1 -0.51186267 8 3 2 -1.20990127 9 4 3 -0.57786899 10 5 1 0.67102887 11 1 2