Re: [R] Multiple colors in a plot title

Leo Espindle wrote: R-Help: Is there a way to use multiple colors in the title of a plot? For instance, to have certain words be red, and certain words be blue? No, you have top typset them word by word. Uwe Ligges thanks in advance, Leo __

Re: [R] Matrix multiplication in R is inaccurate!!

Vikas Rawal wrote: If you multiply a matrix by its inverse you should get an identity matrix. In R, you get an answer that is accurate up to about 16 decimal points? Why can't one get a perfect answer? Because R makes use of a digital computer which cannot work exactly using floating-point

Re: [R] self-written function

Hi! I´m using the function for a data frame, but up to now it only works with single vectors, such as backsin(variable,grouping.factor) Actually, you´re right, the ... in the function definition is not needed...:-) Regards Christoph Erin Hodgess wrote: Hello Christoph! I have a question about

Re: [R] Multiple colors in a plot title

Try something like this: plot(1:10) mtext(c(title1, title2), col = c(2:3), at = c(4, 6)) Leo Espindle wrote: R-Help: Is there a way to use multiple colors in the title of a plot? For instance, to have certain words be red, and certain words be blue? thanks in advance, Leo

Re: [R] clustering

Did you try pam() from package cluster? I'd recommend it as an improvement compared to kmeans(). Martin Maechler, ETH Zurich __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

RE: [R] R-Help : running MIX package

On 28-Jan-05 Jeanhee Hong wrote: Hello all. I am inexperienced with R and am clumsily trying to work through it for specific multiple imputations Id like to run for my thesis.In running the MIX package, I keep getting an error message regarding the use of the prelim.mix command. Error in

[R] R-Help : running MIX package

Hello all. I am inexperienced with R and am clumsily trying to work through it for specific multiple imputations Id like to run for my thesis.In running the MIX package, I keep getting an error message regarding the use of the prelim.mix command. Error in as.integer.default(list(alcohol =

Re: [R] self-written function

Hi! OK, here are some more details on the function: My dataframe consists of several columns of categorical variables (lets call them A,B,C) plus a column with a response variable y (which is arcsine-square root transformed proportions) I am now trying to write a function that automatically

[R] GLM fitting

DeaR R-useRs, I'm trying to fit a logist model with these data: dati y x 1 1 37 2 1 35 3 1 33 4 1 40 5 1 45 6 1 41 7 1 42 8 0 20 9 0 21 10 0 25 11 0 27 12 0 29 13 0 18 I use glm(), having this output: g-glm(y~x,family=binomial,data=dati) Warning messages: 1: Algorithm did not

[R] Conflicts using Rcmdr, nlme and lme4

Hello all! R2.0.1, W2k. All packages updated. I´m heavily dependant on using mixed models. Up til´now I have used lme() from nlme as I have been told to. Together with estimable() from gmodels it works smooth. I also often run Rcmdr, mostly for quick graphics. After using Rcmdr, on reopening

[R] plot

hello, I wonder whether there is a way to plot in R the following: I have two arrays of dimension (200*1) say x1 and x2. I would like to plot on the same graph x1[1:50,]*x2[1:50,] with pch=1 x1[51:100,]*x2[51:100,] with pch=2

Re: [R] Finding runs of TRUE in binary vector

Sean Davis wrote: I have a binary vector and I want to find all regions of that vector that are runs of TRUE (or FALSE). a - rnorm(10) b - a0.5 b [1] TRUE TRUE TRUE FALSE TRUE FALSE TRUE TRUE TRUE TRUE My function would return something like a list: region[[1]] 1,3 region[[2]] 5,5

RE: [R] GLM fitting

You have a perfect separtaion of y by x, i.e. y == (x30) is true for all units. Bendix Carstensen -- Bendix Carstensen Senior Statistician Steno Diabetes Center Niels Steensens Vej 2 DK-2820 Gentofte Denmark tel: +45 44 43 87 38 mob: +45 30 75 87 38 fax: +45 44 43 07 06

[R] latent class model in R

Hi Everyone, I am facing the problem of running the latent class model in R. By mmlcr, the program always shows as following, Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : invalid variable type Does anyone happen to know why this happened? I am

Re: [R] clustering

Hi, EMclust in package mclust fits normal mixtures. Note that if you split your data values into intervals, the resulting distributions conditional on the intervals are not normals, but truncated normals! This is important if you try to check within group normality, unless you have strongly

Re: [R] self-written function

Dear Jari, I had included the test for NA´s just to be sure the function can handle datasets with missing values. Here´s an updated version without this specification: backsin-function(x,y){ backtransf-list() back-((sin(x))^2)*100 backtransf$mback-tapply(back,y,mean)

Re: [R] plot

try this: x1 - as.matrix(rnorm(200)) x2 - as.matrix(rnorm(200)) ### dim(x1) - dim(x2) - c(50,4) matplot(x1,x2, pch=1:4, col=1) Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35,

Re: [R] Conflicts using Rcmdr, nlme and lme4

CG Pettersson wrote: Hello all! R2.0.1, W2k. All packages updated. I´m heavily dependant on using mixed models. Up til´now I have used lme() from nlme as I have been told to. Together with estimable() from gmodels it works smooth. I also often run Rcmdr, mostly for quick graphics. After using

Re: [R] plot

NICOLAS DEIG wrote: hello, I wonder whether there is a way to plot in R the following: I have two arrays of dimension (200*1) say x1 and x2. I would like to plot on the same graph x1[1:50,]*x2[1:50,] with pch=1

RE: [R] plot

if * in x1[1:50,]*x2[1:50,] means times, try : matplot(matrix(x1,ncol=4)*matrix(x2,ncol=4)) if * means vs, try : matplot(matrix(x1,ncol=4),matrix(x2,ncol=4)) -Message d'origine- De : [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] la part de NICOLAS DEIG Envoye : vendredi 28 janvier 2005

[R] Windows Batch File

Hiya, I have read the FAQ for windows and read point 2.10. I setup a dos batch file which contains the following line: C:\progra~1\R\rw2001\bin\Rterm.exe --vanilla I then tried: R CMD BATCH --argument myscript.R result.out then it started R. I thought by running R CMD BATCH executes R

Re: [R] plot

Hello, You search that command : plot(x1,x2,pch=rep(1:4,50)) I think you should read some basic manuals about R, the FAQ, etc ... NICOLAS DEIG wrote: hello, I wonder whether there is a way to plot in R the following: I have two arrays of dimension (200*1) say x1 and x2. I would like to plot on

Re: [R] plot

Hello, You search that command : plot(x1,x2,pch=rep(1:4,50)) Oups, sorry, I meant that : plot(x1,x2,pch=rep(1:4,each=50)) I need more cafeine on morning ! I think you should read some basic manuals about R, the FAQ, etc ... NICOLAS DEIG wrote: hello, I wonder whether

Re: [R] self-written function

Hi Christoph On 28 Jan 2005 at 11:04, Christoph Scherber wrote: Hi! OK, here are some more details on the function: My dataframe consists of several columns of categorical variables (lets call them A,B,C) plus a column with a response variable y (which is arcsine-square root transformed

RE: [R] R-Help : running MIX package

On 28-Jan-05 Ted Harding wrote: On 28-Jan-05 Jeanhee Hong wrote: Hello all. I am inexperienced with R and am clumsily trying to work through it for specific multiple imputations Id like to run for my thesis.In running the MIX package, I keep getting an error message regarding the use of the

[R] xfig

hello all, I know I can export a graphic from R to Xfig from the manuals, but I just cannot figure out it works. does anyone know have to do so? thanks in advance. Nicolas __ R-help@stat.math.ethz.ch mailing list

Re: [R] Windows Batch File

Elaine Han wrote: Hiya, I have read the FAQ for windows and read point 2.10. I setup a dos batch file which contains the following line: C:\progra~1\R\rw2001\bin\Rterm.exe --vanilla This line fires up R, but won't accept any file as input... The FAQ tells you to use path_to_R\bin\Rterm.exe

RE: [R] Conflicts using Rcmdr, nlme and lme4

Dear CG, Doug Bates has already addressed the conflict between the nlme and lme4 packages. The Rcmdr package loads a number of other packages to use functions in them. I frankly don't recall why I need nlme, but something in it must be used somewhere in the Rcmdr. You could get rid of the

RE: [R] Conflicts using Rcmdr, nlme and lme4

Dear CG, An addendum to my previous response: If you do decide to recompile the Rcmdr package to eliminate the dependency on nlme, as I suggested, you'll also have to edit the .onLoad() function in the package (which is in the file startup.R in the source package) to remove nlme from the vector

[R] Double integration

__ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] xfig

NICOLAS DEIG wrote: hello all, I know I can export a graphic from R to Xfig from the manuals, but I just cannot figure out it works. does anyone know have to do so? thanks in advance. Nicolas What about reading ?xfig Use it like any other device: xfig(filname) plot(1:10) dev.off() Uwe Ligges

RE: [R] Error: cannot allocate vector of size... but with a twist

Just a couple of remarks below... From: James Muller Hi, I have a memory problem, one which I've seen pop up in the list a few times, but which seems to be a little different. It is the Error: cannot allocate vector of size x problem. I'm running R2.0 on RH9. My R program is

[R] SparseLogReg make r crash problem

We are currently using runSparseLogreg in Package SparseLogReg, here is our code : library(SparseLogReg) randomdata - matrix(runif(1500,min=0,max=3),ncol=10) class.zero.one-c(rep(0,7),rep(1,3)) # give the ramdom gene expression samples with class label

[R] Begginer with R

Hello all, Im just beggining using R. I have read a couple of introductory documents but they are very general. Is there a document focused on classification? (this is what Ill be using R for). Examples would be just fine. Thanks in advance, Ulises __

RE: [R] GLM fitting

Vito: Please plot your data: y - c(1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0) x - c(37, 35, 33, 40, 45, 41, 42, 20, 21, 25, 27, 29, 18) plot(x, y) You will see that ANY step function between 29 x 33 will describe these observations perfectly. Charles Annis, P.E. [EMAIL PROTECTED] phone:

[R] read.matrix.csr bug (e1071)?

Hello, I would like to read and write sparse matrices using the functions write.matrix.csr() and read.matrix.csr() of the package e1071. Writing is OK but reading back the matrix fails: x - rnorm(100) m - matrix(x, 10) m[m 0.5] - 0 m.csr - as.matrix.csr(m)

Re: [R] read.matrix.csr bug (e1071)?

Don't you want read.matrix.csr not read.matrix? url:www.econ.uiuc.edu/~rogerRoger Koenker email [EMAIL PROTECTED] Department of Economics vox:217-333-4558University of Illinois fax:217-244-6678

RE: [R] GLM fitting

To push the point a bit further, Vito, if you allow just a bit variation in the data by changing one of the y=1 cases to 0 or one of the y=0 cases to 1, then you'll be able to fit the glm model. If these are real-world data and if you still want to describe them, then a deterministic statement

Re: [R] read.matrix.csr bug (e1071)?

On Fri, 28 Jan 2005, roger koenker wrote: Don't you want read.matrix.csr not read.matrix? Sorry for the mistake, the error is caused by the line read.matrix.csr(~/data/sparse.dat) __ R-help@stat.math.ethz.ch mailing list

Re: [R] Conflicts using Rcmdr, nlme and lme4

Douglas Bates wrote: I can guarantee that this is the last time that I will fundamentally redesign the computational methods and data structures used in fitting mixed-effects models. This is the fifth time I have done it from scratch and after this one I am quitting the business. Hope that

Re: [R] SparseLogReg make r crash problem

[EMAIL PROTECTED] wrote: We are currently using runSparseLogreg in Package SparseLogReg, here is our code : library(SparseLogReg) randomdata - matrix(runif(1500,min=0,max=3),ncol=10) class.zero.one-c(rep(0,7),rep(1,3)) # give the ramdom gene expression samples with class label

Re: [R] clustering

Hi, I think truncated normality is what I meant in the last email. The experiments (, which might not be representative) show k-mean and EM gave me comparable results, while EM is a little bit better. (I used Weka for this purpose). I will try them on the real data and also try pam() and if I have

[R] JOB: Software Architect, Insightful Seattle

Insightful is searching for a candidate to fill a technology leadership position whose primary function is to drive the Insightful product and solutions architectures. Utilize your previous experience as you work closely with a seasoned executive team and a world-class development team to help

Re: [R] Conflicts using Rcmdr, nlme and lme4

Kjetil Brinchmann Halvorsen wrote: Douglas Bates wrote: I can guarantee that this is the last time that I will fundamentally redesign the computational methods and data structures used in fitting mixed-effects models. This is the fifth time I have done it from scratch and after this one I am

[R] Line types in the legend function

Hi, I'm using version 2.0.1 of R on a Windows 2000 platform. The legend function has parameters _lty_ and _pch_ for drawing line types and point types in the legend box, but I can't find any way of getting the patterns corresponding to _type='b'_ in the _plot_ function. When you enter _type='b'_

Re: [R] How to generate labels or names?

Dear Eric, try something like N - 20 paste(Lab,1:N,sep=) Regards, Samuel. Eric Rodriguez [EMAIL PROTECTED] wrote: Hi, I'm new to R and I would like to generate labels like data.frame does : V1 V2 V3 I'm trying to generate a N vector with label such as Lab1 Lab2 ... LabN. I guess

[R] extracting from a data.frame

Hi, I am sorry for this simple question, but... How do I extract something from a data.frame? The following is my Problem: I have got a dataframe a with various columns. One of those columns is called V3 and contains elements of the following levels: levels(a$V3) [1] C CA CB CD CD1 CD2

Re: [R] extracting from a data.frame

On Fri, 28 Jan 2005 19:26:14 +0100, dax42 [EMAIL PROTECTED] wrote : Hi, I am sorry for this simple question, but... How do I extract something from a data.frame? The following is my Problem: I have got a dataframe a with various columns. One of those columns is called V3 and contains elements

Re: [R] extracting from a data.frame

Look at ?subset. Sean On Jan 28, 2005, at 1:26 PM, dax42 wrote: Hi, I am sorry for this simple question, but... How do I extract something from a data.frame? The following is my Problem: I have got a dataframe a with various columns. One of those columns is called V3 and contains elements of the

Re: [R] extracting from a data.frame

newdat - subset(a, a$V3 == O) ?subset dax42 wrote: Hi, I am sorry for this simple question, but... How do I extract something from a data.frame? The following is my Problem: I have got a dataframe a with various columns. One of those columns is called V3 and contains elements of the following

Re: [R] Begginer with R

http://www.liacc.up.pt/~ltorgo/DataMiningWithR/ might be a good start point. also the book: MASS (4th Ed) Check some previous lists and you will find some similar topics too. I used R just for like 2 weeks. Not sure of your level of knowledge in data mining. but hope this help. Ed On Fri,

Re: [R] Installing Problems

Hi Ashok, On Thu, Jan 27, 2005 at 11:48:08AM -0500, Ashok Veeraraghavan wrote: I tried installing R on my MAC OS 10.3. After R installation I tried installing BioConductor which requires R. I ran into some problems with Bioconductor. Right now I want to remove (uninstall) all R and

[R] using RODBC

I am trying to bring data into R from an excel spreadsheet in order to perform several statistical tests on it. I was trying to use odbcConnectExcel in the RODBC package. Once I am connected to the excel file, how do I select rows and columns from the file in order to analysis them in R.

[R] Exponential Fits to Distribution Tails

EPA has suggested an exponential fit to the upper 10 percent of a PM10 distribution and using that fit to determine a once per year frequency of occurrence. My question is one of pedigree -- does such a technique have merit and status in the statistics community? Or is there a better technique

[R] Direct sum of matrices

Hi R users: How can I built a direct sum function of matrices in R? I mean A(mxn), B(pxq), C(rxs),... X-ds(A,B,C,...) X = [ A, 0, 0 0, B, 0 0, 0, C] ((m+p+r+...) x (n+q+s+...)) Thank you for your help. Kenneth -- Using Opera's revolutionary e-mail client:

RE: [R] using RODBC

Have you consulted the R Data Import/Export Manual, available from the Rgui via the Help tab? Mike -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jessica Higgs Sent: 28 January 2005 19:24 To: R-help@stat.math.ethz.ch Subject: [R] using RODBC I am

Re: [R] Direct sum of matrices

?bdiag from package:assist may help you a bit. HTH, Danny On Fri, 28 Jan 2005, Kenneth wrote: Hi R users: How can I built a direct sum function of matrices in R? I mean A(mxn), B(pxq), C(rxs),... X-ds(A,B,C,...) X = [ A, 0, 0 0, B, 0 0, 0, C]

RE: [R] using RODBC

Would this function be useful to you. This imports the data from an excel file and gives it to you as a data frame. http://www.maths.lth.se/help/R/.R/library/gdata/html/read.xls.html I have no data yet. It is a capital mistake to theorize before one has data. Insensibly one begins to twist facts

[R] R for CGI

Dear R Users; Perl is the common language to write CGI scripts which handle Forms. My question is that can R be as fast as perl to do the same job(with using CGIwithR package). Is it an optimal solution to connect R directly to a commercial HTML webpages, Sincerely, Sean

Re: [R] Direct sum of matrices

Kenneth krcabrer at unalmed.edu.co writes: : : Hi R users: : : How can I built a direct sum function of matrices in R? : : I mean A(mxn), B(pxq), C(rxs),... : : X-ds(A,B,C,...) : : X = [ A, 0, 0 : 0, B, 0 :0, 0, C] ((m+p+r+...) x (n+q+s+...)) : : Thank you for

[R] Why am I not getting some r-help posts?

At the risk of being flamed for mortal stupidity, I'm trying to figure out what could have possibly changed (at my end almost certainly) that make Brian Ripley's posts (in particular) not show up on my mail server (at least not on the one where I receive my email). All my email comes to a

[R] avoiding loops

Hi again, thanks a lot for the quick answer. I just forgot the comma, always these stupid mistakes... Anyways, as I said before, I have two data.frames containing about 1000 rows and I would like to avoid looping through all of them... In each data.frame are coordinates (x,y,z), so every row

RE: [R] avoiding loops

Doesn't this work? d1 - data.frame(x=rnorm(100), y=rnorm(100), z=rnorm(100)) d2 - data.frame(x=rnorm(100), y=rnorm(100), z=rnorm(100)) dist.12 - sqrt(rowSums((d1 - d2)^2)) Andy From: dax42 Hi again, thanks a lot for the quick answer. I just forgot the comma, always these stupid

Re: [R] Why am I not getting some r-help posts?

Marc R. Feldesman [EMAIL PROTECTED] writes: At the risk of being flamed for mortal stupidity, I'm trying to figure out what could have possibly changed (at my end almost certainly) that make Brian Ripley's posts (in particular) not show up on my mail server (at least not on the one where I

[R] error in gmake CrossCompileBuild

Dear all, I expect that I am on thin ice trying to build this on freeBSD 4.10, however, after tweaking Makefile-rcb v11Oct04 /* MAKE=gmake; MINGW_CROSS = mingw-cross4 */, I encounter: mingw32-gcc -isystem ~/RCrossBuild/cross- tools/mingw32/include -O2 -Wall -pedantic -I../include -I. -

Re: [R] error in gmake CrossCompileBuild

Steve Dutky [EMAIL PROTECTED] writes: tools/snip/include/float.h but does not appear to be defined elsewhere. Any advice greatly appreciated. Not this bit I fear, but: Please take it to the r-devel list; it is way beyond the level of r-help. -- O__ Peter Dalgaard

[R] test comparing spatial point patterns?

I'm working with animal relocation data. These are x-y coordinates of radio-collared coyotes over a period of time.Relocations are far enough apart in time to be considered independent. We want to test if the pattern of space-use has changed from one year to the next. Relocations roughly

RE: [R] avoiding loops

I didn't read the question carefully. You were asking all pairwise distances between the two sets, right? This isn't the most efficient way to do it, but it seems OK for problem of the size you stated: pairdist - function(m1, m2) { + idx1 - rep(1:nrow(m1), each=nrow(m2)) + idx2 -

Re: [R] Why am I not getting some r-help posts?

At 01:59 PM 1/28/2005, Peter Dalgaard wrote: Marc R. Feldesman [EMAIL PROTECTED] writes: At the risk of being flamed for mortal stupidity, I'm trying to figure out what could have possibly changed (at my end almost certainly) that make Brian Ripley's posts (in particular) not show up on my

[R] Database Connection Problem with RMySQL package

Folks, I failed to create a connection to the database under MySQL DBMS in the R system via RMySQL's method dbConnection(...). My setup is as follows: Microsoft Windows XP 5.1.2600 MySQL 4.1.9 R 2.0.1 DBI 0.1-8 RMySQL 0.5-5 Both of DBI and RMySQL packages were downloaded from bell lab: