[R] Help is sought for a small problem

Dear there, I tried today within my R 1.7.0 under Windows XP the command update.packages() and then I answered y in all subsequent questions. From the software responses I understood that everythin went OK from VR to Rcmdr, but in survival it said Error in file(filem r) : unable to open

Re: [R] Help is sought for a small problem

Dipl. Ing. Kurt Neumann wrote: Dear there, I tried today within my R 1.7.0 under Windows XP the command update.packages() and then I answered y in all subsequent questions. From the software responses I understood that everythin went OK from VR to Rcmdr, but in survival it said Are you sure

[R] rank with ties

Hi, Is there a function like rank but that solves the ties by randomly assigning a value (doesn't average ranks of ties). This is what I actually need: I want to make NA all elements of each column in an array that are ranked in a position larger that rankmax for each column. # Say I've got an

[R] curves with shaded areas?

Hi, I want to make a plot with abline where the area below or above the curve is shaded. I can't find any documentation on that. Can anybody help me with that? thanks, Juffi __ [EMAIL PROTECTED] mailing list

[R] barplot

is this the right mailinglist for my trivial barplot() problem? mfg wolfgang __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] curves with shaded areas?

See ?polygon, ?rect for how to shade areas. If you make a plot with abline, surely it has no curves (only straight lines)? On Tue, 22 Jul 2003, Johannes Fuernkranz wrote: I want to make a plot with abline where the area below or above the curve is shaded. I can't find any documentation on

RE: [R] Confidence Band for empirical distribution function

Dear Leif, If you look at the definition of ks.test, you'll find the lines pkstwo - function(x, tol = 1e-06) { if (is.numeric(x)) x - as.vector(x) else stop(Argument x must be numeric) p - rep(0, length(x)) p[is.na(x)] - NA IND - which(!is.na(x) (x 0)) if

RE: [R] curves with shaded areas?

polygon() together with par(usr) would be your friend here. As I do not know if you want horisontal/vertical lines or any line, here is a general example: # Generate the data x - seq(from=1, to=3.5*pi, length=100) y - sin(x) # Creates an empty plot of right size plot(x,y, type=n) # Get

RE: [R] curves with shaded areas?

De : [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] la part de Johannes Fuernkranz [...] I want to make a plot with abline where the area below or above the curve is shaded. I can't find any documentation on that. Can anybody help me with that? [...] You have to play with par(usr). See the

Re: [R] curves with shaded areas?

Prof Brian Ripley wrote: See ?polygon, ?rect for how to shade areas. Thank you. This covers most but not all of my problems. :-) If you make a plot with abline, surely it has no curves (only straight lines)? Right, sorry. I was a bit over-simplifying my problem, because I thought this can be

Re: [R] curves with shaded areas?

On Tue, 22 Jul 2003, Johannes Fuernkranz wrote: Prof Brian Ripley wrote: See ?polygon, ?rect for how to shade areas. Thank you. This covers most but not all of my problems. :-) If you make a plot with abline, surely it has no curves (only straight lines)? Right, sorry. I was a

RE: [R] R and C++ compared with only C++

I must say that I agree with Andrew. I recently did just that: took an R package with C code and turn it into a stand-alone C program. Granted I'm far from being a C expert, but not only that the code didn't run any faster (because all the heavy computations were already done in C anyway), but

Re: [R] Calling R from C

Duncan Temple Lang's article In Search of C/C++ Fortran Routines in R News Vol 1/3, September 2001 is definitely worth reading (even though it's mainly about calling C from R): http://cran.r-project.org/doc/Rnews/Rnews_2001-3.pdf and you should have a look at Writing R Extensions :

RE: [R] generate a series of fucntion

You need to be careful if you take the approach below. Try it: g1 - list(function(x) log(x), function(x) log(x+1), function(x) x, function(x) x^2) g2 - lapply(g1, function(g) { + function(x,t) exp(-t[1]-t[2]*g(x)-t[3]*g(1+x)) + } + ) lapply(g2, function(f) get(g,

Re: Rcmdr problems (was:Re: [R] Minor nuisance with rw1071)

Dear Kjetil, At 11:04 PM 7/21/2003 -0400, kjetil brinchmann halvorsen wrote: On 18 Jul 2003 at 9:54, John Fox wrote: Sorry for being so late coming back on this, I was travelling the weekend. I tried to remove as much as possible from the workspace (.RData file), but nothing created by Rcmdr.

[R] variable names

Dear helpers I want to use rpart several times in a loop to build a classification tree. My problem is that rpart needs a formula as argument and for that the variables need to have names and this doesn't happen in my case. Every iteration in the loop has a different dataset with several

[R] Conditional Statements for Graphing

Dear List I have math test scores for male and female students where gender is a dummy code (female =1). I also have a variety of other demographic variables. However to begin, I want to create a very simple stripchart where female math scores are a blue circle and male scores are a red

Re: [R] variable names

On Tue, 22 Jul 2003, Luis Miguel Almeida da Silva wrote: Dear helpers I want to use rpart several times in a loop to build a classification tree. My problem is that rpart needs a formula as argument and for that the variables need to have names and this doesn't happen in my case. Every

RE: [R] variable names

Dear Luis, You might want to have a look at Bill Venables. Programmer's niche. R News, 2(2):24-26, June 2002 which you can find at http://cran.r-project.org/doc/Rnews/ or look into the manual R Language Definition, chapter Computing on the language. Assuming that in your case the variable to be

RE: [R] variable names

I didn't noticed that fact. I've already found a way to do that x - 1:40 colnames(df.treino) - paste(Ncp,x,sep=.) and this generates names that I can relate with the variables. Thanks anyway The problem is that I use rpart in a loop and the class labels are in the last column. For the above

RE: [R] variable names

On Tue, 22 Jul 2003, Luis Miguel Almeida da Silva wrote: I didn't noticed that fact. I've already found a way to do that x - 1:40 colnames(df.treino) - paste(Ncp,x,sep=.) and this generates names that I can relate with the variables. Thanks anyway The problem is that I use rpart in

RE: [R] variable names

Read Bill Venables' column in R News issue 2/2, page 24. Andy From: Luis Miguel Almeida da Silva [mailto:[EMAIL PROTECTED] I didn't noticed that fact. I've already found a way to do that x - 1:40 colnames(df.treino) - paste(Ncp,x,sep=.) and this generates names that I can relate

[R] greek in main title

Hello, I have written a function that demonstrates the CLT by generating samples following the exponential distribution, calculating the means, plotting the histogram, and drawing the limiting normal curve as an overlay. I have the title of each histogram state the sample size and rate (1/theta)

[R] Making a group membership matrix

Hi Helpers: I have a factor object that has 314k entries of 39 land cover types. (This object can be coerced to characters neatly should that be easier to work with.) length(foo) [1] 314482 foo[1:10] [1] Montane Chaparral BarrenRed Fir Red Fir [5] Red Fir Red

Re: [R] greek in main title

Hi, change the name of function argument theta to something else say thta and then use the following heading - expression(paste(Exp w/ , theta, =, thta, , n = , n[i])) Hope it will help... Mahbub. --- [EMAIL PROTECTED] wrote: Hello, I have written a function that demonstrates the

Re: [R] greek in main title

I was not able to get theta italicized, but I otherwise got something matching my understanding of what you wanted. Consider the following: plotExpNormal - function (theta) { layout(matrix(c(1:4),2,2,byrow=TRUE)) par(bg = cornsilk) n - c(1,10,25,50) for(i in 1:4) { xbar - rep(0,

Re: [R] Making a group membership matrix

Hi, The resulting matrix will be *really* large, so make sure you have enough RAM. You might reduce this by dropping all the levels of the factor that have zero counts. As far as the solution, it is a one-liner, but not really a crypic one: model.matrix(~foo-1) (assuming that you haven't

[R] animal models and lme

-- Hi R users, I'm trying to fit an animal model using lme and cannot fit the corelation matrix which describes the additive genetic correlation between individuals. I have created a test data set (pedigree) where the first column specifies the id of each individual, and the second column

Re: [R] greek in main title

Hi -- your solution looses the ability to assign the values of theta and n[i] in the simulation... you get thta and n_i in the titles. In a message dated 7/22/2003 10:44:59 AM Eastern Standard Time, [EMAIL PROTECTED] writes: Hi, change the name of function argument theta to something else

Re: [R] animal models and lme

On Wed, 23 Jul 2003, Jarrod Hadfield wrote: I can't get your message included, for some reason. Your lme model has no random effects: that is not allowed. Either include a random effect or (more likely) use gls or lm.gls (MASS) if you really just want a GLS fit. -- Brian D. Ripley,

Re: [R] greek in main title

Here is another solution. I'm not sure what you are trying to accomplish with font.main=6, though. See ?par for details: have a closer look at the description for font option. Cheers, Jerome theta - 2 layout(matrix(c(1:4),2,2,byrow=TRUE)) par(bg = cornsilk) n - c(1,10,25,50) for(i

Re: [R] plot character

On Tue, 2003-07-22 at 11:59, Luis Miguel Almeida da Silva wrote: Dear helpers Is it possible to plot with a desired character? For example tr for training error and te for test error. I tried plot(x,y,pch=tr) but only appears t in the plot One option is to use text(x, y,

Re: [R] plot character

Have you considered ?text? hope this helps. spencer graves Luis Miguel Almeida da Silva wrote: Dear helpers Is it possible to plot with a desired character? For example tr for training error and te for test error. I tried plot(x,y,pch=tr) but only appears t in the plot

Re: [R] plot character

On Tue, 22 Jul 2003 17:59:24 +0100, Luis Miguel Almeida da Silva [EMAIL PROTECTED] wrote : Dear helpers Is it possible to plot with a desired character? For example tr for training error and te for test error. I tried plot(x,y,pch=tr) but only appears t in the plot Use text() instead:

[R] Matrix manipulation - dropping rows based on a query

I have a matrix with some very funky data. I want to drop the columns with a mean 1 (for instance). #Example: my.mat - matrix(rnorm(100, 0, 3), nrow = 10) #I know I can get the row means: apply(my.mat,1, mean) I think I need to get a vector that has those rows 1 my.mat[-rows1,] The actual

Re: [R] Matrix manipulation - dropping rows based on a query

tst - apply(my.mat, 2, mean) 1 my.mat[,tst] Will return only columns with mean 1 Jim James W. MacDonald Affymetrix and cDNA Microarray Core University of Michigan Cancer Center 1500 E. Medical Center Drive 7410 CCGC Ann Arbor MI 48109 734-647-5623 [EMAIL PROTECTED] 07/22/03 01:53PM I

[R] libblas.so.3

I am using a Pentium 4 and R 1.7.1 I removed libblas.so.3 to make sure that my programs are using ATLAS ( maching dependent implementation of BLAS ). R is not working. Any suggestion? __ [EMAIL PROTECTED] mailing list

Re: [R] Matrix manipulation - dropping rows based on a query

James MacDonald wrote: tst - apply(my.mat, 2, mean) 1 ... and rowMeans(my.mat) 1 will be faster, but perhaps not working for the requested more complicated query. Uwe Ligges my.mat[,tst] Will return only columns with mean 1 Jim James W. MacDonald Affymetrix and cDNA Microarray Core

Re: [R] libblas.so.3

Francisco J Molina wrote: I am using a Pentium 4 and R 1.7.1 I removed libblas.so.3 to make sure that my programs are using ATLAS ( maching dependent implementation of BLAS ). R is not working. Any suggestion? First suggestion: Please tell us some details, in particular why you think R is not

Re: [R] libblas.so.3

I am using linux rehat 9 When I try to run R I get usr/lib/R/bin/R.bin: error while loading shared libraries: libblas.so.3: cannot open shared object file: No such file or directory ( because I removed it ) I thought of compiling R from source but I have read: R currently uses only level 1

Re: [R] libblas.so.3

Francisco J Molina wrote: I am using linux rehat 9 When I try to run R I get usr/lib/R/bin/R.bin: error while loading shared libraries: libblas.so.3: cannot open shared object file: No such file or directory ( because I removed it ) I thought of compiling R from source but I have read: R

[R] Curious warning in R for OS X w/Xwindows

I recently reformatted my TiBook and reinstalled everything from the system up. After installing Apple's version of XWindows and the XWindows version of R 1.7.0, I ran the two and got a strange warning: - R : Copyright 2003, The R Development Core Team Version 1.7.0 Patched (2003-05-14) R

Re: [R] Curious warning in R for OS X w/Xwindows

Bear F. Braumoeller [EMAIL PROTECTED] writes: I recently reformatted my TiBook and reinstalled everything from the system up. After installing Apple's version of XWindows and the XWindows version of R 1.7.0, I ran the two and got a strange warning: - R : Copyright 2003, The R

Re: [R] Curious warning in R for OS X w/Xwindows

On Tuesday, July 22, 2003, at 06:24 PM, Peter Dalgaard BSA wrote: OS X is not quite Unix/Linux, so I don't really know, but I'd suspect that somehow the environment variable that describes the terminal capabilities (usually called TERM) got unset or not set properly for the terminal emulator in

[R] Processing a large number of files

I maintain the Devore5 package which contains the data sets from the 5th edition of Jay Devore's text Probability and Statistics for Engineering and the Sciences. The 6th edition has now been published and it includes several new data sets in exercises and examples. In addition, some exercises

Re: [R] Boosting, bagging and bumping. Questions about R tools andpredictions.

Hi, If you want to learn the theory of boosting and bagging, and other classification techniques, then you will want to refer to Hastie, Tibshirani and Friedman's The Element of Statistical Learning: Data Mining, Inference, and Prediction by Springer. It is the best book I have seen in these

[R] Boosting,bagging and bumping. Questions about R tools and predictions.

I'm interested in further understanding the differences in using many classification trees to improve classification rates. I'm also interested in finding out what I can do in R and which methods will allow prediction. Can anybody point me to a citation or discussion? Specifically, I want to

RE: [R] Boosting,bagging and bumping. Questions about R tools and predictions.

http://www.boosting.org/publications.html I found some of the papers on this page useful in understanding the concepts you refer to. I will leave it to the better informed members of the group to talk about the packages that relate to this field. -Original Message- From: [EMAIL

[R] pls.pcr compared to Unscrambler

Dear R-helpers, Has anybody ever tried to compare pls regression outputs from the pls.pcr R-package developped by Wehrens with outputs from the Unscrambler software developped by CAMO company ? I find very different outputs and wonder if this comes from differences between methods/algorithms

RE: [R] animal models and lme

Hi, You should look at Pinheiro and Bates (2000) Mixed-effects models in S and S-Plus. It describes how to format the correlation matrix to pass to functions lme and gls. Basically, the correlation matrix has to be one of the corStruct classes, probably corSymm for your example. So in the call

[R] using getBioC()

Hello, I am trying to install R/Bioconductor on a G4 Mac running OS X. I have successfully installed R so that a command window opens, but installation of the downloaded Bioconductor package is giving me trouble. After copying/pasting the Bioconductor installation script in to the window

[R] Package submission---HyperbolicDist

I have just uploaded a package for the hyperbolic distribution to CRAN. From the help file: This library provides a collection of functions for working with the hyperbolic distribution. Included are the density function, distribution function, quantiles, random number generation and fitting