similar to: dataframes and factors

Hi! The postscript device in 0.90 does not produce symbols (pch=) as it should. E.g.: postscript(file="test.ps",width=8,height=6.4,horizontal=FALSE) plot(1:10,1:10,pch=1:5) legend(1,9,letters[1:5],pch=1:5) dev.off() The result looks quite different compared to the X-windows version. The symbol 1 (pch=1) is much to small, while the other symbols seem to big. Is it a bug? How can I

Dear list, I want to perform a biplot, using customized titels for the x and y axis. Setting xlab="" and ylab="" resulted in an error, e.g.: > data(USArrests) > biplot(princomp(USArrests),xlab="",ylab="") Error in biplot.default(t(t(scores[, choices])/lam), t(t(x$loadings[, : length of dimnames[1] not equal to array extent > How do I

Hi everyone I'm trying to generate a vector of the form a <- c(1,2,3,4,1,2,3,1,2,1) where n = 5 in a general way, where n can be any positive integer. I've run out of ideas. Does anyone have any suggestions? Thanks Karen Karen Kotschy Centre for Water in the Environment University of the Witwatersrand Johannesburg Tel: 011 716-2218

Hi everyone I'm trying to generate a vector of the form a <- c(1,2,3,4,1,2,3,1,2,1) where n = 5 in a general way, where n can be any positive integer. I've run out of ideas. Does anyone have any suggestions? Thanks Karen Karen Kotschy Centre for Water in the Environment University of the Witwatersrand Johannesburg Tel: 011 716-2218

Hi! Since installation of R 0.90 there are problems with the size of text and numbers also in the x11-device (I reported about similar problems with the postscript device recently). For numbers and letters a different size is used: plot(1:10,1:10,type="n") text(2,8,"F5F5F5F5F5") text(2,7,"F5F5F5F5F5",cex=2) The same problem is also found in the postscript

Hi! Since installation of R 0.90 there are problems with the size of text and numbers also in the x11-device (I reported about similar problems with the postscript device recently). For numbers and letters a different size is used: plot(1:10,1:10,type="n") text(2,8,"F5F5F5F5F5") text(2,7,"F5F5F5F5F5",cex=2) The same problem is also found in the postscript

Dear all, first, I would like to thank R-core team for making so a great software available. Thank you all, really! Using R (1.0.0) on a Win98 platform, I need an x-server to run xgobi (through R xgobi interface). I've been using a MI/X server but trying to upgrade I noticed that MI/X is no longer available for free. Similarly, eXceed is available at evaluation version (for a limited time

Dear R-developpers, currently the function window shows the following behaviour, if start or end is out of range: > x <- ts(1:14,start=c(1950,2),frequency=4) > y <- window(x,start=c(1950,1),end=c(1953,4)) Warning messages: 1: start value not changed in: window.default(x, ...) 2: end value not changed in: window.default(x, ...) Instead of this I would like it to fill with NAs the

From: "Wolfgang Koller" <koller2@fgr.wu-wien.ac.at> Organization: University of Economics, Vienna Date: Tue, 23 Nov 1999 16:29:22 CET (+0100) MIME-Version: 1.0 Content-type: text/plain; charset=US-ASCII Content-transfer-encoding: 7BIT CC: r-bugs@biostat.ku.dk, r-help@stat.math.ethz.ch References: <153DF9230E8@fgr.wu-wien.ac.at> X-mailer: Pegasus Mail

Dear R-Users! My aim is to produce boxplots without the outliers included in the plot. I started to write a function that looks something like: myboxplot <-function(x,fa) { bpdata <- boxplot(x~fa,plot=FALSE) bpnames <- names(bpdata) for (JJ in bpnames) { command <- paste("bpdata$",JJ,"$out <- numeric(0)",sep=""); eval(command) }

Hello! When I know the vector of the variance of the disturbances (i.e. the structure of heteroskedasticity), say Var(u_{i})=v_{i}, what is the weights I should use as argument to lm(): M <- lm(y~x,weigths=1/v) or M <- lm(y~x,weights=1/(v^0.5)) ??? In the help pages I did not find a clear answer to this question, so please could someone help me! Thanks, Wolfgang Koller

================================================= This is an automated summary of the status of the R-bugs repository. Note that this may be neither complete nor perfectly correct at any given instance: Not all bugs are reported, and some reported bugs may have been fixed, but the repository not yet updated. Some bug fixes are difficult to verify because they pertain to specific hardware or

Dear R-Users and Developpers, Currently R does not include functions to compute the skewness and kurtosis. I programmed it myself in the following way, but probably *real* programmers/statisticians can do that better: mykurtosis <- function(x) { m4 <- mean((x-mean(x))^4) kurt <- m4/(sd(x)^4)-3 kurt } myskewness <- function(x) { m3 <- mean((x-mean(x))^3) skew <-

Full_Name: Wolfgang Koller Version: 1.0.0 OS: linux Submission from: (NULL) (137.208.7.48) > data(UKLungDeaths) > ts.union(mdeaths, fdeaths,dframe=TRUE) Error in names<-.default(*tmp*, value = nmsers) : names attribute must be the same length as the vector > cbind(mdeaths, fdeaths,dframe=TRUE) Error in names<-.default(*tmp*, value = nmsers) : names attribute

Thanks to all who replied. Mostly the suggestion was to add the argument horizontal= FALSE to the postscript function. However, I was trying to generate a landscape graph. The "fault" lies with gs as Peter Dalgaard pointed out. Under Windows 95, gsview has an option to view landscape or portrait and therefore this problem does not arise on that platform.

>>>>> "WoKo" == Wolfgang Koller <koller2@fgr.wu-wien.ac.at> writes: WoKo> Hi! The postscript device in 0.90 does not produce symbols WoKo> (pch=) as it should. E.g.: WoKo> postscript(file="test.ps",width=8,height=6.4,horizontal=FALSE) WoKo> plot(1:10,1:10,pch=1:5) WoKo> legend(1,9,letters[1:5],pch=1:5) WoKo>

Dear Thomas, I want to use evaluate effect of Age using restricted cubic form in the FGR function as Fgr.crr <- FGR(Hist(time, event) ~ rcs(Age_years), data=dat) It provides error. " Error in parse(text = termtext, keep.source = FALSE): .... 1: response ~ rcs(Age_years Do I need to change any of the R code? Regards Amalraj Raja The University of Aberdeen is a charity

Dear R-helpers, I am looking for a segmentation package that gives some "tree identifier" as output for every observation in the data set (my response variable is binary). I have skimmed through "rpart", "ada" and "adabag": The output "trees" gives you the formula, but I have to run several thousand segmentations on different data sets and it

Hello list, I am writing a package which builds on a function ("foo") of another package. I only need that particular function. What is the state of the art to do this? Do I need a NAMESPACE file in my package? Currently I still don't have one. Or can I do this with the Imports field in my DESCRIPTION file? Reading chapter 1 of the "Writing R Extensions" manual I do

Dear Geskus, I want to develop a prediction model. I followed your paper and analysed thro' weighted coxph approach. I can develop nomogram based on the final model also. But I do not know how to do internal validation of the model and subsequently obtain calibration plot. Is it possible to use Wolbers et al Epid 2009 approach 9 (R code for internal validation and calibration) . It is