Re: [R] boxplot issue

2010-10-15 Thread David Winsemius
A) you hijacked another thread. On Oct 15, 2010, at 9:50 AM, Jonas Josefsson wrote: Hi! I am trying to produce a graph which shows overlap in latitude for a number of species. I have a dataframe which looks as follows species1,species2,species3,species4. minlat

Re: [R] p-values from coxph?

2010-10-15 Thread David Winsemius
On Oct 15, 2010, at 9:21 AM, Öhagen Patrik wrote: Dear List, I each iteration of a simulation study, I would like to save the p- value generated by coxph. I fail to see how to adress the p-value. Do I have to calculate it myself from the Wald Test statistic? No. And the most important

Re: [R] interaction contrasts

2010-10-15 Thread David Winsemius
Should you need to do it again, you may want to look at the relevel function. I suppose that would meet the definition of some versions of on the fly but once I have a model, rerunning with a different factor leveling is generally pretty painless. -- David. On Oct 15, 2010, at 9:09 AM,

[R] [R-pkgs] solaR: version 0.20

2010-10-16 Thread David Winsemius
From: Oscar Perpiñan Lamigueiro oscar.perpi...@upm.es Date: October 8, 2010 4:25:34 AM EDT To: r-packa...@r-project.org Subject: [R-pkgs] solaR: version 0.20 Hello, I'd like to announce the availability of the version 0.20 of the solaR package. It provides a set of calculation methods of solar

Re: [R] delete data row

2010-10-16 Thread David Winsemius
IRD; There is a danger in applying logical tests of equality to floating point numbers. It may be safer to use all.equal or zapsmall in the construction of your tests. all.equal( (2^(0.5))^2 , 2) [1] TRUE (2^(0.5))^2 == 2 [1] FALSE -- David. On Oct 16, 2010, at 1:30 PM, Joshua Wiley

Re: [R] reduce the size of points in plot???

2010-10-16 Thread David Winsemius
On Oct 16, 2010, at 12:12 PM, Hongwei Dong wrote: Hi, R users, Can anyone tell me how I can change the size of points in my plot? For example: x - c(1,3,6,9,12) y - c(1.5,2,7,8,15) plot(x,y,pch=20) How do I reduce the size of those points? plot(x,y,pch=20, cex=0.2) plot(x,y,pch=20,

Re: [R] delete data row

2010-10-16 Thread David Winsemius
at 12:12 PM, David Winsemius dwinsem...@comcast.net wrote: IRD; There is a danger in applying logical tests of equality to floating point numbers. It may be safer to use all.equal or zapsmall in the construction of your tests. all.equal( (2^(0.5))^2 , 2) [1] TRUE (2^(0.5))^2 == 2 [1] FALSE

Re: [R] delete data row

2010-10-17 Thread David Winsemius
On Oct 17, 2010, at 3:56 PM, William Dunlap wrote: I had been thinking of: x - c(1, (2^(0.5))^2 , 3, 5, (2^(0.5))^2 , 3, 1) y - 2 x[-which(zapsmall(x-y) == 0)] [1] 1 3 5 3 1 Using which() to convert logicals into integer subscripts is almost always unnecessary and often wrong. At one

Re: [R] read.spss warning message (Unrecognized record type 7, subtype 18 encountered in system file)

2010-10-18 Thread David Winsemius
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] How to read only ten rows from a SAS dataset (read.ssd)?

2010-10-19 Thread David Winsemius
/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https

Re: [R] readLines: how to make a data.frame?

2010-10-19 Thread David Winsemius
-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] Error: object 'short' not found

2010-10-19 Thread David Winsemius
and the the ^^^ will not be aligned with the word short. You should fix the settings of you mail client in both these areas. Viki [[alternative HTML version deleted]] David Winsemius, MD West Hartford, CT

Re: [R] Chron object in time series plot

2010-10-19 Thread David Winsemius
://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing

Re: [R] Chron object in time series plot

2010-10-19 Thread David Winsemius
On Oct 19, 2010, at 12:19 PM, Manta wrote: David Winsemius wrote: You seen to be under the mistaken impression that the internal representation of DateTime classes of 08:00 would be 8. Since the internal representation of time is in seconds, the even number hours would be at integer

Re: [R] How to read only ten rows from a SAS dataset (read.ssd)?

2010-10-19 Thread David Winsemius
internal variable named something like _N_ that is a line number. Perhaps ... , if=_N_ = 10, OR: -- using SAS to output a smaller file. OR: --- adding a SAS line to the output returned within the read.ssd function. -- David. -J 2010/10/19 David Winsemius dwinsem...@comcast.net

Re: [R] Chron object in time series plot

2010-10-19 Thread David Winsemius
mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r

Re: [R] Problems with a specific calculate.

2010-10-19 Thread David Winsemius
://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list

Re: [R] scatter.smooth() fitted by loess

2010-10-19 Thread David Winsemius
-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] Gini Coefficient

2010-10-19 Thread David Winsemius
coefficient and add them together? Once again i am sorry if this is completely the wrong place to ask such a question. Peter David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] How to select not continous rows?

2010-10-20 Thread David Winsemius
On Oct 20, 2010, at 9:25 AM, skan wrote: Hello How can I select several not continuous rows ? If I wanted to select rows 1 to 7 I'll write mydata[,1:7] But what if I need to select rows 1 to 5 and 10 to 15? mydata[ , c(1:5, 10:15) ] -- David Winsemius, MD West Hartford, CT

Re: [R] number format, writing 1e-5 instead of 0.00001

2010-10-20 Thread David Winsemius
= 3) 0.1 [1] 0.1 0.001 [1] 1e-07 0.01 [1] 0.01 Other options for output: ?format ?sprintf -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] CI using ci.numeric

2010-10-20 Thread David Winsemius
-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help

Re: [R] (no subject)

2010-10-20 Thread David Winsemius
(mustart) : Value -1717986918 out of range (0, 1) Look more carefully at your 28th case. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] How to select not continous rows?

2010-10-20 Thread David Winsemius
On Oct 20, 2010, at 10:01 AM, David Winsemius wrote: On Oct 20, 2010, at 9:25 AM, skan wrote: Hello How can I select several not continuous rows ? If I wanted to select rows 1 to 7 I'll write mydata[,1:7] But what if I need to select rows 1 to 5 and 10 to 15? mydata[ , c(1:5, 10:15

Re: [R] remove black square from factor plot

2010-10-20 Thread David Winsemius
is produced by each species, there is a small black square in center position of the box. Is it possible to remove this. Have tried to include bxp(par(medpch=NA)) without success. Thanks! Jonas Josefsson David Winsemius, MD West Hartford, CT __ R

Re: [R] create a list fails

2010-10-20 Thread David Winsemius
])), total.patientnums.trt1=dat2[ ,2], total.patientnums.trt23=dat2[ ,2], num.countstrt1=dat2[ ,5], num.countstrt23=dat2[ ,6] ) === Dr. Jim Maas David Winsemius, MD West Hartford, CT __ R-help

Re: [R] Changing sign on absolute numbers 0 problems

2010-10-20 Thread David Winsemius
function: nonneg - function(x) { 0 + !sign(x) } nonneg(0) [1] 1 Any help is appreciated Thank you sadz ps please email me if you need more info -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] Changing sign on absolute numbers 0 problems

2010-10-20 Thread David Winsemius
On Oct 20, 2010, at 3:19 PM, David Winsemius wrote: On Oct 20, 2010, at 11:47 AM, Sadz A wrote: Hi, I am trying to do some calculations turning DMS data to decimal degrees using the formula (D+(M/60)+(S/3600)), some of the D's involve -ve numbers, the easiest way to do the calculation

Re: [R] question about masked object

2010-10-20 Thread David Winsemius
: namespace:Hmisc base::units function (x) UseMethod(units) environment: namespace:base Thank you and kind regards, gigi [[alternative HTML version deleted]] -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list

Re: [R] loading workspace- getting annoying

2010-10-20 Thread David Winsemius
also skipped loading my utility functions from .Rprofile as well as skipping the workspace .Rdata file. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] rounding up (always)

2010-10-20 Thread David Winsemius
, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented

Re: [R] All other variables in upper scope arg for stepAIC

2010-10-21 Thread David Winsemius
21.276 3.2365 + y11 0.149400 21.126 5.0956 + y31 0.025849 21.250 5.2122 Start: AIC=1 get(paste(y, i, sep = )) ~ 1 Df Sum of SqRSSAIC none 19.027 1.0023 + y11 0.227457 18.799 2.7618 + y21 0.023117 19.004 2.9780 -- David Winsemius, MD

Re: [R] Efficient nested loops

2010-10-21 Thread David Winsemius
the max of the last slice of the variable. Thanks in advance! -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

Re: [R] Different time between date() and Sys.date()

2010-10-21 Thread David Winsemius
, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self

Re: [R] 3D-scatterplots - high quality rendering?

2010-10-21 Thread David Winsemius
simply create a collection of images that I animate elsewhere? require(rgl) ?movie3d -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] Efficient nested loops

2010-10-21 Thread David Winsemius
/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https

Re: [R] Changing sign on absolute numbers 0 problems

2010-10-21 Thread David Winsemius
will be appreciated. Thank you sadz From: David Winsemius dwinsem...@comcast.net To: David Winsemius dwinsem...@comcast.net Cc: Sadz A sadz_a1...@yahoo.co.uk; r-help@r-project.org Sent: Wed, 20 October, 2010 20:45:39 Subject: Re: [R] Changing sign on absolute numbers 0 problems On Oct 20, 2010, at 3

Re: [R] exact pattern match in grep for column headers

2010-10-21 Thread David Winsemius
-- View this message in context: http://r.789695.n4.nabble.com/grep-tp3004422p3005885.html Sent from the R help mailing list archive at Nabble.com. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] data.frame query

2010-10-21 Thread David Winsemius
/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https

Re: [R] Different time between date() and Sys.date()

2010-10-21 Thread David Winsemius
setting and changing it in my system does the trick: Sys.time() [1] 2010-10-21 19:48:42 CEST -- David. Do you know where can I find a list of the timezone ? Message du 21/10/10 15:07 De : David Winsemius A : omerle Copie à : r-help@r-project.org Objet : Re: [R] Different time between

Re: [R] How to access values in s4 method

2010-10-21 Thread David Winsemius
3 dim(ovrl) [1] 3 3 -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal

Re: [R] printing a variable during a loop

2010-10-21 Thread David Winsemius
(0,0.025,0.25,0.5,0.75,0.975,1)) print(i) } I'd like to have the value of i printed for each loop (step). As I could see the values of i are shown on screen only after all the work is done. Thanks in advance for any suggestion. Best regards, Antonio -- David Winsemius, MD West Hartford, CT

Re: [R] previous business day

2010-10-21 Thread David Winsemius
On Oct 21, 2010, at 9:38 PM, Li, Jing Yi wrote: How to get the previous business day in R? I saw some post about using functions in timeSeries package before but can not find it anymore. ??holiday Perhaps you should look in the tis package. Thanks! David Winsemius, MD West Hartford

Re: [R] previous business day

2010-10-21 Thread David Winsemius
}} __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT

Re: [R] R step-by-step execution

2010-10-22 Thread David Winsemius
On Oct 22, 2010, alais wrote: |Moreover can you please tell me how I can concatenate variables and strings so to print some debugging messages / The value of x is + x + and the value of y is +y. Commas, not + cat(The value of x is , x , and the value of y is , y) -- David Winsemius

Re: [R] superscript characters in title with '+'

2010-10-22 Thread David Winsemius
^2+) 'onsets'), xlab = 'sec') David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented

Re: [R] issue with Matrix package

2010-10-22 Thread David Winsemius
de la classe derivedDefaultMethod any idea on how I could solve this issue ? Update, first. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] lm looking for weights outside of the user-defined function

2010-10-22 Thread David Winsemius
of formula. Why is it looking outside the function for the object that has just been defined inside the function? David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] lm looking for weights outside of the user-defined function

2010-10-22 Thread David Winsemius
the rules, I just follow them. I agree that one might guess that to be the search order, but it is not what is documented. -- David Winsemius, MD West Hartford, CT Dimitri On Fri, Oct 22, 2010 at 9:15 AM, David Winsemius dwinsem...@comcast.net wrote: On Oct 22, 2010, at 9:01 AM, Dimitri

Re: [R] superscript characters in title with '+'

2010-10-22 Thread David Winsemius
-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

Re: [R] Confusing error statement

2010-10-22 Thread David Winsemius
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org

Re: [R] lm looking for weights outside of the user-defined function

2010-10-22 Thread David Winsemius
On Oct 22, 2010, at 12:17 PM, William Dunlap wrote: -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of David Winsemius Sent: Friday, October 22, 2010 6:25 AM To: Dimitri Liakhovitski Cc: r-help Subject: Re: [R] lm looking

Re: [R] wait for graph to finish plotting

2010-10-22 Thread David Winsemius
On Oct 22, 2010, at 12:50 PM, n...@aleblanc.cotse.net wrote: Tal Galili tal.gal...@gmail.com writes: I suspect that using dev.copy2eps Is not going to help you here. Please try again using: pdf(...) # Check: ?pdf for(i in something) { plot(things) } dev.off() But give pdf() the

Re: [R] read.table input array

2010-10-22 Thread David Winsemius
On Oct 22, 2010, at 8:40 PM, Balpo wrote: Hello Jim How can I ensure this reading it from the file? The thing is, I cannot use the textConnection(), because the real file has millions of rows. The advice being offered is not that you use textConnection. That is just for illustration.

Re: [R] contour on a simplex

2010-10-23 Thread David Winsemius
If Dennis' interpretation of your request is correct then Harrell rms/ Hmisc packages provide a mechanism for restricting plotting of contours to the region within a perimeter. He has a perimeter function and the bplot function takes a perim argument. The reason I did not offer it earlier

Re: [R] help

2010-10-23 Thread David Winsemius
On Oct 23, 2010, at 4:13 AM, Liviu Andronic wrote: Hello On Sat, Oct 23, 2010 at 10:43 AM, Mastaki Kambale jkmast...@hotmail.com wrote: Dears R I am a self taught novice user of R. I begin to understand its philosophy and some basics like objects, vectors, arrays, lists etc...I still

Re: [R] convert wind direction from degrees to basic compass directions

2010-10-23 Thread David Winsemius
On Oct 23, 2010, at 1:40 PM, will phillips wrote: Hello, I have a data set that includes a predictor variable wind direction. This variable is in degrees with 0=North. I've put this predictor into a linear model and its coefficient is far from significant. I was thinking about

Re: [R] Long model formulae

2010-10-23 Thread David Winsemius
On Oct 23, 2010, at 6:56 PM, Gabor Grothendieck wrote: On Sat, Oct 23, 2010 at 9:51 PM, James Hirschorn james.hirsch...@hotmail.com wrote: What is a good way to enter a very long model formula. For example: y ~ Input.2 + Input.3 + ... + Input.1000 (assuming the corresponding dataframe has

Re: [R] If Statement Help

2010-10-23 Thread David Winsemius
On Oct 23, 2010, at 7:44 PM, Jason Kwok wrote: Thanks Jorge. It works. Is there a way to keep the actual price in the price column instead of TRUE/FALSE but filtering on when price100? Huh? When I use subset I get what you ask for: subset(x, Price 100) Price 2010-10-12 101

Re: [R] Contour Plot on a non Rectangular Grid

2010-10-24 Thread David Winsemius
On Oct 24, 2010, at 4:30 AM, Lorenzo Isella wrote: Dear All, I would like to plot a scalar (e.g. a temperature) on a non- rectangular domain (or even better: I would simply like to be able to draw a contour plot on an arbitrary 2D domain). I wonder if there is any tool to achieve that

Re: [R] Contour Plot on a non Rectangular Grid

2010-10-24 Thread David Winsemius
On Oct 24, 2010, at 6:12 AM, Lorenzo Isella wrote: As to the domain of the function, at least in case (1), that should arise from the collected data points in (x,y) if the sampling is dense enough. And that is precisely what you get from the perimeter function. The earlier Design

Re: [R] Contour Plot on a non Rectangular Grid

2010-10-24 Thread David Winsemius
On Oct 24, 2010, at 9:30 AM, Lorenzo Isella wrote: Hi, And thanks for helping. I am anyway a bit puzzled, since case (1) is not only a matter of interpolation. Probably the point I did not make clear (my fault) is that case (1) in my original email does not refer to an irregular grid

Re: [R] Contour Plot on a non Rectangular Grid

2010-10-25 Thread David Winsemius
On Oct 25, 2010, at 3:41 AM, Lorenzo Isella wrote: On 10/25/2010 01:32 AM, David Winsemius wrote: You were advised to look at rms. Why have you dismissed this suggestion? Using your data setup below and packaging into a dataframe. require(rms) ddf - datadist(xysf - as.data.frame(xys

Re: [R] How to save R file into specific type

2010-10-25 Thread David Winsemius
On Oct 25, 2010, at 1:07 AM, zhiji19 wrote: Hello everyone Can you please teach me how to save my homework as .R file? ?savehistory # at least on a Mac (If you wanted the whole console session with output, then it would be select-all, copy, paste into a text editor, save as a text

Re: [R] Artifacts in filled.contour+pdf

2010-10-25 Thread David Winsemius
On Oct 25, 2010, at 6:50 AM, Mario Valle wrote: Dear all, I'm using R 2.12.0 on Windows 7 (32bits) I created a filled contour from the attached data using the following code: load('bug.RData') pdf('bug.pdf', width=14, height=7) filled.contour(o4$x, o4$y, o4$z, color=rainbow, xlim=c(4,18),

Re: [R] divide column in a dataframe based on a character

2010-10-25 Thread David Winsemius
On Oct 25, 2010, at 8:56 PM, Daisy Englert Duursma wrote: Hello, If I have a dataframe: example(data.frame) zz- c (aa_bb ,bb_cc ,cc_dd,dd_ee,ee_ff,ff_gg,gg_hh,ii_jj,jj_kk,kk_ll) ddd - cbind(dd, group = zz) and I want to divide the column named group by the _, how would I do this?

Re: [R] (no subject)

2010-10-26 Thread David Winsemius
On Oct 25, 2010, at 11:17 PM, Arsalan Fathi wrote: hello. how are you? my name is arsalan. I'm from iran. i want to write the program that done random walk(one variable and two variable). please help me. You should learn to use a search engine that is specific to R. Here's an example:

Re: [R] Setting constraints in the glm package

2010-10-26 Thread David Winsemius
On Oct 26, 2010, at 7:50 AM, Viechtbauer Wolfgang (STAT) wrote: The constraint b1=b2 in a model such as b0 + b1 x1 + b2 x2 + b3 x3 implies that b0 + b1 (x1 + x2) + b3 x3, so just add x1 and x2 (call this x12) and fit the model b0 + b1 x12 + b3 x3 and you have imposed the constraint that

Re: [R] Forcing results from lm into datframe

2010-10-26 Thread David Winsemius
On Oct 26, 2010, at 8:08 AM, Small Sandy (NHS Greater Glasgow Clyde) wrote: Hi I need some help getting results from multiple linear models into a dataframe. Let me explain the problem. I have a dataframe with ejection fraction results measured over a number of quartiles and grouped

Re: [R] Zoom in in a plot

2010-10-26 Thread David Winsemius
On Oct 26, 2010, at 8:30 AM, Alaios wrote: in a simple plot. When i do plot is it possible to zoom in or out or this is not possible at all? Zoom? Do you mean restrict the region plotted to specific ranges? xlim and ylim arguments provide that facility. -- David. Best Regards Alex

Re: [R] Setting constraints in the glm package

2010-10-26 Thread David Winsemius
On Oct 26, 2010, at 9:27 AM, David Smith wrote: Many thanks for the help. You could express your thanks by including context the next time you present a follow-up (as requested in the Posting Guide). Only a minority of readers view this list on Nabble, so we don't see the web

Re: [R] Forcing results from lm into datframe

2010-10-26 Thread David Winsemius
thought maybe I sould see a cbind() call in there, and I suppose this section as.call(c(expression(data.frame), x ... may have that effect -- David. Why doesn't fitsdf - as.data.frame(t(fits)) work? Sandy Small From: David Winsemius [dwinsem

Re: [R] superscript characters in title with '+'

2010-10-26 Thread David Winsemius
On Oct 26, 2010, at 11:15 AM, Gavin Simpson wrote: On Fri, 2010-10-22 at 15:39 +0200, Claudia Beleites wrote: On 10/22/2010 03:15 PM, DrCJones wrote: snip / Being a chemist, it seemed natural to me to put the i after the concentration brackets into a subscript - though you didn't say

Re: [R] Cox Proportional Models and Haplotypes

2010-10-26 Thread David Winsemius
On Oct 26, 2010, at 8:57 AM, sr500 wrote: Hello, I was wondering if anyone knew of a function that fits haplotype data into a cox proportional hazard model. I have computed my Haplotype frequencies using the haplo.stats package. I have also been using the haplo.glm function but this is

Re: [R] ggplot - unwanted sorted X values

2010-10-27 Thread David Winsemius
as is typically done at time of factor creation. You probably need to work with factor levels. (All guesswork in the absence of a reproducible example.) -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] coxph linear.predictors

2010-10-27 Thread David Winsemius
in this area will get expert review and correction. Your comments are appreciated. Stephen Bond -Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: Wednesday, October 27, 2010 1:15 PM To: Bond, Stephen Cc: r-help@r-project.org Subject: Re: [R] coxph

Re: [R] How do I read multiple rows of different lengths?

2010-10-27 Thread David Winsemius
://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] (no subject)

2010-10-27 Thread David Winsemius
On Oct 27, 2010, at 2:20 PM, Czerminski, Ryszard wrote: I would appreciate any advice about how to add rectangular grid to the plot generated by levelplot() similarly as grid() function adds grid to the regular plot() ?panel.grid Best regards, Ryszard -- David Winsemius, MD West

Re: [R] How do I read multiple rows of different lengths?

2010-10-27 Thread David Winsemius
and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

Re: [R] How do I read multiple rows of different lengths?

2010-10-27 Thread David Winsemius
preprocess with readLines if you want an exact match of width or maybe take a guess with ... , colClasses = rep(numeric, 250) to create a sufficiently wide set of columns to hold everything. -- David Winsemius, MD West Hartford, CT __ R

Re: [R] Fwd: as.list

2010-10-27 Thread David Winsemius
-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT

Re: [R] how to add rectangular grid to the plot generated by levelplot()?

2010-10-27 Thread David Winsemius
plot() ?panel.grid Best regards, Ryszard -- David Winsemius, MD West Hartford, CT David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] bugs when using R-2.12.0 Please HELP!

2010-10-28 Thread David Winsemius
not a member). A word which has less implication of software deficiencies would be problems. This is really about BioConductor and it has its own mailing list. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https

Re: [R] runtime on ising model

2010-10-28 Thread David Winsemius
and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

Re: [R] runtime on ising model

2010-10-28 Thread David Winsemius
On Oct 28, 2010, at 12:20 PM, David Winsemius wrote: On Oct 28, 2010, at 11:52 AM, Michael D wrote: Mike, I'm not sure what you mean about removing foo but I think the method is sound in diagnosing a program issue and the results speak for themselves. I did invert my if statement

Re: [R] wait for graph to finish plotting

2010-10-28 Thread David Winsemius
: Could you please paste the exact code you are using? (the one with the pdf and dev.off, outside the loop ) Contact On Sat, Oct 23, 2010 at 4:53 PM, n...@aleblanc.cotse.net wrote: David Winsemius dwinsem...@comcast.net writes: On Oct 22, 2010, at 12:50 PM, n

Re: [R] Clustering

2010-10-28 Thread David Winsemius
package. Can anyone suggest a more appropriate package to use for such a large dataset? David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

Re: [R] Clustering

2010-10-29 Thread David Winsemius
of the same event, hence the time conditon. Hope this clarifies the problem. It clarifies it to the extent that it show how much more you will need to further clarify. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list

Re: [R] strftime vs strptime ??

2010-10-29 Thread David Winsemius
. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained

Re: [R] Clustering

2010-10-29 Thread David Winsemius
. Those would be character values and, if you want to do calculations, would obviously need to be coerced to numeric. Thanks, Doug -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo

Re: [R] Clustering

2010-10-29 Thread David Winsemius
On Oct 29, 2010, at 12:08 PM, David Winsemius wrote: On Oct 29, 2010, at 11:37 AM, dpender wrote: Apologies for being vague, The structure of the output is as follows: Still no code? $ cluster1 : Named num [1:131] 3.05 2.71 3.26 2.91 2.88 3.11 3.21 -1 2.97 3.39 ... ..- attr

Re: [R] Differenciate numbers from reference for rows

2010-10-30 Thread David Winsemius
-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented

Re: [R] Clustering

2010-10-30 Thread David Winsemius
On Oct 30, 2010, at 7:49 AM, dpender wrote: David Winsemius wrote: On Oct 29, 2010, at 12:08 PM, David Winsemius wrote: On Oct 29, 2010, at 11:37 AM, dpender wrote: Apologies for being vague, The structure of the output is as follows: Still no code? I am using the Clusters

Re: [R] Differenciate numbers from reference for rows

2010-10-30 Thread David Winsemius
I am missing) these results are concatenated in a string, and then evaluated, a step which I do get. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] solving equation consisting of factorials

2010-10-30 Thread David Winsemius
David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained

Re: [R] Using names function

2010-10-30 Thread David Winsemius
month: max(diff(msales)) I get the numeric value displayed when I execute the above function, but not the month. cumsum(msales) displays the results under each month, but not max. How can I get both the month and value displayed? diff(msales)[which.max(diff(msales))] Mar 450 -- David

Re: [R] for loop

2010-10-30 Thread David Winsemius
example applying split to the builtin airquality dataframe. The plyr package also provides functions on dataframes. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Confidence interval for response variable in mixed effects models

2010-10-30 Thread David Winsemius
http://finzi.psych.upenn.edu/R/Rhelp02/archive/76742.html http://finzi.psych.upenn.edu/R/Rhelp02/archive/81237.html http://finzi.psych.upenn.edu/R/Rhelp02/archive/82567.html -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing

Re: [R] for loop

2010-10-30 Thread David Winsemius
the results creating a list of plots. -- DAvid. Thanks,m -Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: Saturday, October 30, 2010 8:24 PM To: Matevž Pavlič Cc: r-help@r-project.org Subject: Re: [R] for loop On Oct 30, 2010, at 2:07 PM, Matevž Pavlič

Re: [R] two group cox model

2010-10-30 Thread David Winsemius
. It seems to create a solid and dashed line in the legend when I use it. ?par # 1=solid (default), 2=dashed, Many thanks in advance. Cheba David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https

  1   2   3   4   5   6   7   8   9   10   >