similar to: generate random number

How can I access a function created in one data directory in another? I have been sourcing the text file that contains the function each time I go to another data dir. This is cumbersome, and I suspect there is a way to make functions global. Thanks in advance, S. David White Dept of Speech and Hearing Science Ohio State University Columbus, Ohio

Hi all, I have a problem which I am biting my teeth into unsuccessfully: (x <- data.frame(S1=c("a","b","d","F"),N1=c(2,4,6,9),N2=c(6,NA,0,6))) S1 N1 N2 1 a 2 6 2 b 4 NA 3 d 6 0 4 F 9 6 > is.na(x) S1 N1 N2 1 FALSE FALSE FALSE 2 FALSE FALSE TRUE 3 FALSE FALSE FALSE 4 FALSE FALSE FALSE No I want to be able to do: 1) Setting

Hello, my group is investigating the development of a statistical program in support of some other work we are doing. Stats is not really our area of expertise and we would like to leverage existing code where possible. From your site, I have gathered that R is pretty capable and well used. We are currently developing in python and targetting a primarily windows environment. We would like to

Hi Group, I wonder if R has the function CensorReg which takes care of left truncation. I know it's in S. Thank you so much! Kel __________________________________________________ Do You Yahoo!? Get personalized email addresses from Yahoo! Mail - only $35 a year! http://personal.mail.yahoo.com/ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help

Hello, I'm trying to plot simple data containing numbers/year. I use plot(numbers~year, type="l"), but the problem is that there is no data for several years. S I would like the line to be interrupted for the period without data, and then continue afterwards. How can I do this ? Also, I've been using par(new=TRUE) to put two curves (same years, different numbers) on the same

Hi, I am almost new to R. Would any one show me (please include an example) how to clip a contour or image figure so that I only see the image or lines inside a selected area? Thank you very much Han -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

Hi, I input this command on a data matrix: > apply( traffic, 2, function(z){z-mean(z)}) and got the following error messages: Error in sum(..., na.rm = na.rm) : invalid "mode" of argument Can anyone tell me what is wrong here? Running R1.2.2 on Windows ME. Cheers, Kevin ------------------------------------------------- Ko-Kang Kevin Wang Statistical Analysis Division

Hello all, to obtain a list of all objects in all search paths, I've found the following to work: > biglist <- sapply(1:length(search()), objects) This more obvious one, however, does not work: > biglist <- sapply(search(), objects) Error in pos.to.env(pos) : invalid argument Still, search() gives [1] ".GlobalEnv" "package:ctest" "Autoloads"

Everyone, I have an array of integer values which are located on a uniform 2-D grid. I want to plot the integer values at the node locations. The closest I have come is with the following code: for (i in 1:ny) { for (j in 1:nx) { ncell <- nx*(i - 1) + j ch <- as.character(ncount[ncell]) tx <- j*dx - dx/2 ty <- ny*dy - (i-1)*dy - dy/2 points(tx, ty, pch = ch, col =

Dear r-plus users, i would like to use outer in the following case outer(x,y,FUN="fun") i suppose that my function fun is of the following form: fun<-function(x,y) { if(y>x) return(x+y) if(y<=x) return(0) } My problem is that the command outer(x,y,FUN="fun") return me a null matrix instead of an upper triangular matrix. Is somebody have a solution ? Thanks for your

Hi, I've got a matrix of two columns, say: > d <- t(array(1:20,dim=c(2,10))) > d [,1] [,2] [1,] 1 2 [2,] 3 4 [3,] 5 6 [4,] 7 8 [5,] 9 10 [6,] 11 12 [7,] 13 14 [8,] 15 16 [9,] 17 18 [10,] 19 20 now I need those values from column #2 where column #1 is equal to certain values. eg.: > i <- c(13,5,19) The

Hello, first, I'd like to congratulate the core team to the new R version 1.1. I think it's a great update, with glimpses into an even greater future (tcltk!). While playing around with the new functions (on Win 95), I found the following: As the documentation states, when I generate a plot and save it with recordPlot, I can regenerate it by printing the variable: > plot(1:10) >

Currently (R-1.4.1 as well as R-devel, according to http://stat.ethz.ch/R-alpha/R-devel/library/base/html/par.html), the html version of help(par) shows [...] lab A numerical vector of the form c(x, y, len) which modifies the way that axes are annotated. The values of x and y give the (approximate) number of tickmarks on the x and y axes and len specifies the label size. The default

Hello, does anyone know of a quick way to create a data frame "like" another, but with more rows? What I'd like to do is this: if mydata is a data.frame like a b c 1 TRUE yes 2 FALSE no 3 TRUE yes I'd like to get mydata2 with the same column names and column types, but without the values and with more rows. All I could think of was to manually do

Hello, I'd like to propose an extension to the function summary.aov. In Splus (2000, I don't know about other versions), summary.aov allows a parameter ssType to be set to 1 or 3 (defaults to 1) to choose the type of Sums of Squares. I know I can get Type III SS in R with drop1(model), but including the functionality into summary.aov would, in my opinion, - yield a more usable table

This is only partly an R question, but this seems a good place to ask. Like many, if not most, statistics programs, R has no post hoc tests for Kruskal Wallis or Friedman, and as far as I can make out they are not available in any of the R packages. I would be grateful if someone could point me in the right direction if I am wrong. As I said this is not uncommon. Parametric Anova procedures

On Thu, Jun 22, 2000 at 05:01:32PM +0200, p.dalgaard@biostat.ku.dk wrote: > I think this is pretty certainly a bug, so I'm cc'ing this to r-bugs. > > -p > > Kaspar Pflugshaupt <pflugshaupt@geobot.umnw.ethz.ch> writes: > > As the documentation states, when I generate a plot and save it with > > recordPlot, I can regenerate it by printing the variable:

On Thu, Jun 22, 2000 at 05:01:32PM +0200, p.dalgaard@biostat.ku.dk wrote: > I think this is pretty certainly a bug, so I'm cc'ing this to r-bugs. > > -p > > Kaspar Pflugshaupt <pflugshaupt@geobot.umnw.ethz.ch> writes: > > As the documentation states, when I generate a plot and save it with > > recordPlot, I can regenerate it by printing the variable:

Hello, is there a quick way of estimating Weibull parameters for some data points that are assumed to be Weibull-distributed? I guess I'm just too lazy to set up a Maximum-Likelihood estimation... ...but maybe there is a simpler way? Thanks for any hint (and yes, I've read help(Weibull) ;) Kaspar Pflugshaupt -- Kaspar Pflugshaupt Geobotanical Institute ETH Zurich, Switzerland

Hello, I've fitted a parametric survival model by > survreg(Surv(Week, Cens) ~ C(Treatment, srmod.contr), > data = poll.surv.wo3) where srmod.contr is the following matrix of contrasts: prep auto poll self home [1,] 1 1 1.0000000 0.0 0 [2,] -1 0 0.0000000 0.0 0 [3,] 0 -1 0.0000000 0.0 0 [4,] 0 0 -0.3333333 1.0 0 [5,] 0 0