similar to: ACL and user list on windows XP

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

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

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

expert: Modeling Without Data Using Expert Opinion Expert opinion is a technique to do statistical modeling when data is scarse (e.g. accidents in nuclear plants) or even absent, at least for the analyst (e.g. confidential settlements in liability insurance). Opinions on the distribution of the decision random variable is sought from experts in the field. The experts give their opinion in the

expert: Modeling Without Data Using Expert Opinion Expert opinion is a technique to do statistical modeling when data is scarse (e.g. accidents in nuclear plants) or even absent, at least for the analyst (e.g. confidential settlements in liability insurance). Opinions on the distribution of the decision random variable is sought from experts in the field. The experts give their opinion in the

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

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

Dear R users, I have been doing some eof analysis using princomp function, then the eof results were rotated with varimax and promax functions. Those functions are working fine. However, after their uses, I cannot find how to obtain the temporal evolution and the variance of the obtain components, as varimax and promax does not have the arguments $scores and $sdev. If any body, know how to

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

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

hi. i''ve got a couple migrations where i''ve added foreign keys and tried to update them via the migration but my items don''t save. i''ve used save! in the hopes that an exception gets raised but no luck. the migrations run without errors but the data isn''t in the db. can anyone point out what i might be doing wrong? thanks! code below. class

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:

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

Dear R-help, I was comparing SAS (I do not know what version it is) and R (version 2.6.0 (2007-10-03) on Linux) survival analyses with time-dependent covariates. The results differed significantly so I tried to understand on a short example where I went wrong. The following example shows that even when argument 'method' in R function coxph and argument 'ties' in SAS procedure

Hi, Can someone tell me how to create .rda data files in R so that they can be used contributed package? Thanks Deli Wang [[alternative HTML version deleted]]

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

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