similar to: New package: expert 1.0-0

Don't use aggregate's simplify=TRUE when FUN() produces return values of various dimensions. In your case, the shape of table(subset)'s return value depends on the number of levels in the factor 'subset'. If you make B a factor before splitting it by C, each split will have the same number of levels (2). If you split it and then let table convert each split to a factor, one

To whom it may concern, I am currently writing a program where I need to use function rep. The results I get are quite confusing. Given two vectors A and B, I want to replicate a[1] b[1] times, a[2] b[2] times and so on. All the entries of vector B are positive integers. My problem comes from the fact that if I sum up all the elements of B, I get a certain value x(for example 10000). And if

Hi Alejo, According to my knowledge the two plots are different because in the first one a point belongs to a group depending on its group in the data whereas in the second plot a point belongs to the group predicted by the linear discriminant analysis. I hope somebody will correct me if I am wrong. Alain Alejo C.S. wrote: > Hi Alain, this is the code: > > > library(MASS) >

Dear R users, When I use aggregate with table as FUN, I get what I would call a strange behaviour if it involves numerical vectors and one "level" of it is not present for every "levels" of the "by" variable: --------------------------- > df <- data.frame(A=c(1,1,1,1,0,0,0,0),B=c(1,0,1,0,0,0,1,0),C=c(1,0,1,0,0,1,1,1)) >

Hello, I am using the recode() function in Rcmdr and the result is not what I expect so I am almost sure I did something wrong but what... > test <- data.frame(x=1:10) > library(car) > recode(test$x,'1:5=0 ; else=1', as.factor.result=TRUE) [1] 0 0 0 0 0 1 1 1 1 1 Levels: 0 1 BUT > library(Rcmdr) # recode from the car package is now masked Now I recode test$x

Hi, I drew an ellipse with the package ellipse. Now I would like to know if a point is inside the ellipse. Is any R functions to do it without computing the equation of the ellipse manually? Thanks. For example, if I do "plot(ellipse(0.8), type = 'l')", I would like to know if (0,1) belongs to the drawn ellipse. Regards, Alain -- Alain Guillet Statistician and Computer

Hi, I want to make my R program run in batch under Windows XP. To do so, I create a bat file with the command RCMD BATCH --vanilla program.R program.out and I use the bat file with the scheduled task of Windows XP. Then I log off. It works up to the log off of another user on the same computer with R-2.9.1 but this problem doesn't appear with R-1.9.1 on the same machine. Is anything

Dear useRs, I have written a very simple function to compute some probabilities on words (function is below). The function includes a which() statement applied to a vector of characters (word.split): sapply (word.split, function(x) which(letters==x)). This statement worked as expected when used outside the global function : > word <- "hello" > (word.split <-

hi! after i always changed to windows for just one stupid programm (from my university - important) which didn't work on wine, i got sick of it. i tried it once more with the newest version hoping it might be work and it already got farther then the last time (which was like december?!?!) but it doesn't work. The program seems to be ToolBook based (http://www.toolbook.com) and has no

Hello, I'm currently using the segmented package of M.R. Muggeo to fit a two-slope segmented regression. I would like to constrain a null-left-slope, but I cannot make it. I followed the explanations of the package (http://dssm.unipa.it/vmuggeo/segmentedRnews.pdf) to write the following code : fit.glm <- glm(y~x) fit.seg <- segmented(fit.glm, seg.Z=~x,psi=0.3) fit.glm

Dear Rusers, I tried to implement a function comparing mean scores between one subject (the patient) and a group a control subjects. The function returns attended results, but I also obtained the following warning : Warning message: the condition has length > 1 and only the first element will be used in: if (substr(fp, 1, 1) == "<") fp else paste("=", fp) Maybe

Dear R-devel, googling for the single letter R yields R-home page as the firt hit, which is extremly nice. By clicking on the figure you get the code of the "Winner of the R Homepage graphics competition 2004." By copy/pasting in your R console it doesn't work because it is impossible to install the mva package. This is my point. As documented in ONEWS the reason is simple:

Hi, I try to setup a simple file sharing with some ACL for a group of users. My Samba 3.0.14 server is running on Fedora Core 4 using pam_ldap. User authentication works correctly with ldap and my nsswitch is correctly configured to use the ldap as "extended list" of users and groups. For exemple, when I type getent passwd or group, I have the complete list of system users and ldap

Dear R specialists, I have a multivariate statistics question that I want to submit to the R community (which conveys a very good statistical knowledge). I need to perform an optimal scaling based on a discrete variable and an ordinal variable. The discrete variable, Area, defines a geographical area. The ordinal variable, EducationLevel, describes the education level of individuals (the

Hello, I have a question about addTaskCallback It is said that the argument visible allows to know wether the result of the top-level evaluation was printed or not. Nevertheless, in the following example, I encountered a problem, as it seems not all printed objects are visible... Could someone tell me where is the tip? TIA

Hi All, I would like to use Kennard-Stone algorithm for splitting a dataset. > mydata <- read.csv(url("http://www.ats.ucla.edu/stat/r/dae/binary.csv ")) > library("soil.spec") > ken.sto(mydata,per.n=0.3) Error in ken.sto(mydata, per.n = 0.3) : subscript out of bounds I found that other people run into this problem as well:

Hi, I created a principal component plot using the first two principal components. I used the function princomp() to calculate the scores. now, I would like to superimpose an ellipse representing the center and the 95% confidence interval of a series of points in my plot (as to illustrate the grouping of my samples). I looked at the ellipse() function in the ellipse package but can't get it

Hi to all R/S gurus, I am developping a library to export R/S objects to HTML (there was developments of this library past year, by myself and Mathieu Ros, but we never finished it). For that, I wrote HTML.* functions handled by the generic method HTML. At this step, allmost all objects exportations are possible: matrix, atomic, data.frame, functions and even plots... but I still don't

Dear fellow users Are there any Vietnamese language resources for beginners of R? If so, I would be interested in hearing from people who have had experience with them and which are better (if there is more than one). I am involved with an aid project in Vietnam, and would like to move the scientists involved from using Excel for 'analysis' to R. Thanks .... Peter Alspach The

Hello, After an upgrade to RedHat 9 and the supplied samba-[|commoan-|client-]2.2.7a-8.9.0 rpms we've discovered a problem with smbmount. It sometimes invokes itself twice and then hangs. Sample script to test: #!/usr/bin/perl # testing Samba problem print "$$\n"; `/usr/bin/smbmount //cmc-apps/settlements ./t-samba -o