[R] font size

2006-03-26 Thread linda.s
Dear Friends, How to make the font size I changed in the R console (Mac Machine) permanent? Now I had to set the font size every time I opened the console. Thanks, L. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] load huge image

2006-03-26 Thread Gottfried Gruber
hello, i have run around 65000 regressions and stored them in a list. then i stored the session with save.image on my hard disk. the file is almost 1GB. when i now want to load the image it took tons of time. even after 12h of loading it was not done, although the saving was done fairly fast.

Re: [R] Regexp subexpression

2006-03-26 Thread Dieter Menne
Gabor Grothendiek wrote: read.table(textConnection(sub(pat, '\\1 \\2', patid)), as.is = TRUE) V1 V2 1 ALAN 334 2 AzD 44 3 NA Looks like you are suffering from a creativity push. I like that one, with no success I had tried the variant sub(pat, c(\\1, \\2') patid)) before.

Re: [R] R gets slow

2006-03-26 Thread pau carre
Hello, I am still having the same problem, even removing the objects (with rm()). I checked the memory allocated by R and it is constant (about 4.5%), so I have no idea of what is happening... Moreover in few minutes (15) the models are generated 5 times slower... restarting the server is the only

Re: [R] MenuRead() Question

2006-03-26 Thread Keith Chamberlain
Dear List-mates, I think the difficulty I'm having is localized to the MenuType() call made from within MenuRead(). I'm not used to seeing text operations such as ^ or [], so am having trouble understanding what's going on. I'm interested in understanding the regular expression:

Re: [R] R-help Digest, Vol 37, Issue 26

2006-03-26 Thread isaac . martin
Mi nueva dirección de correo es: [EMAIL PROTECTED] New e-mail address: [EMAIL PROTECTED] __ 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

[R] how to remove the framebox out the 3D scatter plot?

2006-03-26 Thread Jinsong Zhao
Hellow, I plot a 3D scatter plot using cloud function in lattice package. Although, I have read the help of cloud, and panel.cloud, I did not find a way to remove the framebox out the 3D scatter plot. Any help? Thanks in advance! Best wishes, Jinsong

[R] write read object from mysql

2006-03-26 Thread Gottfried Gruber
hi, i want to store an lm-object in mysql and use the following command: dbSendQuery(con,paste(INSERT INTO tablename (id,name,result,result_txt,rsquared) VALUES (,x,,',z2,',',serialize(z1,NULL,ascii=F),',',serialize(z1,NULL,ascii=T),',,summary(z1)$r.squared,),sep=)) z1 is the lm-object z2 is

[R] new to the list - problems with non-solid lines in eps export

2006-03-26 Thread Gregor Volberg
Dear all, I am new to this list and, unfortunately, could not provide help to anyone as yet. I hope I can do so in the future, though! Until then, I am grateful for helpful hints from you more experienced users. For use in an upcoming publication, I generated an eps figure from my x11 window

Re: [R] Newbie clustering/classification question

2006-03-26 Thread Sean Davis
Mark A. Miller wrote: My laboratory is measuring the abundance of various proteins in the blood from either healthy individuals or from individuals with various diseases. I would like to determine which proteins, if any, have significantly different abundances between the healthy and

Re: [R] new to the list - problems with non-solid lines in eps export

2006-03-26 Thread Uwe Ligges
Gregor Volberg wrote: Dear all, I am new to this list and, unfortunately, could not provide help to anyone as yet. I hope I can do so in the future, though! Until then, I am grateful for helpful hints from you more experienced users. For use in an upcoming publication, I generated an eps

Re: [R] VR (looking for MASS) on OS X

2006-03-26 Thread Uwe Ligges
John Christie wrote: Hey, when did I lose MASS? Anyway, I found the VR package but she's not compiling on OS X. Any tips? Yes, I have the developer tools and Fortran. Could someone send me a compiled package? Package bundle VR (and hence package MASS) is included in release

[R] range and bwplot

2006-03-26 Thread vincent david
Hi, is there an equivalent to the 'range' option of the boxplot function to be found in the bwplot function of the trellis package? regards Vincent [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list

Re: [R] how to remove the framebox out the 3D scatter plot?

2006-03-26 Thread Uwe Ligges
Jinsong Zhao wrote: Hellow, I plot a 3D scatter plot using cloud function in lattice package. Although, I have read the help of cloud, and panel.cloud, I did not find a way to remove the framebox out the 3D scatter plot. Any help? You can set te argument par.box =

Re: [R] range and bwplot

2006-03-26 Thread Uwe Ligges
vincent david wrote: Hi, is there an equivalent to the 'range' option of the boxplot function to be found in the bwplot function of the trellis package? Please read the help files. ?bwplot points you to ?panel.bwplot which explains that you should use argument coef. Uwe Ligges regards

Re: [R] range and bwplot

2006-03-26 Thread Sundar Dorai-Raj
vincent david wrote: Hi, is there an equivalent to the 'range' option of the boxplot function to be found in the bwplot function of the trellis package? regards Vincent Hi, Vincent, From ?panel.bwplot, you can supply your own function for stats. For your case, you would need:

Re: [R] object size vs. file size

2006-03-26 Thread Duncan Murdoch
On 3/25/2006 10:16 PM, Gabor Grothendieck wrote: You can place functions in lists or environments and pass the environment to the function and have it look there first. That way you can have different versions of a function with the same name. 1. Here is an example using lists: A - list(f

[R] merge two vectors of the same length into one?

2006-03-26 Thread Zhiqiang Ye
Dear all, I have a question about merge two character vectors of the same length into one. For example: X= c('A', 'B', 'C', 'A', 'D', ...) # 26 possible alphabet values with length 1000 Y=c('B', 'A', 'C', 'A', 'C', ) # 26 possible aphabet values with

Re: [R] merge two vectors of the same length into one?

2006-03-26 Thread Sarah Goslee
What about this: X - c(A, B, C, A, D, A) Y - c(B, A, C, A, C, B) XY - paste(X, Y, sep=) sort(table(XY), decreasing=TRUE) XY AB AA BA CC DC 2 1 1 1 1 On 3/26/06, Zhiqiang Ye [EMAIL PROTECTED] wrote: Dear all, I have a question about merge two character vectors of the same

Re: [R] merge two vectors of the same length into one?

2006-03-26 Thread Duncan Murdoch
On 3/26/2006 11:26 AM, Zhiqiang Ye wrote: Dear all, I have a question about merge two character vectors of the same length into one. For example: X= c('A', 'B', 'C', 'A', 'D', ...) # 26 possible alphabet values with length 1000 Y=c('B', 'A', 'C', 'A',

Re: [R] MenuRead() Question

2006-03-26 Thread Philippe Grosjean
Hi Keith, If you want to define a Tk menu using MenuRead() and a text file to define your menu, you have to start it with Tk.. That way, MenuRead() recognizes that it is a Tk menu. So, rewrite your menu definition file as: menu.txt - in RHOME $Tk.KSesnMain |$MenuTest ||Objects

Re: [R] merge two vectors of the same length into one?

2006-03-26 Thread Zhiqiang Ye
AB AA BA CC DC 2 1 1 1 1 Yeah, thank you two for your help. BTW, can I change the output like this way: AB 2 AA 1 BA 1 CC 1 DC 1 Thanks a lot! -- Zhiqiang Ye __ R-help@stat.math.ethz.ch mailing list

Re: [R] MenuRead() Question

2006-03-26 Thread Keith Chamberlain
Dear Gabor Philipe, Thank you for clarifying regular expressions, and how tk menus are demarcated in menus.txt. Regular expressions are beginning to make sense for me, and the menu is now being added properly to the tk window. regexpr(^[$]Tk[.].+/, menu)... If I understand better, the

[R] Problem with names() in a plot after ordering a data.frame. Syntax or stupidity?

2006-03-26 Thread John Kane
I am a complete newbie in R . Using R 2.2.0 Windows XP This started as a simple exercise to see if I could produce a simple Cleveland style dotchart with a line from the '2' axis to the dot rather than a continuous line as do dotchart and dotchart2. At least I could not find that

Re: [R] MenuRead() Question

2006-03-26 Thread Prof Brian Ripley
On Sun, 26 Mar 2006, Keith Chamberlain wrote: Dear Gabor Philipe, Thank you for clarifying regular expressions, and how tk menus are demarcated in menus.txt. Regular expressions are beginning to make sense for me, and the menu is now being added properly to the tk window.

Re: [R] Problem with names() in a plot after ordering a data.frame. Syntax or stupidity?

2006-03-26 Thread Gabor Grothendieck
Try this or suitable variation: with(provs[order(provs$Area),], { idx - seq(Pcode) plot(Area, idx, col = 2, yaxt = n) segments(0, idx, Area, idx, col = 2) axis(2, idx, as.character(Pcode)) }) On 3/26/06, John Kane [EMAIL PROTECTED] wrote: I am a complete newbie

[R] addition using binary

2006-03-26 Thread Fred J.
Dear R users I looked around for a package which can help me with the task of binary arithmetic, the closest I found is sfsmisc with digitsBase function, which may help a little in this task, but still wondering how to get this problem solved. example: given the addition operation of 2

Re: [R] merge two vectors of the same length into one?

2006-03-26 Thread Peter Dalgaard
Zhiqiang Ye [EMAIL PROTECTED] writes: AB AA BA CC DC 2 1 1 1 1 Yeah, thank you two for your help. BTW, can I change the output like this way: AB 2 AA 1 BA 1 CC 1 DC 1 Just feed it to cbind() -- O__ Peter Dalgaard Øster Farimagsgade 5, Entr.B

Re: [R] R gets slow

2006-03-26 Thread Liaw, Andy
I'm just taking wild guesses here, and I have no clue what could be the problem. - The timing you mentioned, I assume, is `wall clock' time? Have you checked whether the CPU time also increase, or just the wall clock time? - Perhaps you can use Rprof() to check which part of the computation

Re: [R] R gets slow

2006-03-26 Thread pau carre
Hello, I finally solved the problem. The solution was just remove the attach operations. When attach is called a copy of the data is generated and since I know it is not possible to delete it with rm as it is not associated to any object. Pau 2006/3/26, pau carre [EMAIL PROTECTED]: Hello, I am

Re: [R] Problem with names() in a plot after ordering a data.frame. Syntax or stupidity?

2006-03-26 Thread Deepayan Sarkar
On 3/26/06, John Kane [EMAIL PROTECTED] wrote: I am a complete newbie in R . Using R 2.2.0 Windows XP This started as a simple exercise to see if I could produce a simple Cleveland style dotchart with a line from the '2' axis to the dot rather than a continuous line as do dotchart

Re: [R] R gets slow

2006-03-26 Thread Peter Dalgaard
pau carre [EMAIL PROTECTED] writes: Hello, I finally solved the problem. The solution was just remove the attach operations. When attach is called a copy of the data is generated and since I know it is not possible to delete it with rm as it is not associated to any object. You can detach()

Re: [R] dot diagram

2006-03-26 Thread Liaw, Andy
hist() returns an object that contains information for you to construct what you want, I believe, so you should be able do something like: obj - hist(..., plot=FALSE) and contruct your dot diagram based on obj. See the `Value' section of ?hist. Andy From: Jinsong Zhao Hello, I hope to

[R] generating random events (dd-mm-yyyy)

2006-03-26 Thread Carsten Steinhoff
Hello, I have to generate a couple of random datasets to validate an actuarial software. Each set has to contain at least two columns: (1) Date of the loss events occurrence; (2) Severity of loss I'm generating the dates on a monthly basis by, for example: ts(matrix(rpois(36,10)),

Re: [R] generating random events (dd-mm-yyyy)

2006-03-26 Thread jim holtman
try this: data.frame(dates = format(seq(ISOdate(2003,1,1), by='month', length=36), + format='%d.%m.%Y'), date=rpois(36,10)) dates date 1 01.01.2003 11 2 01.02.20038 3 01.03.20037 4 01.04.20037 5 01.05.20035 6 01.06.2003 12 7 01.07.20039 8 01.08.2003

Re: [R] VR (looking for MASS) on OS X

2006-03-26 Thread John Christie
OK, that worked great. But, oddly enough, in the Mac OS X version of R the packages are not in the distribution, and VR does not come up as a binary in the GUI (only as source). Anyway, it didn't come up for me. I did try reinstalling. On Mar 26, 2006, at 9:20 AM, Uwe Ligges wrote:

Re: [R] post hoc comparison following with multcomp?

2006-03-26 Thread Spencer Graves
Your use of the multcomp package looks plausible to me (though I certainly did not check every detail). Do you have any specific concerns? [Also, have you reviewed vignette(Rmc)?] And have you considered the plotting, e.g.:

Re: [R] R package for computing state path using Viterbi algorithm

2006-03-26 Thread Spencer Graves
Have you tried www.bioconductor.org? If you are studying nucleatide sequences, I would think you might get better help from Bioconductor. If you've already tried that (and you haven't received replies I've missed), it looks like you've found an opportunity to develop new

Re: [R] lme convergence

2006-03-26 Thread Spencer Graves
Dear Dr. Chaudhary: Because your example is not self contained (i.e., your data.frame family is not something distributed with R), I can only guess. Have you tried things like the following: with(family, table(cengirth, familyid)) If each level of

Re: [R] NLME Covariates

2006-03-26 Thread Spencer Graves
What you think about the following: set.seed(1) DF0.3 - data.frame(X=c(a,a, b), y=rnorm(3)) lme(y~1, random=~1|X, data=DF0.3) Linear mixed-effects model fit by REML Data: DF0.3 Log-restricted-likelihood: -2.148692 Fixed: y ~ 1 (Intercept) -0.4261464 Random effects:

[R] Clustering question \ dist(datmat)

2006-03-26 Thread kumar zaman
Hello everybody. I am trying to cluster circular data (data points which are angles), thus i can not use the dist function in mclust to generate my distance matrix, I am using the function Dij = 0.5*( 1 - cos(theta_i - theta_j)). The thing is hclust will not accept this distance matrix, i

Re: [R] Clustering question \ dist(datmat)

2006-03-26 Thread Gabor Grothendieck
A distance matrix must be of class dist. Try hclust(as.dist(df)) On 3/26/06, kumar zaman [EMAIL PROTECTED] wrote: Hello everybody. I am trying to cluster circular data (data points which are angles), thus i can not use the dist function in mclust to generate my distance matrix, I am using

[R] Matrix elements

2006-03-26 Thread Akkineni,Vasundhara
Hello all, I have a matrix object, xx given as below: [,1] [,2] [1,] 0.02237883 0.601572660 [2,] -0.39880918 0.126498958 [3,] 0.20269214 0.567953402 [4,] -0.45750812 -0.031193600 [5,] -0.30666134 -0.084819484 [6,] -0.37718928 0.078675868 [7,] -0.25432685

Re: [R] Matrix elements

2006-03-26 Thread Liaw, Andy
Something like this should work: which(rowSums(xx 0) == ncol(xx)) Andy From: Akkineni,Vasundhara Hello all, I have a matrix object, xx given as below: [,1] [,2] [1,] 0.02237883 0.601572660 [2,] -0.39880918 0.126498958 [3,] 0.20269214 0.567953402 [4,]

Re: [R] Matrix elements

2006-03-26 Thread Christos Hatzis
Since both columns being positive implies that the row wise sum and row wise product are both positive, you can try the following: x - matrix(rnorm(20),ncol=2) x[apply(x, 1, FUN=function(x) ifelse(sum(x)0 prod(x)0,TRUE,FALSE)),] -Christos Hatzis -Original Message- From: [EMAIL

Re: [R] Matrix elements

2006-03-26 Thread Christos Hatzis
To follow up on my previous posting, if you just need the row numbers try this instead: row(x)[apply(x, 1, FUN=function(x) ifelse(sum(x)0 prod(x)0,TRUE,FALSE)),1] -Christos -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Akkineni,Vasundhara Sent:

Re: [R] Clustering question \ dist(datmat)

2006-03-26 Thread kumar zaman
Dear Gabor and all ; I know this will work; but i already have a distance matrix calculated using my distance measure Dij = 0.5 * ( 1 - cos(theta_i - theta_j)), if i do hclust(as.dist(df)) then i am taking distance another time for a matrix df which is supposed to be a distance matrix,

Re: [R] Matrix elements

2006-03-26 Thread Gabor Grothendieck
Try this: which(xx[,1] 0 xx[,2] 0) On 3/26/06, Akkineni,Vasundhara [EMAIL PROTECTED] wrote: Hello all, I have a matrix object, xx given as below: [,1] [,2] [1,] 0.02237883 0.601572660 [2,] -0.39880918 0.126498958 [3,] 0.20269214 0.567953402 [4,]

Re: [R] Clustering question \ dist(datmat)

2006-03-26 Thread Liaw, Andy
as.dist() does _not_ recompute the distances if given a matrix. It simply takes the lower triangular portion of the distance matrix given and attach some attributes about the original dimension. I don't think you need to object to that. Andy From: kumar zaman Dear Gabor and all ; I

Re: [R] Clustering question \ dist(datmat)

2006-03-26 Thread kumar zaman
Liaw and Gabor : Thank you a bunch, you both are right, i just doubled check my data it just put things the way hclust likes. Thank you again Ahmed Florida Liaw, Andy [EMAIL PROTECTED] wrote: as.dist() does _not_ recompute the distances if given a matrix. It simply takes the

Re: [R] generating random events (dd-mm-yyyy)

2006-03-26 Thread Carsten Steinhoff
Hi Jim, that's not exactly what I search. In your solution a sequence of dates is generated and a poisson-frequency separatly. But I have to create for example a poisson distributed sequence. E.g.: January 2003 = 3 Events; February 2003 = 4 Events 01.01.2003 02.01.2003 03.01.2003 01.02.2003

Re: [R] Matrix elements

2006-03-26 Thread Uwe Ligges
Gabor Grothendieck wrote: Try this: which(xx[,1] 0 xx[,2] 0) ... or more generally for an arbitrary number of columns: which(apply(xx 0, 1, all)) Uwe Ligges On 3/26/06, Akkineni,Vasundhara [EMAIL PROTECTED] wrote: Hello all, I have a matrix object, xx given as below:

Re: [R] data frame as X in linear model lm() ?

2006-03-26 Thread Prof Brian Ripley
I guess you are looking for print(summary(lm(y ~ ., data=x.d))) since '.' refers to all the columns of 'data' except perhaps the response. On Sat, 25 Mar 2006, ivo welch wrote: Dear R wizards: This must have an obvious solution, but I am stumped. I can run a linear regression giving

Re: [R] load huge image

2006-03-26 Thread Martin Maechler
Gottfried == Gottfried Gruber [EMAIL PROTECTED] on Sun, 26 Mar 2006 10:27:35 +0200 writes: Gottfried hello, i have run around 65000 regressions and Gottfried stored them in a list. then i stored the session Gottfried with save.image on my hard disk. the file is Gottfried