Re: [R] kruskal-wallis, stratified

2010-04-13 Thread Heinz Tuechler
See the thread stratified Wilcoxon available? at http://tolstoy.newcastle.edu.au/R/help/05/08/11143.html Heinz At 11:21 13.04.2010, Kay Cichini wrote: hello everyone, can anybody tell me if there is a kruskal-wallis, or another non-parametric test, that can deal with multiple samples that

Re: [R] kruskal-wallis, stratifiedhttp://n4.nabble.com/forum/NewNode.jtp?tpl=replynode=1838232

2010-04-13 Thread Heinz Tuechler
Sorry for not being precise enough. Here http://tolstoy.newcastle.edu.au/R/help/05/08/11177.html you should find the attachment http://tolstoy.newcastle.edu.au/R/help/att-11177/KW.strat.2005.R I used it, and it seems to work. In some cases some elements of weight may become Inf. Heinz At

Re: [R] Can we get rid of bar charts with error bars?

2009-12-03 Thread Heinz Tuechler
Frank, the example on http://biostat.mc.vanderbilt.edu/DynamitePlots is nice, and I agree with you. Just one minor question: would it be possible to mention as An article with nice dot plots a paper, which is freely available? Heinz At 14:56 03.12.2009, Frank E Harrell Jr wrote: Bar charts

Re: [R] Building static HTML help pages in R 2.10.x on Windows

2009-12-27 Thread Heinz Tuechler
Instead of an answer, may I add question c) can someone state that it is impossible to generate static HTML help pages under Windows? At 21:40 22.12.2009, Steve Rowley wrote: I upgraded to R2.10.1pat and discovered, along with everybody else, that static HTML pages are no longer the default.

Re: [R] Applying function to parts of a matrix based on a factor

2010-01-13 Thread Heinz Tuechler
If your matrix were a data.frame, it could work like this: df - data.frame(age=1:100, sex=rep(1:2, 50)) with(df, by(age, sex, mean)) without the lapply, sapply etc. family. h At 18:16 13.01.2010, Doran, Harold wrote: with(yourdataframe, tapply(age,sex,mean)) -Original Message- From:

Re: [R] merging issue.........

2010-01-13 Thread Heinz Tuechler
Did you consider to look at the help page for merge? h At 22:01 13.01.2010, karena wrote: hi, I have a question about merging two files. For example, I have two files, the first file is like the following: id trait1 110.2 211.1 39.7 610.2 78.9 10 9.7 11 10.2 The

Re: [R] BMDP and SAS (was R in clinical trials)

2010-02-22 Thread Heinz Tuechler
Once I suggested to BMDP to introduce a module-statement that would direct the syntax to the specified module (1L, 2L, ...), so that all syntax could reside in one job, but they did not like that idea. Heinz At 14:55 19.02.2010, Terry Therneau wrote: I used both BMDP and SAS in my earlier

Re: [R] Output mean/median survival time from survfit

2008-02-05 Thread Heinz Tuechler
Maybe this thread is of use for you. How to access results of survival analysis Xiaochun Li (06 May 2006) http://tolstoy.newcastle.edu.au/R/help/06/05/26713.html Heinz At 21:28 04.02.2008, Xing Yuan wrote: Hi all, Does anybody know how to output the mean/median survival time from survfit? Thank

Re: [R] Problem with cut

2008-02-22 Thread Heinz Tuechler
At 15:22 22.02.2008, Henrique Dallazuanna wrote: Try this: grep(330, levels(cc), value=T) Could you please explain in a little more detail, how this answers the original question? I would have expected 330 to fall into (313,330] category. Can you please advice what do I do wrong? Thank you

Re: [R] correlation between categorical data

2009-06-21 Thread Heinz Tuechler
At 07:40 21.06.2009, J Dougherty wrote: [...] There are other ways of regarding the FET. Since it is precisely what it says - an exact test - you can argue that you should avoid carrying over any conclusions drawn about the small population the test was applied to and employing them in a

[R] survSplit with data.frame containing a Surv object

2009-07-13 Thread Heinz Tuechler
Dear All, since years I am struggling with Surv objects in data.frames. The following seems to have to do with it. See below the modified example from the help page of survSplit. The original works, as expected. If, however, a Surv object is added to the data.frame, each record gets doubled.

Re: [R] survSplit with data.frame containing a Surv object

2009-07-13 Thread Heinz Tuechler
At 20:18 13.07.2009, Charles C. Berry wrote: On Mon, 13 Jul 2009, Heinz Tuechler wrote: Dear All, since years I am struggling with Surv objects in data.frames. The following seems to have to do with it. See below the modified example from the help page of survSplit. The original works

Re: [R] Two envelopes problem

2008-08-26 Thread Heinz Tuechler
Mark My experience was similarly frustrating. Maybe formulating the problem a bit differently could help to clarify it. State it like this: Someone chooses an amount of money x. He puts 2x/3 of it in one envelope and x/3 in an other. There is no assumption about the distribution of x. If you

Re: [R] Mode Vs Class

2008-04-09 Thread Heinz Tuechler
Congratulation Bill for this very clear and useful explanation. Heinz At 14:58 08.04.2008, [EMAIL PROTECTED] wrote: 'mode' is a mutually exclusive classification of objects according to their basic structure. The 'atomic' modes are numeric, complex, charcter and logical. Recursive objects have

Re: [R] Trend test for survival data

2008-04-22 Thread Heinz Tuechler
Dear Markus! Since I did not see an answer yet, my suggestion is to use coxph with the groups variable numerically coded as the only independent variable. Heinz At 13:39 21.04.2008, Markus Kreuz wrote: Hello, is there a R package that provides a log rank trend test for survival data in =3

Re: [R] survival curves for time dependent covariates (was consultation)

2009-05-13 Thread Heinz Tuechler
At 14:50 12.05.2009, Terry Therneau wrote: *I´m writing to ask you how can I do Survivals Curves using Time-dependent *covariates? Which packages I need to Install?* This is a very difficult problem statistically. That is, there are not many good ideas for what SHOULD be done. Hence,

[R] Where to find a changelog for the survival package

2009-05-20 Thread Heinz Tuechler
Dear All, since some days I try to use the versions 2.35-4 of the survival package instead of versions 2.31, I had installed until now. Several changes in print.survfit, plot.survfit and seemingly in the structure of ratetabels effect some of my syntax files. Is there somewhere a

Re: [R] Where to find a changelog for the survival package

2009-05-20 Thread Heinz Tuechler
Thank you Richie. I had seen this before, but my impression is that it's not up to date. I gave a wrong version number in my previous post. I changed from 2.34-1 to 2.35-4. For example, the plot.survfit function lost it's legend parameters, but I don't see this in the changelog. Thanks again,

Re: [R] Changelog for the survival package

2009-05-21 Thread Heinz Tuechler
Dear Terry, first of all, thank you for your immense work. At the moment, I don't have a small reproducible example for the ratetable difficulty I have. I will work on it. Maybe the error message I get is of some information to you. Error in match.ratetable(m[, rate], ratetable) : Data

Re: [R] survfit, summary, and survmean (was Changelog for survival package)

2009-05-22 Thread Heinz Tuechler
Dear Terry, sorry that I did not see this change, and thank you for it. It is very useful. Heinz At 15:28 22.05.2009, Terry Therneau wrote: Further I appreciate your new function survmean(). At the moment it seems to be intended as internal, and not documented in the help. The

Re: [R] Different results in calculating SD of 2 numbers

2008-01-16 Thread Heinz Tuechler
At 11:32 16.01.2008, Jim Lemon wrote: (Ted Harding) wrote: On 16-Jan-08 08:45:04, Martin Maechler wrote: RM == Ron Michael [EMAIL PROTECTED] on Wed, 16 Jan 2008 00:14:56 -0800 (PST) writes: RM Hi all, RM Can anyone tell me why I am getting different results in calculating SD

[R] How to keep attributes when dropping factor levels?

2009-03-19 Thread Heinz Tuechler
Dear All, to drop unused factor levels two ways are outlined in R-help. In both cases a label attribute is lost. The same happens, when using car:::recode. Is there a simple way to avoid losing attributes? Thanks, Heinz ## example ff - factor(substring(statistics, 1:10, 1:10),

[R] factor, as.factor and levels

2009-04-08 Thread Heinz Tuechler
Dear All, to my surprise as.factor does not accept a levels argument. Maybe I did not read the documentation well enough. See the example below. I wanted to use ch1 as factor in the newdata argument of survfit, so I assumed that I could write as.factor(ch1, levels=ch1), since the order

Re: [R] factor, as.factor and levels

2009-04-08 Thread Heinz Tuechler
does not accept levels as an argument. use the first form that you have factor(ch1, levels=ch1) On Wed, Apr 8, 2009 at 7:36 AM, Heinz Tuechler tuech...@gmx.at wrote: Dear All, to my surprise as.factor does not accept a levels argument. Maybe I did not read the documentation well enough. See

Re: [R] factor, as.factor and levels

2009-04-08 Thread Heinz Tuechler
if it is a factor first. Here is the definition: ) as.factor function (x) if (is.factor(x)) x else factor(x) environment: namespace:base You can always list out what the function does to get a better understanding of how it works. On Wed, Apr 8, 2009 at 8:16 AM, Heinz Tuechler tuech...@gmx.at

Re: [R] (senza oggetto)

2009-04-16 Thread Heinz Tuechler
At 11:10 16.04.2009, giuseppef...@libero.it wrote: Dear all, I have a database x,y,value imported in R with read.table: dati- read.table(dati.dat) value is a categorical data (land use) and i want to plot in the same colour the same land use. It is possible with R. Thanks a lot

[R] Umlaut read from csv-file

2008-11-06 Thread Heinz Tuechler
Dear All! Reading character strings containing an umlaut from a csv-file I find a (to me) surprising behaviour in R 2.8.0, that I did not notice in R 2.7.2. A comparison by == results in FALSE, while grep does find the aggreement. See the example below. The crucial line is x==div 1-2

Re: [R] Umlaut read from csv-file

2008-11-06 Thread Heinz Tuechler
there is a problem in the use of the CHARSXP cache: if I save the session then x0 == x becomes true when I reload it, even though the encodings remain different. I've found the immediate cause and will change this in R-patched shortly. On Thu, 6 Nov 2008, Heinz Tuechler wrote: Dear All! Reading character

[R] Encoding() and strsplit()

2008-11-06 Thread Heinz Tuechler
Dear All, Encoding() goes beyond my understanding. See the example. I would expect from reading the help for Encoding() that strsplit preserves the encoding for each resulting element, but for simple letters it gets lost. Also it seems that an Encoding() cannot be declared for simple

Re: [R] Encoding() and strsplit()

2008-11-07 Thread Heinz Tuechler
), including levels of factors to latin1? Heinz Tüchler ASCII characters are not marked as Latin-1 nor UTF-8. On Fri, 7 Nov 2008, Heinz Tuechler wrote: Dear All, Encoding() goes beyond my understanding. See the example. I would expect from reading the help for Encoding() that strsplit preserves

Re: [R] Umlaut read from csv-file

2008-11-07 Thread Heinz Tuechler
At 13:34 07.11.2008, Peter Dalgaard wrote: Heinz Tuechler wrote: Dear Prof.Ripley! Thank you very much for your attention. In the given example Encoding(), or the encoding parameter of read.csv solve the problem. I hope your patch will solve also the problem, when I read a spss file

Re: [R] Mismatch in logical result?

2008-11-07 Thread Heinz Tuechler
Maybe this? http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f At 11:23 07.11.2008, Shubha Vishwanath Karanth wrote: Content-Type: text/plain Content-Disposition: inline Content-length: 569 Hi R, I have certain checkings, which gives FALSE,

Re: [R] Umlaut read from csv-file

2008-11-07 Thread Heinz Tuechler
At 16:52 07.11.2008, Prof Brian Ripley wrote: On Fri, 7 Nov 2008, Peter Dalgaard wrote: Heinz Tuechler wrote: Dear Prof.Ripley! Thank you very much for your attention. In the given example Encoding(), or the encoding parameter of read.csv solve the problem. I hope your patch will solve also

Re: [R] Umlaut read from csv-file

2008-11-08 Thread Heinz Tuechler
immagine, would be some function, which marks every character string in the whole structure of a data.frame, including all attributes, as latin1. On Sat, 8 Nov 2008, Heinz Tuechler wrote: At 16:52 07.11.2008, Prof Brian Ripley wrote: On Fri, 7 Nov 2008, Peter Dalgaard wrote: Heinz Tuechler

[R] attr.all.equal() and all.equal(attributes(), attributes())

2008-11-09 Thread Heinz Tuechler
Dear All! If I try to compare the attributes of two objects, I find a surprising behaviour of attr.all.equal(). With identical attributes I receive the answert NULL. If the attributes differ, the answer is as expecxted and differences are shown. all.equal(attributes(), attributes()) instead

Re: [R] Compare objects

2008-11-09 Thread Heinz Tuechler
At 13:26 09.11.2008, Leon Yee wrote: Hi, friends Is there any functions for object comparing? For example, I have two list objects, and I want to know whether they are the same. Since the the components of list are not necessary atomic, this kind of comparison should be recursive. Does this

Re: [R] Umlaut read from csv-file

2008-11-09 Thread Heinz Tuechler
At 06:25 09.11.2008, Prof Brian Ripley wrote: On Sat, 8 Nov 2008, Heinz Tuechler wrote: At 08:01 08.11.2008, Prof Brian Ripley wrote: We have no idea what you understood (you didn't tell us), but the help says encoding: character vector. The encoding(s) to be assumed when 'file

Re: [R] attr.all.equal() and all.equal(attributes(), attributes())

2008-11-09 Thread Heinz Tuechler
At 14:24 09.11.2008, Peter Dalgaard wrote: Heinz Tuechler wrote: Dear All! If I try to compare the attributes of two objects, I find a surprising behaviour of attr.all.equal(). With identical attributes I receive the answert NULL. If the attributes differ, the answer is as expecxted

Re: [R] survival::survfit,plot.survfit

2009-02-27 Thread Heinz Tuechler
At 15:28 26.02.2009, Terry Therneau wrote: plot(survfit(fit)) should plot the survival-function for x=0 or equivalently beta'=0. This curve is independent of any covariates. This is not correct. It plots the curve for a hypothetical subject with x= mean of each covariate. Does this

[R] How to write a Surv object to a csv-file?

2008-12-19 Thread Heinz Tuechler
Dear All, trying to write a data.frame, containing Surv objects to a csv-file I get Error in dimnames(X) - list(dn[[1L]], unlist(collabs, use.names = FALSE)) : length of 'dimnames' [2] not equal to array extent. See example below. May be, I overlooked something, but I expected that also

Re: [R] How to write a Surv object to a csv-file?

2008-12-19 Thread Heinz Tuechler
I tried the compare package, I had no success comparing data.frames containing Surv objects. Thanks again Heinz At 22:31 19.12.2008, David Winsemius wrote: On Dec 19, 2008, at 2:04 PM, Heinz Tuechler wrote: Dear All, trying to write a data.frame, containing Surv objects to a csv-file I

Re: [R] How to write a Surv object to a csv-file?

2008-12-20 Thread Heinz Tuechler
Dear Charles, yes, your solution does what I need. Maybe, it offers also a way to use the compare package with Surv objects. Thank you, Heinz At 23:30 19.12.2008, Charles C. Berry wrote: On Fri, 19 Dec 2008, Heinz Tuechler wrote: Dear David! Thank you for your response. I like csv files

Re: [R] How many attributes are there of a variable?

2009-09-07 Thread Heinz Tuechler
Peng, based on a suggestion, Frank made years ago (18.7.2006), I use one attribute that contains all further attributes, I want to assign to variables. It's necessary to create your own class and subsetting method, so that this attribute does not get lost. Together with some functions I use

Re: [R] What determines the unit of POSIXct differences?

2009-09-11 Thread Heinz Tuechler
differences in secs [1] 86400 3542400 On Fri, Sep 11, 2009 at 7:50 AM, Heinz Tuechler tuech...@gmx.at wrote: Dear All, what determines if a difference between POSIXct objects gets expressed in days or seconds? In the following example, it's sometimes seconds, sometimes days. as.POSIXct('2009

[R] How to assess object names within a function in lapply or l_ply?

2009-09-28 Thread Heinz Tuechler
Dear All, to produce output of several columns of a data frame, I tried to use lapply and also l_ply. In both cases, I would like to print a header line containing also the name of the respective column in the data frame. For example, I would like the following lapply(data.frame(a=1:3,

Re: [R] How to assess object names within a function in lapply or l_ply?

2009-09-28 Thread Heinz Tuechler
, Henrique Dallazuanna wrote: You can use names insteed: DF - data.frame(a=1:3, b=2:4) lapply(names(DF), function(x){ print(x) DF[x] }) On Mon, Sep 28, 2009 at 8:22 AM, Heinz Tuechler tuech...@gmx.at wrote: Dear

Re: [R] How to assess object names within a function in lapply or l_ply?

2009-09-28 Thread Heinz Tuechler
Dallazuanna wrote: Heinz, Try this: lapply(DF, function(x)names(DF)[as.numeric(gsub([^0-9], , deparse(substitute(x]) On Mon, Sep 28, 2009 at 8:43 AM, Heinz Tuechler tuech...@gmx.at wrote: Thank you, Henrique, my example was simplified. In a more complexe function I want to use the objects

Re: [R] How to assess object names within a function in lapply or l_ply?

2009-09-28 Thread Heinz Tuechler
Hadley, many thanks for your answer and for the enormous work you put into plyr, a really powerful package. For now, I will solve my problem with a variable label attribute, I usually attach to columns in data frames. I asked the list, because I thought, I am overlooking something trivial,

Re: [R] How to assess object names within a function in lapply or l_ply?

2009-09-28 Thread Heinz Tuechler
At 18:17 28.09.2009, hadley wickham wrote: many thanks for your answer and for the enormous work you put into plyr, a really powerful package. For now, I will solve my problem with a variable label attribute, I usually attach to columns in data frames. I asked the list, because I thought, I

Re: [R] Different result of multiple regression in R and SPSS

2011-07-20 Thread Heinz Tuechler
At 19.07.2011 18:50 -0700, Spencer Graves wrote: On 7/19/2011 4:04 PM, Bert Gunter wrote: On Tue, Jul 19, 2011 at 3:45 PM, David Winsemiusdwinsem...@comcast.net wrote: On Jul 19, 2011, at 6:29 PM, J. wrote: Thanks for the answer. # However, I am still curious about

Re: [R] value.labels

2011-08-11 Thread Heinz Tuechler
At 11.08.2011 21:50 +0300, Zeki Çatav wrote: Prş, 2011-08-11 tarihinde 19:27 +0200 saatinde, Uwe Ligges yazdı: On 11.08.2011 19:22, David Winsemius wrote: On Aug 11, 2011, at 11:42 AM, Uwe Ligges wrote: On 11.08.2011 16:10, zcatav wrote: Hello R people, I have a data.frame.

Re: [R] value.labels

2011-08-11 Thread Heinz Tuechler
At 12.08.2011 09:11 +1200, Rolf Turner wrote: On 12/08/11 09:59, Heinz Tuechler wrote: At 11.08.2011 21:50 +0300, Zeki Çatav wrote: Prş, 2011-08-11 tarihinde 19:27 +0200 saatinde, Uwe Ligges yazdı: On 11.08.2011 19:22, David Winsemius wrote: On Aug 11, 2011, at 11:42 AM, Uwe Ligges

Re: [R] value.labels

2011-08-11 Thread Heinz Tuechler
At 12.08.2011 11:05 +1200, Rolf Turner wrote: On 12/08/11 11:34, Heinz Tuechler wrote: At 12.08.2011 09:11 +1200, Rolf Turner wrote: On 12/08/11 09:59, Heinz Tuechler wrote: At 11.08.2011 21:50 +0300, Zeki Çatav wrote: Prş, 2011-08-11 tarihinde 19:27 +0200 saatinde, Uwe Ligges yazdÄ

Re: [R] value.labels

2011-08-11 Thread Heinz Tuechler
At 12.08.2011 01:53 +0100, Heinz Tuechler wrote: At 12.08.2011 11:05 +1200, Rolf Turner wrote: On 12/08/11 11:34, Heinz Tuechler wrote: At 12.08.2011 09:11 +1200, Rolf Turner wrote: On 12/08/11 09:59, Heinz Tuechler wrote: At 11.08.2011 21:50 +0300, Zeki Çatav wrote: Prş, 2011-08-11

[R] retain class after merge

2011-08-19 Thread Heinz Tuechler
Dear All, is there a simple way to retain the class attribute of a column, if merging two data.frames? When merging the example data.frames form help(merge) I am unable to keep the class attribute as set before merging (see below). Two columns are assigned new classes before merge (myclass1,

Re: [R] Feature request: rating/review system for R packages

2011-03-20 Thread Heinz Tuechler
It's unclear to me, why the rating/review system should relate to entire packages. Would it not be more informative, if single specific functions would be rated and reviewed? I would like to see if + is rated better than -, or if more difficulties are reported for * than for /. I could then

Re: [R] Why do we have to turn factors into characters for various functions?

2010-12-12 Thread Heinz Tuechler
At 12.12.2010 00:48 +0200, Tal Galili wrote: Hello dear R-help mailing list, My question is *not* about how factors are implemented in R (which is, if I understand correctly, that factors keeps numbers and assign levels to them). My question *is* about why so many functions that work on factors

Re: [R] Why do we have to turn factors into characters for various functions?

2010-12-13 Thread Heinz Tuechler
Hello Petr, don't want to convince you. If you like the following: x - factor(1:4, labels=c(one, two, three, four)) y - factor(3:5, labels=c(three, four, five)) data.frame(character=c(as.character(x), as.character(y)), numeric=c(x, y)) character numeric 1 one 1 2 two

[R] survexp - unable to reproduce example

2010-12-20 Thread Heinz Tuechler
Dear All, when I try to reproduce an example of survexp, taken from the help page of survdiff, I receive the error message Error in floor(temp) : Non-numeric argument to mathematical function . It seems to come from match.ratetable. I think, it has to do with character variables in a

[R] survexp - example produces error

2010-12-31 Thread Heinz Tuechler
Dear All, reposting, because I did not find a solution, maybe someone could check the example below. It's taken from the help page of survdiff. Executing it, gives the error Error in floor(temp) : Non-numeric argument to mathematical function best regards, Heinz library(survival) ##

Re: [R] survexp - example produces error

2010-12-31 Thread Heinz Tuechler
wrote: On Dec 31, 2010, at 10:21 , Heinz Tuechler wrote: Dear All, reposting, because I did not find a solution, maybe someone could check the example below. It's taken from the help page of survdiff. Executing it, gives the error Error in floor(temp) : Non-numeric argument to mathematical

Re: [R] survexp - example produces error

2010-12-31 Thread Heinz Tuechler
[61] 669 29 619 595 89 16 544 20 514 95 481 444 427 79 333 [76] 396 109 369 206 185 339 264 164 179 130 108 30 10 Browse[1] Q At 31.12.2010 13:46 +0100, peter dalgaard wrote: On Dec 31, 2010, at 10:21 , Heinz Tuechler wrote: Dear All, reposting, because

Re: [R] survexp - example produces error

2010-12-31 Thread Heinz Tuechler
=old.stringsAsFactors) it seems to work. Heinz At 31.12.2010 15:53 +0100, Heinz Tuechler wrote: Dear Peter, Dear All, a further attempt led me to an answer. If I set options(stringsAsFactors=TRUE), which I usually have set to FALSE, no error occurs. I am, however not happy with this solution. Heinz

Re: [R] CSV value not being read as it appears

2011-01-14 Thread Heinz Tuechler
At 14.01.2011 07:09 -0800, Peter Ehlers wrote: On 2011-01-14 02:09, bgr...@dyson.brisnet.org.au wrote: Brian, Thanks. My response to David follows. I should add that this problem has never occurred previously as far as I know (I have now checked the previous report I was sent): This problem

[R] big speed difference in source btw. R 2.15.2 and R 3.0.2 ?

2013-10-29 Thread Heinz Tuechler
Dear All, is it known that source works much faster in R 2.15.2 than in R 3.0.2 ? In the example below I observe e.g. for a data.frame with 10^7 rows the following timings: R version 2.15.2 Patched (2012-11-29 r61184) length: 1e+07 user system elapsed 62.040.22 62.26 R version

Re: [R] big speed difference in source btw. R 2.15.2 and R 3.0.2 ?

2013-10-30 Thread Heinz Tuechler
All was run on the identical machine in independent sessions. I did not restart Windows. I also tried 32bit R 3.0.2 and it seemed slightly faster than 64bit. Using Process Explorer v15.23 (http://technet.microsoft.com/de-de/sysinternals/bb896653) my impression was that R 3.0.2 manages memory

Re: [R] big speed difference in source btw. R 2.15.2 and R 3.0.2 ?

2013-10-30 Thread Heinz Tuechler
Best thanks for confirming my impression. I use dump for storing large data.frames with a number of attributes for each variable. save/load is much faster, but I am unsure, if such files will be readable by R versions years later. What format/functions would you suggest for data

Re: [R] big speed difference in source btw. R 2.15.2 and R 3.0.2 ?

2013-10-31 Thread Heinz Tuechler
memory. Thank you for the hint. It appears to me that source() in R-devel performs at about the same speed as in R 2.15.2. Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com -Original Message- From: Heinz Tuechler [mailto:tuech...@gmx.at] Sent: Wednesday, October 30, 2013 1:43 PM

Re: [R] Comparing Cox model with Competing Risk model

2013-03-16 Thread Heinz Tuechler
Dear Terry, as soon as the vignette is ready, I would be very happy, to know about it. Will you send a note to r-help, or will it be announced in some other way? best regards, Heinz On 08.03.2013 15:12, Terry Therneau wrote: -- begin included message -- I have a competing risk data where

Re: [R] correlation between categorical data

2015-01-28 Thread Heinz Tuechler
comment inline David Winsemius wrote on 24.01.2015 21:08: On Jan 23, 2015, at 5:54 PM, JohnDee wrote: Heinz Tuechler wrote At 07:40 21.06.2009, J Dougherty wrote: [...] There are other ways of regarding the FET. Since it is precisely what it says - an exact test - you can argue that you

Re: [R] Probable Error in fmsb package

2015-12-31 Thread Heinz Tuechler
Anindya Sankar Dey wrote/hat geschrieben on/am 30.12.2015 07:35: Hi All, The fmsb package has a function called Variance Inflation Factor and it states the definition of the function as follows:- "To evaluate multicolinearity of multiple regression model, calculating the variance inflation

Re: [R] p values from GLM

2016-04-02 Thread Heinz Tuechler
Bert Gunter wrote on 01.04.2016 23:46: ... of course, whether one **should** get them is questionable... http://www.nature.com/news/statisticians-issue-warning-over-misuse-of-p-values-1.19503#/ref-link-1 This paper repeats the common place statement that a small p-value does not necessarily

Re: [R] post_processor in rmarkdown not working

2017-09-07 Thread Heinz Tuechler
Are you sure that you want to read in the output_file in text <- readLines(output_file, warn = FALSE)? best regards, Heinz Thierry Onkelinx wrote/hat geschrieben on/am 06.09.2017 11:41: Dear all, I'm trying to write a post_processor() for a custom rmarkdown format. The goal of the

Re: [R] boot.stepAIC fails with computed formula

2017-08-23 Thread Heinz Tuechler
It seems that if you build the formula as a character string, and postpone the "as.formula" into the lm call, it works. instead of frm1 <- as.formula(paste(trg,"~1")) use frm1a <- paste(trg,"~1") and then strt <- lm(as.formula(frm1a),dat) regards, Heinz Stephen O'hagan wrote/hat geschrieben

Re: [R] SE for all fixed factor effect in GLMM

2018-12-29 Thread Heinz Tuechler
maybe qvcalc https://cran.r-project.org/web/packages/qvcalc/index.html is useful for you. Marc Girondot via R-help wrote/hat geschrieben on/am 30.12.2018 05:31: Dear members, Let do a example of simple GLMM with x and G as fixed factors and R as random factor: (note that question is the same

Re: [R] Value Labels: SPSS Dataset to R

2020-02-07 Thread Heinz Tuechler
Maybe it helps searching at https://rseek.org/ for "SPSS to R transition value labels". In particular https://cran.r-project.org/web/packages/expss/vignettes/labels-support.html seems useful, as well as https://www.r-bloggers.com/migrating-from-spss-to-r-rstats/ best regards, Heinz Jim Lemon

Re: [R] My dream ...

2020-05-12 Thread Heinz Tuechler
Abby Spurdle wrote/hat geschrieben on/am 12.05.2020 10:38: In my opinion the advantage of computers is not Artificial Intelligence, but rather Artificial Patience (most AI that I have seen is really doing a bunch of what I would consider to be boring, really fast so people don't have to). Leave

Re: [R] Some seemingly odd behavior of survfit in library survival

2020-03-17 Thread Heinz Tuechler
Jeff Newmiller wrote/hat geschrieben on/am 17.03.2020 05:39: The coxph function appears to rely on finding the name of the data argument in the environment in which the formula was created. The lm function does not have this problem. Oh, and df is the name of the F distribution density

Re: [R] How to convert European short dates to ISO format?

2020-06-10 Thread Heinz Tuechler
maybe isoDates <- as.Date(oriDates, format = "%d/%m/%y") Heinz Luigi Marongiu wrote/hat geschrieben on/am 10.06.2020 10:20: Hello, I have been trying to convert European short dates formatted as dd/mm/yy into the ISO 8601 but the function as.Dates interprets them as American ones (mm/dd/yy),

Re: [R] svglite with multiple files

2020-12-09 Thread Heinz Tuechler
hat people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, Dec 9, 2020 at 3:33 AM Heinz Tuechler wrote: Dear All, while svg() (package grDevices) can produce several files, svglite() (packag

[R] svglite with multiple files

2020-12-09 Thread Heinz Tuechler
Dear All, while svg() (package grDevices) can produce several files, svglite() (package svglite) is limited to one file/page only (as documented in the respective help page). Is there a simple solution to make svglite() work like svg() to produce several files? Of course one could call svglite()

Re: [R] Inappropriate color name

2020-11-19 Thread Heinz Tuechler
inline - David Wright wrote on 19.11.2020 12:39: Appropriation of Indian Red as 'Chestnut' (or other alternative) will be viewed by some as 'making appropriate' the label for a colour, and no doubt by other groups as cultural theft by excising reference to its origin. Seems the best option is

Re: [R] Defining partial list of variables

2021-01-05 Thread Heinz Tuechler
What about the Cs()-function in Hmisc? library(Hmisc) Cs(a,b,c) [1] "a" "b" "c" Steven Yen wrote/hat geschrieben on/am 05.01.2021 13:29: Thanks Eric. Yes, "unlist" makes a difference. Below, I am doing not regression but summary to keep the example simple. > set.seed(123) >

Re: [R] Defining partial list of variables

2021-01-05 Thread Heinz Tuechler
see below Steven Yen wrote/hat geschrieben on/am 05.01.2021 08:14: I constantly define variable lists from a data frame (e.g., to define a regression equation). Line 3 below does just that. Placing each variable name in quotation marks is too much work especially for a long list so I do that

Re: [R] Translation of the charter

2021-11-02 Thread Heinz Tuechler
Greg, Greg Minshall wrote/hat geschrieben on/am 02.11.2021 08:57: Heinz, x <- c("a","b","c") lettersnum <- 1:length(letters[]) names(lettersnum) <- letters[] lettersnum[x] lettersnum[x] a b c 1 2 3 i'm not sure if the following is obviously better, but one might do b <- match(a, a)

Re: [R] Translation of the charter

2021-11-02 Thread Heinz Tuechler
Alice wrote/hat geschrieben on/am 31.10.2021 07:33: Dear members, How to translate the charter to the underline inter? I tried this: x <- c("a","b","c") as.numeric(x) [1] NA NA NA Warning message: NAs introduced by coercion It didn't work. Sorry for my newbie questions. B.R.