[R] Text wrapping in R

2010-10-24 Thread Johannes Huesing
I am about to give an introduction to R to some clinical data managers used to SAS. There is already a lot of material in printed form and on the web that paves the way. What I haven't found so far are text wrapping capabilities in setting tables in raw text as in SAS PROC REPORT. At the moment

Re: [R] combining plots (curve + Plot functions)

2010-11-01 Thread Johannes Huesing
mms...@comcast.net mms...@comcast.net [Mon, Nov 01, 2010 at 04:11:53AM CET]: Hello,   What I really want to do is to add a rejection region in the form of a long rectangle to a density plot I have drawn.  I am getting  2 plots.  How can I add rectangle to first plot?  see code below.

Re: [R] stacking effectively - without loops

2010-11-02 Thread Johannes Huesing
Dimitri Liakhovitski dimitri.liakhovit...@gmail.com [Tue, Nov 02, 2010 at 09:57:04PM CET]: Hello! I have 2 vectors: x-letters[1:5] y-1:3 Is there a way - without loops - to create a data frame such that we repeat the whole y within each level of x so that it looks like this: a 1

Re: [R] setting attributes

2010-11-05 Thread Johannes Huesing
David Winsemius dwinsem...@comcast.net [Fri, Nov 05, 2010 at 05:18:23AM CET]: [...] If you were trying for something else, then please explain in plain English rather than simply showing code that throws an error. Well, code never hurts. But it should be accompanied with a hint towards what

Re: [R] Sweave question

2010-11-13 Thread Johannes Huesing
Ralf B ralf.bie...@gmail.com [Sat, Nov 13, 2010 at 10:03:49PM CET]: It seems that Sweave is supposed to be used from Latex and R is called during the LaTeX compilation process whenever R chunks appear. This is not how it works. In the first page of

Re: [R] diverting output from nested loops

2010-11-24 Thread Johannes Huesing
Lara Poplarski larapoplar...@gmail.com [Wed, Nov 24, 2010 at 06:14:58PM CET]: for (K in ...){ ... create object ... topDirName - as.character(paste(topName, K, sep=)) topDirMake - dir.create(topDirName) You may want to insert setwd(topDirName) here. for (L in ...) {

Re: [R] Statistical mailing list

2010-08-15 Thread Johannes Huesing
Tal Galili tal.gal...@gmail.com [Thu, Jul 29, 2010 at 03:42:28PM CEST]: This is the new and exiting kid on the block: http://stats.stackexchange.com/ Sorry for being late to this thread, but while looking somewhat exciting, this does not seem to be a mailing list. -- Johannes Hüsing

Re: [R] Combining two variables in text

2010-09-28 Thread Johannes Huesing
trb1 thomasrbol...@yahoo.co.uk [Tue, Sep 28, 2010 at 12:12:01PM CEST]: [...] Also, how can I combine several variables in text i.e. suppose I have a - one b - two then what expression can I place in title(main=...) in terms of a and b to give an output title of one, two i.e. a then b

[R] Sweave and LaTeX beamer class

2010-09-30 Thread Johannes Huesing
I am failing to uncover Sweave chunks step by step using the LaTeX beamer class. The following minimal example: \documentclass{beamer} \usepackage{Sweave} \begin{document} \begin{frame}[fragile] In the year \uncover2-{25}\uncover3-{\Sexpr{5*5}} \uncover4-{ echo=TRUE, print=TRUE= 5*5*101 @ }

Re: [R] Sweave and LaTeX beamer class

2010-10-01 Thread Johannes Huesing
Charles C. Berry cbe...@tajo.ucsd.edu [Fri, Oct 01, 2010 at 01:24:43AM CEST]: [...] \begin{uncoverenv}4- echo=TRUE, print=TRUE= 5*5*101 @ %def \end{uncoverenv} Thanks! Seems I was not familiar enough with the beamer documentation. -- Johannes Hüsing There is something

Re: [R] Non-Parametric Adventures in R

2010-10-03 Thread Johannes Huesing
Jamesp james.jrp...@gmail.com [Sat, Oct 02, 2010 at 11:27:09PM CEST]: [...] 1) I was thinking I'd have to go through each nominal variable (i.e. table(X$race) ), but I think I have it figured

Re: [R] Mapping the coordinates!

2010-10-09 Thread Johannes Huesing
Mehdi Zarrei gagzar...@yahoo.com [Sun, Oct 10, 2010 at 06:11:23AM CEST]: Hello, I have a series of coordinates (latitudes and longitudes) each one/several associated to a code (from 1 to 28). I used function points (latitude,

Re: [R] Reading in and writing out one line at a time

2010-05-21 Thread Johannes Huesing
Jan van der Laan djvanderl...@gmail.com [Fri, May 21, 2010 at 11:22:35AM CEST]: Perhaps you mean something like sapply or apply? Neither am I certain what he means exactly, but I would use apply and cousins only if I don't care in which sequence the elements are processed. His problem

Re: [R] GPS finding and mark interesting POI within table with large GPS points collection

2010-05-21 Thread Johannes Huesing
Eugeniusz Kałuża eugeniusz.kal...@polsl.pl [Fri, May 21, 2010 at 04:47:21PM CEST]: Dear R users, # #There is table containing 1000 (lat, lon, time) GPS positions, wchich

Re: [R] Learning R

2009-11-30 Thread Johannes Huesing
Esmail esmail...@gmail.com [Mon, Nov 30, 2009 at 01:07:15PM CET]: ps: Just checked, 'R in a Nutshell': http://oreilly.com/catalog/9780596801717 release date dec 2009/jan 2010 You can't really judge a book by its cover, and in this case not even the cover is presented to us. But

Re: [R] Sampling with replacement

2010-06-16 Thread Johannes Huesing
Somnath Somnath somnath700...@gmail.com [Wed, Jun 16, 2010 at 07:27:32PM CEST]: Thanks for all those reply. Is there any general rule to determine how many samples I would get from a population of size n, I draw a sample of size m (m may be greater than n) if sample is drawn with replacement?

[R] Read code from character string

2010-06-17 Thread Johannes Huesing
Dear expRts, I have a character string, say a - print(9**2). How do I execute the contents of the string, parsed as R code? Do I have to open a connection and use cat(a), and parse it at the other end? Cheers Johannes -- Johannes Hüsing There is something fascinating about

Re: [R] Read code from character string

2010-06-18 Thread Johannes Huesing
Peter Langfelder peter.langfel...@gmail.com [Thu, Jun 17, 2010 at 09:38:45PM CEST]: eval(parse(text=print(9**2))) cheers, I overlooked the text option. On Thu, Jun 17, 2010 at 12:32 PM, Johannes Huesing johan...@huesing.name wrote: Dear expRts, I have a character string, say a - print

Re: [R] Read code from character string

2010-06-18 Thread Johannes Huesing
Greg Snow greg.s...@imail.org [Fri, Jun 18, 2010 at 04:57:03PM CEST]: You should also look at fortune(106) and think about possible other solutions to your overall objective. I am not installing fortune solely for this purpose but I understand that anything which smells like macro expansion is

Re: [R] Read code from character string

2010-06-19 Thread Johannes Huesing
Greg Snow greg.s...@imail.org [Sat, Jun 19, 2010 at 12:44:48AM CEST]: There are a lot of other reasons to install the fortunes package that just the one fortune, there is much wisdom, some wit, (and then there are mine) throughout the package. Sorry, but I am subscribing to r-help already.

Re: [R] Read code from character string

2010-06-20 Thread Johannes Huesing
Gabor Grothendieck ggrothendi...@gmail.com [Sat, Jun 19, 2010 at 11:19:35AM CEST]: On Sat, Jun 19, 2010 at 4:56 AM, Johannes Huesing johan...@huesing.name wrote: A rule would be a data structure containing the ID of the rule, the rule in human readable language, an expression evaluating

Re: [R] Replacing elements of a list over a certain threshold

2010-06-22 Thread Johannes Huesing
Jim Hargreaves ja...@ipec.co.uk [Mon, Jun 21, 2010 at 12:34:01PM CEST]: Dear List, I have a list of length ~1000 filled with numerics. I need to replace the elements of this list that are above a certain numerical threshold with the value of the threshold. e.g example=list(1, 2, 3, 4, 5,

Re: [R] Basic question - more efficient method than loop?

2010-06-28 Thread Johannes Huesing
GL pfl...@shands.ufl.edu [Mon, Jun 28, 2010 at 05:46:13PM CEST]: I'm guessing there's a more efficient way to do the following using the index features of R. Appreciate any thoughts 1st thought: ifelse() for (i in 1:nrow(dbs1)){ if(dbs1$Payor[i] %in% Payor.Group.Medicaid)

Re: [R] Basic question - more efficient method than loop?

2010-06-28 Thread Johannes Huesing
Johannes Huesing johan...@huesing.name [Mon, Jun 28, 2010 at 06:31:20PM CEST]: [...] eval(parse(paste(Payor.Group, lst[[1]], sep=.))), eval(parse(text=paste(Payor.Group, lst[[1]], sep=.))), -- Johannes Hüsing There is something fascinating about science

Re: [R] Novel (Maybe?) Visualizations

2010-02-10 Thread Johannes Huesing
D. Dashcle d.dasc...@hotmail.com [Wed, Feb 10, 2010 at 01:13:58AM CET]: I'm interested in using R's plotting capabilities to try to generate a graph showing the relationship/pairing frequency of words appearing in a block of unstructured text. I don't have a specific algorithm or

Re: [R] Reading from Google Docs

2009-07-09 Thread Johannes Huesing
Gabor Grothendieck ggrothendi...@gmail.com [Wed, Jul 08, 2009 at 10:31:30PM CEST]: [...]    [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] HI

2008-09-16 Thread Johannes Huesing
Amit Patel [EMAIL PROTECTED] [Tue, Sep 16, 2008 at 03:32:01PM CEST]: Does anyone know an easy way to convert all the zero values in a imported csv table into NA's Depends on the data structure you gave your imported table. In a single numeric vector (named, say, vec), the syntax is

[R] Accessing methods and extending S4 classes of existing packages

2008-06-29 Thread Johannes Huesing
Dear all, I am trying to understand how to access S4 methods after loading a package, using the online documentation of getMethod and friends. This is what I have been trying: library(coin) findMethods(ApproxNullDistribution) list() Warning message: In findMethods(ApproxNullDistribution) :

Re: [R] Accessing methods and extending S4 classes of existing packages

2008-06-29 Thread Johannes Huesing
Johannes Huesing [EMAIL PROTECTED] [Sun, Jun 29, 2008 at 08:13:27AM CEST]: Dear all, I am trying to understand how to access S4 methods after loading a package, using the online documentation of getMethod and friends. I am just realizing that the lack of online documentation on S4 methods

Re: [R] Accessing methods and extending S4 classes of existing packages

2008-06-29 Thread Johannes Huesing
Johannes Huesing [EMAIL PROTECTED] [Sun, Jun 29, 2008 at 08:13:27AM CEST]: Dear all, I am trying to understand how to access S4 methods after loading a package, using the online documentation of getMethod and friends. This is what I have been trying: library(coin) findMethods

Re: [R] Sunset in Dortmund

2008-07-02 Thread Johannes Huesing
Marc Schwartz [EMAIL PROTECTED] [Wed, Jul 02, 2008 at 06:22:40PM CEST]: For those attending useR! 2008 next month, something to look forward to: http://www.gnome-look.org/content/show.php/show.php?content=84455 THIS is a sunset over Dortmund: http://www.23hq.com/mBusolt/photo/2230816 Deep

Re: [R] Sunset in Dortmund

2008-07-03 Thread Johannes Huesing
Marc Schwartz [EMAIL PROTECTED] [Thu, Jul 03, 2008 at 05:57:01PM CEST]: Lots of train fans would pay extra for a hotel with such a view At the peril of taking this thread further astray, I agree, having been one in my youth growing up in and around New York City in the 60's and 70's.

Re: [R] Interpreting messages when building packages

2008-07-07 Thread Johannes Huesing
The earlier messages suggest that you have automatically produced help files (from package.skeleton) and haven't edited them. I'd suggest moving them to a temporary directory until you're ready to edit them properly. ok, that get's me a little bit closer to a minimal working example as

Re: [R] Interpreting messages when building packages

2008-07-07 Thread Johannes Huesing
Duncan Murdoch [EMAIL PROTECTED] [Mon, Jul 07, 2008 at 09:24:07PM CEST]: You didn't say how you were doing the install, but the directory name looks like one that was used by R CMD check (.../dynalc.Rcheck/dynalc). No, that was an attempt to install during a call to R CMD check. --

Re: [R] Interpreting messages when building packages

2008-07-07 Thread Johannes Huesing
Duncan Murdoch [EMAIL PROTECTED] [Mon, Jul 07, 2008 at 09:57:50PM CEST]: Does a plain install succeed? No, the error messages look familiar: * Installing to library '/usr/local/lib/R/site-library' * Installing *source* package 'dynalc' ... ** libs WARNING: no source files found chmod: Zugriff

Re: [R] Interpreting messages when building packages

2008-07-07 Thread Johannes Huesing
This message goes away when I empty the data directory. The problem must be around there. Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : line 12 did not have 19 elements -- Johannes Hüsing There is something fascinating about science.

Re: [R] Another packaging question

2008-07-13 Thread Johannes Huesing
Uwe Ligges [EMAIL PROTECTED] [Sat, Jul 12, 2008 at 11:48:38AM CEST]: Johannes Huesing wrote: I am still trying to build a package. At the moment I am stuck with a file not found error message when processing R code from the tests subdirectory. What would be the accurate relative path

[R] initialize a factor vector

2008-07-13 Thread Johannes Huesing
What is the least surprising way of initializing a factor with predefined levels and with length 0? as.factor(c(eins, zwei, drei))[FALSE] does the job but looks a bit weird. -- Johannes Hüsing There is something fascinating about science. One gets

Re: [R] initialize a factor vector

2008-07-13 Thread Johannes Huesing
Gavin Simpson [EMAIL PROTECTED] [Sun, Jul 13, 2008 at 08:18:37PM CEST]: On Sun, 2008-07-13 at 18:47 +0200, Johannes Huesing wrote: [...] as.factor(c(eins, zwei, drei))[FALSE] does the job but looks a bit weird. [...] factor(levels = c(one,two,three)) factor(0) Levels: one two three

Re: [R] shapiro wilk normality test

2008-07-13 Thread Johannes Huesing
Frank E Harrell Jr [EMAIL PROTECTED] [Sun, Jul 13, 2008 at 08:07:37PM CEST]: (Ted Harding) wrote: On 13-Jul-08 13:29:13, Frank E Harrell Jr wrote: [...] A large P-value means nothing more than needing more data. No conclusion is possible. Please read the classic paper Absence of

Re: [R] shapiro wilk normality test

2008-07-13 Thread Johannes Huesing
Ted Harding [EMAIL PROTECTED] [Sun, Jul 13, 2008 at 10:59:21PM CEST]: On 13-Jul-08 19:53:47, Johannes Huesing wrote: Frank E Harrell Jr [EMAIL PROTECTED] [Sun, Jul 13, 2008 at 08:07:37PM CEST]: (Ted Harding) wrote: On 13-Jul-08 13:29:13, Frank E Harrell Jr wrote: [...] A large P-value

Re: [R] Landscape mode for pdf

2008-07-18 Thread Johannes Huesing
[EMAIL PROTECTED] [EMAIL PROTECTED] [Fri, Jul 18, 2008 at 05:04:15PM CEST]: [...] Generally, this communication is for informational purposes only and it is not intended as an offer or solicitation for the purchase or sale of any financial instrument or as an official confirmation of any

Re: [R] initialize a factor vector

2008-07-19 Thread Johannes Huesing
I have found that factor vectors are shy beasts. Consider: empty - factor(levels=c(eins, zwei)) value - factor(eins, levels=c(eins, zwei)) c(empty, value) [1] 1 empty[1] - value empty [1] eins Levels: eins zwei I could not exactly predict this behaviour, but ?c says: The output type

Re: [R] Discretize continous variables....

2008-07-19 Thread Johannes Huesing
Frank E Harrell Jr [EMAIL PROTECTED] [Sat, Jul 19, 2008 at 08:03:01PM CEST]: But even then why condition on incomplete information when complete information is available? I.e., why compute Pr(Y=1 | Xx) in place of Pr(Y=1 | X=x)? Because regulatory bodies demand it? Being employed in a

Re: [R] Discretize continous variables....

2008-07-20 Thread Johannes Huesing
Frank E Harrell Jr [EMAIL PROTECTED] [Sun, Jul 20, 2008 at 12:20:28AM CEST]: Johannes Huesing wrote: Because regulatory bodies demand it? [...] And how anyway does this relate to predictors in a model? Not at all; you're correct. I was mixing the topic of this discussion up with another

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

2008-07-23 Thread Johannes Huesing
It seems you have accidentally hit a surgeons' mailing list where all you wanted was some advice how to use this scalpel on your body. Sorry if we can't be of any help without intimidating you with unrelated and pompous terms -- like coagulation. Michal J. Figurski [EMAIL PROTECTED] [Wed, Jul 23,

Re: [R] Object-oriented programming in R for Java programmers?

2008-07-27 Thread Johannes Huesing
Werner Wernersen [EMAIL PROTECTED] [Sun, Jul 27, 2008 at 01:55:26PM CEST]: [...] Is there any brief tutorial on oo programming in R especially for people who have done oo in Java or C++ before? That would be really helpful. How S4 methods work highlights the differences between R's

Re: [R] Monty Hall simulation

2008-05-11 Thread Johannes Huesing
cirrus74 [EMAIL PROTECTED] [Sun, May 11, 2008 at 03:44:46AM CEST]: Is it possible to simulate the Monty Hall problem using R? If so, could someone please show me how? Thanks for any help rendered. The kind of simulation, as any thinking about this seeingly paradoxical situation, depends on

Re: [R] Gamma

2009-05-16 Thread Johannes Huesing
Kon Knafelman konk2...@hotmail.com [Sat, May 16, 2009 at 08:49:44AM CEST]: Hey, I am having trouble graphing the following function ???2??(n/2)/[???n ??? 1??((n ??? 1)/2 for the values of n between 2 and 50. i know that ??(n) = (n-1)!, which in R is factorial(n-1) When i type that

Re: [R] How to write loop

2009-06-13 Thread Johannes Huesing
Alex Roy alexroy2...@gmail.com [Sun, Jun 14, 2009 at 06:43:52AM CEST]: Hi Jiim, Thanks . I want to do the following: 1. each time I need to drop one column, say first column 1 from matrix X. 2 then take out row 1 of the remainning matrix and that row becomes response (y) 3. do lasso

[R] libRlapack.so not found

2009-03-21 Thread Johannes Huesing
Whenever I try to load the Matrix package, I get the following error message: libRlapack.so: cannot open shared object file: No such file or directory A file with that name is indeed not on the hard disk. I am using the R version which comes with Ubuntu Hardy Heron LTS. Here is the output of

Re: [R] libRlapack.so not found

2009-03-22 Thread Johannes Huesing
Johannes Huesing johan...@huesing.name [Sat, Mar 21, 2009 at 08:12:05AM CET]: Whenever I try to load the Matrix package, I get the following error message: libRlapack.so: cannot open shared object file: No such file or directory A file with that name is indeed not on the hard disk. http

Re: [R] Selecting closest values

2009-03-22 Thread Johannes Huesing
P_M pmart...@broadpark.no [Sun, Mar 22, 2009 at 06:48:36PM CET]: Hi I have a table with ID (1 to 183) and Location (144 to -22). My problem is that I want to select the 10 ID's that are closest in Location to ID 1, ID 2 and so on. Also, some ID have the same Location. Say, if 11 ID's are

Re: [R] generalized Poisson regression

2009-03-27 Thread Johannes Huesing
Arafat Ud Zaman azam...@isrt.ac.bd [Tue, Mar 24, 2009 at 06:26:59PM CET]: Dear sir, I want to know about R command for parameter estimation for generalized Poisson regression. Do you mean the method of Consul and Shoukri? I once tried Generalized Poisson parameter estimation by their

Re: [R] Selecting elements from a vector (a simple question with weird results)

2009-03-28 Thread Johannes Huesing
Tal Galili tal.gal...@gmail.com [Sat, Mar 28, 2009 at 06:48:36PM CET]: Hello people. I wish to reorder a simple vector of numbers by another vector of the order (and then do the same, but with a data frame rows) I try this (which doesn't work) : aa - c(3, 1 ,2 ) aa[aa] [1] 2 3 1 To

[R] Frailty models and omnibus test

2009-03-29 Thread Johannes Huesing
This is very possibly not a question on R. I was under the impression that the argument that gives rise to Fisher's LSD method in ANOVA works in other situations with three-way comparisons too, given that formal logic works the same (if the omnibus test rejects, only two of the three groups may

Re: [R] Separating variables in read.table

2009-04-17 Thread Johannes Huesing
Richardson, Patrick patrick.richard...@vai.org [Fri, Apr 17, 2009 at 09:31:35PM CEST]: If I have a table (we'll call it, test) containing seven columns (as below): i x1 x2 x3 x4 x5 y 0 1.125 0.232 7.160 0.0859 8.905 1.5563 7 0.920 0.268 8.804

Re: [R] Binomial simulation

2009-04-17 Thread Johannes Huesing
beetle2 samandbren...@aapt.net.au [Fri, Apr 17, 2009 at 11:28:56PM CEST]: Hi Guy's I was wondering if someone could point me in the right direction. dbinom(10,1,0.25) I am using dbinom(10,1,0.25) to calculate the probabilty of 10 judges choosing a certain brand x times. dbinom returns

Re: [R] Binomial simulation

2009-04-19 Thread Johannes Huesing
beetle2 samandbren...@aapt.net.au [Sun, Apr 19, 2009 at 03:01:43AM CEST]: I'm thinking I will just use: results - rbinom(1000, 10, .25) d = sum(results == 0 ) df = (d/1000) df And do each individually ok, at least you are trying some things out for yourself. In fact, along with reading

Re: [R] get() versus getAnywhere()

2009-04-19 Thread Johannes Huesing
Rolf Turner r.tur...@auckland.ac.nz [Sun, Apr 19, 2009 at 10:01:34PM CEST]: On 19/04/2009, at 9:45 PM, Duncan Murdoch wrote: On 18/04/2009 8:47 PM, Rolf Turner wrote: On 17/04/2009, at 10:21 PM, Duncan Murdoch wrote: Benjamin Tyner wrote: Many thanks Duncan. Perhaps this merits a more

Re: [R] (no subject)

2008-11-18 Thread Johannes Huesing
Salas, Andria Kay [EMAIL PROTECTED] [Tue, Nov 18, 2008 at 04:10:06PM CET]: I need help with another problem I am having that deals with the generation of vectors that I asked about yesterday. I now need to have each value in the vector (all values either 1 or -1) have a probability p that it

Re: [R] Welcome to the R-help mailing list

2009-01-29 Thread Johannes Huesing
Bishwa S Koirala bkoir...@unm.edu [Thu, Jan 29, 2009 at 08:56:17PM CET]: You must know your password to change your options (including changing the password, itself) or to unsubscribe. It is: It is definitely a good idea to change your password ASAP. -- Johannes Hüsing There

Re: [R] Problems in Recommending R

2009-02-12 Thread Johannes Huesing
Andrew Choens andy.cho...@gmail.com [Thu, Feb 12, 2009 at 04:56:33PM CET]: * Quick-R - http://www.statmethods.net/ * The R Inferno - www.burns-stat.com/pages/Tutor/R_inferno.pdf * Rseek.org As an aside: Last time I tried, rseek.org yielded no results when searching

Re: [R] Adjusting confidence intervals for paired t-tests of multiple endpoints

2009-02-28 Thread Johannes Huesing
Erich Studerus, Psychiatrische Uni-Klinik erich.stude...@bli.uzh.ch [Sat, Feb 28, 2009 at 09:23:05PM CET]: Thanks, but as far as I know the multcomp package and the tukeyHSD function can only handle one dependent variable, whereas I have 60 dependent variables. Couldn't you re-phrase your

Re: [R] R on netbooks et al?

2009-03-06 Thread Johannes Huesing
chaogai chao...@xs4all.nl [Thu, Mar 05, 2009 at 07:04:19PM CET]: I'm having similar experiences on my Acer Aspire One. Everything will work good. Only thing that takes a lot of time is compiling R if you are in the habit of doing so. On the Fedora version that came with my Acer Aspire One, I

Re: [R] standard error for median

2009-03-07 Thread Johannes Huesing
Ralph Scherer scherer...@googlemail.com [Sat, Mar 07, 2009 at 10:34:28PM CET]: Dear Prof. Ripley, thank you for your fast answer. But which book do you mean? I can't find an MASS book. Try library(MASS) citation(package=MASS) -- Johannes Hüsing There is something

Re: [R] Examples of advanced data visualization

2008-11-29 Thread Johannes Huesing
Tom Backer Johnsen [EMAIL PROTECTED] [Fri, Nov 28, 2008 at 10:13:04PM CET]: Hans W. Borchers wrote: [...] Please answer to my e-mail address. In case enough interesting material comes up, I will enter a summary here. It is nice that you are willing to summarize whatever appears, but

[R] Formula parsing and updating

2008-12-04 Thread Johannes Huesing
Hi all, I can't come over a problem with formula. Suppose I have a coxmod model with the following formula: somemod$formula Surv(lebzeit, tot == 1) ~ sex + (alter = 65) + diff3k + zelltyp_k_c + q_nuc_3k + kar_k80_g80 + stadium and I want to drop the stadium explanatory variable from the

Re: [R] Formula parsing and updating

2008-12-05 Thread Johannes Huesing
Charles C. Berry [EMAIL PROTECTED] [Fri, Dec 05, 2008 at 01:21:41AM CET]: On Thu, 4 Dec 2008, Johannes Huesing wrote: somemod$formula Surv(lebzeit, tot == 1) ~ sex + (alter = 65) + diff3k + zelltyp_k_c + q_nuc_3k + kar_k80_g80 + stadium You need to give us something to reproduce. I get

Re: [R] R and Scheme

2008-12-09 Thread Johannes Huesing
Stavros Macrakis [EMAIL PROTECTED] [Wed, Dec 10, 2008 at 04:59:25AM CET]: So I conclude that what is really meant by R semantics are based on Scheme semantics is R has functions as first-class citizens and a correct implementation of lexical scope, including upwards funarg. One other thing

Re: [R] ATT Researchers and the New York Times

2009-01-08 Thread Johannes Huesing
stephen sefick ssef...@gmail.com [Fri, Jan 09, 2009 at 03:16:11AM CET]: It has worked wonders for me over the last years. On Thu, Jan 8, 2009 at 8:58 PM, Robert Wilkins irishhac...@gmail.com wrote: And by the way, ARE YOU GUYS EVER GOING TO FIX your mailing list platform? It is extremely

Re: [R] R in the NY Times

2009-01-10 Thread Johannes Huesing
Bert Gunter gunter.ber...@gene.com [Sat, Jan 10, 2009 at 08:31:03PM CET]: [...] Perhaps an argument is that certain code might not get written at all if it were not proprietary. Device drivers might be an example. Device drivers are not an example. Linux is ubiquitous _despite_ device

Re: [R] Best R text editors?

2009-09-10 Thread Johannes Huesing
Martin Maechler maech...@stat.math.ethz.ch [Wed, Sep 02, 2009 at 09:17:42AM CEST]: PaCo == p connolly p_conno...@slingshot.co.nz on Wed, 02 Sep 2009 12:19:31 +1200 writes: PaCo On Mon, 31-Aug-2009 at 08:25PM +1000, Jim Lemon wrote: PaCo [...] [...] PaCo | Emacs still

Re: [R] two questions for R beginners

2010-02-27 Thread Johannes Huesing
Dieter Menne dieter.me...@menne-biomed.de [Fri, Feb 26, 2010 at 08:39:14AM CET]: Patrick Burns wrote: * What were your biggest misconceptions or stumbling blocks to getting up and running with R? (This derives partly from teaching) [...] The concept of environment. With

[R] OT: Import HTML, was Re: How good is R at making publication quality tables?

2010-03-18 Thread Johannes Huesing
Jim Lemon j...@bitwrit.com.au [Wed, Mar 17, 2010 at 11:43:21PM CET]: [...] All of the word processing packages that I use can directly import HTML documents, and the result has been good enough for any application so far. I have only MS Word as work and failed to import HTML so far

Re: [R] Coding of categorical variables for logistic regression?

2010-03-28 Thread Johannes Huesing
Ravi Kulkarni ravi.k...@gmail.com [Sun, Mar 28, 2010 at 07:15:52PM CEST]: [...] I am trying to do a logistic regression and have one predictor variable (x) that is ratio and two predictor variables (y and z) that are categorical. These have three levels each which I have called High, Medium

Re: [R] Confusing concept of vector and matrix in R

2010-04-01 Thread Johannes Huesing
Rolf Turner r.tur...@auckland.ac.nz [Tue, Mar 30, 2010 at 09:48:12PM CEST]: [...] Who designed Excel? Among others, Joel Spolsky. Is that an appeal to authority? http://www.joelonsoftware.com/articles/fog20.html -- Johannes Hüsing There is something fascinating about

Re: [R] R as a non-functional language

2011-03-20 Thread Johannes Huesing
Russ Abbott russ.abb...@gmail.com [Sun, Mar 20, 2011 at 06:46:04AM CET]: I'm afraid I disagree. I don't think there is much disagreement. You have been shown that there is a functional way to rephrase your example. I suppose you don't disagree that the structure() expression has no side

Re: [R] Problem applying McNemar's - Different values in SPSS and R

2010-12-28 Thread Johannes Huesing
Marc Schwartz marc_schwa...@me.com [Tue, Dec 28, 2010 at 06:30:59PM CET]: On Dec 28, 2010, at 11:05 AM, Manoj Aravind wrote: Hi friends, I get different values for McNemar's test in R and SPSS. Which one should i rely on when the p values differ. [...] The SPSS test appears to be

Re: [R] Problem applying McNemar's - Different values in SPSS and R

2010-12-28 Thread Johannes Huesing
Marc Schwartz marc_schwa...@me.com [Tue, Dec 28, 2010 at 07:14:49PM CET]: [...] An old question of mine: Is there any reason not to use binom.test() other than historical reasons? (I meant in lieu of the McNemar approximation, sorry if some misunderstanding ensued). I may be missing the

Re: [R] Problem applying McNemar's - Different values in SPSS and R

2010-12-29 Thread Johannes Huesing
Marc Schwartz marc_schwa...@me.com [Wed, Dec 29, 2010 at 03:28:56PM CET]: On Dec 29, 2010, at 6:48 AM, Manoj Aravind wrote: Thank you Marc :) It Certainly helped me to get the exact value of P. How to understand when to apply mcnemar.exact or just mcnemar.test? [...] Generally

Re: [R] Problem applying Chi-square in R and Cochran's Recommendations

2010-12-29 Thread Johannes Huesing
Manoj Aravind aravin...@gmail.com [Wed, Dec 29, 2010 at 03:59:16PM CET]: Sir, I have a problem here while applying chisquare test to the following Data ( below the subject of this mail) ...when I wanted to test the significance using three different free statistical packages, here R, EpiInfo

[R] Meaning of pterms in survreg object?

2011-01-10 Thread Johannes Huesing
I am trying to model survival data with a Weibull distribution using survreg. Units are clustered two apiece, sometimes receiving the same treatment and

Re: [R] Meaning of pterms in survreg object?

2011-01-27 Thread Johannes Huesing
Is there a special mailing list for the survival package where my question would have been more appropriate? Johannes Huesing johan...@huesing.name [Mon, Jan 10, 2011 at 09:26:37PM CET]: I am trying to model survival data with a Weibull distribution

Re: [R] Meaning of pterms in survreg object?

2011-01-27 Thread Johannes Huesing
David Winsemius dwinsem...@comcast.net [Thu, Jan 27, 2011 at 10:08:00PM CET]: You got a perfectly sensible reply from Thereau, the author of the package, a day after your posting and then failed to respond to his questions. I'm not sure what more you are expecting. My apologies. Next time

Re: [R] (no subject)

2011-06-02 Thread Johannes Huesing
Ungku Akashah kasla...@yahoo.com [Thu, Jun 02, 2011 at 10:31:53AM CEST]: Error in lda.default(x, grouping, ...) : variables 1 3 5 8 10 15 17 20 27 29 34 appear to be constant within groups Please do read the posting guide before posing a question here. The message would have benefitted

Re: [R] Meaning of pterms in survreg object?

2011-01-28 Thread Johannes Huesing
Johannes Huesing johan...@huesing.name [Mon, Jan 10, 2011 at 09:26:37PM CET]: I am trying to model survival data with a Weibull distribution using survreg. Units are clustered two apiece, sometimes receiving

Re: [R] Data manipulation

2011-02-12 Thread Johannes Huesing
mathijsdevaan mathijsdev...@gmail.com [Sat, Feb 12, 2011 at 03:00:18PM CET]: Hi, I have a dataset with info on individuals (B) that have been involved in projects (A) during multiple years (C). The dataset contains three columns: A, B, C. Example: A B C 1 1 a 1999 2 1 b 1999 3

[R] R coredumps when calling functions from rms package

2014-10-06 Thread Johannes Huesing
My wife is trying to run some bootstrap model validation. She is using code to the effect of library(rms) residual - sasxport.get(~/R_boot/residual.xpt) c1.odat - somers2(residual$risk.5, residual$caco)[C]# c - statistic for original risk estimate in full sample # fitboot -

Re: [R] Checking if a matrix exists/is defined

2014-10-19 Thread Johannes Huesing
Steven Yen sye...@gmail.com [Sun, Oct 19, 2014 at 06:17:48AM CEST]: Hello Can someone help me with the following, specifically in judging whether a matrix exists. exists(my.matrix) I have trouble with the first line below. In this case, matrix obj$hessian exists and is 74 x 74. I receive

Re: [R] Counting consecutive events in R

2015-05-14 Thread Johannes Huesing
I normally use rle() for these problems, see ?rle. for instance, k - rbinom(999, 1, .5) series - function(run) {

[R] problem understanding grid coordinate systems

2015-07-06 Thread Johannes Huesing
According to R Graphics by Paul Murrell, the coordinates that were used for the most recent lattics plot can be retrieved with native units. I have difficulties to access these coordinates. The following code renders the following results: library(grid) library(lattice) pl - xyplot(1:10 ~

Re: [R] Inquiry: R library for bioinformatics development process

2015-07-28 Thread Johannes Huesing
Faris Zaidi Bin Mohd Nor farisza...@mardi.gov.my [Wed, Jul 29, 2015 at 05:38:59AM CEST]: Hi, I would like to know about the R library. Is there any library for bioinformatic development software purpose? Have you looked at the Bioconductor project? -- Johannes Hüsing There

[R] problem understanding grid coordinate systems

2015-07-13 Thread Johannes Huesing
According to R Graphics by Paul Murrell, the coordinates that were used for the most recent lattics plot can be retrieved with native units. I have difficulties to access these coordinates. The following code renders the following results: library(grid) library(lattice) pl - xyplot(1:10 ~

Re: [R] problem understanding grid coordinate systems

2015-07-15 Thread Johannes Huesing
Paul Murrell p...@stat.auckland.ac.nz [Wed, Jul 15, 2015 at 11:12:33PM CEST]: ... downViewport(plot_01.panel.1.1.vp) ... This works like a charm. Thank you! Upon reading more of the documentation, I found that current.vpTree() shows all names of active viewports. -- Johannes Hüsing

Re: [R] no results

2015-11-10 Thread Johannes Huesing
Alaa Sindi [Tue, Nov 10, 2015 at 08:47:23PM CET]: Hi All, I am not getting any summary results and I do not have any error. what would be the problem? sem=mlogit.optim ( LL , Start, method = 'nr', iterlim = 2000, tol = 1E-05, ftol = 1e-08, steptol = 1e-10,

[R] problem understanding grid coordinate systems

2015-07-08 Thread Johannes Huesing
According to R Graphics by Paul Murrell, the coordinates that were used for the most recent lattics plot can be retrieved with native units. I have difficulties to access these coordinates. The following code renders the following results: library(grid) library(lattice) pl - xyplot(1:10 ~

Re: [R] Rayleigh Distribution

2015-10-02 Thread Johannes Huesing
Jomy Jose [Fri, Oct 02, 2015 at 05:45:58PM CEST]: Thank you,how to add correct reference line to this Rayleigh Q-Q plot. ?abline -- Johannes Hüsing http://derwisch.wikidot.com Threema-ID: VHVJYH3H __

Re: [R] Randomness tests

2015-09-24 Thread Johannes Huesing
Giorgio Garziano [Tue, Sep 22, 2015 at 10:24:25PM CEST]: Hi, to test randomness of time series whose values can only be +1 and -1, are all following randomness tests applicable or only a part of ? I don't know the details of all those tests but the general

Re: [R] Extracing Unique Rows based on 2 Column

2015-11-29 Thread Johannes Huesing
Ragia Ibrahim [Mon, Nov 30, 2015 at 04:55:08AM CET]: Dear group, kindly,  I have a data frame, as follows:  Measure_id i j value      rank I want to select distinct rows based on two coulmn ( Measure_id  and i ) I didn't get your example code to run but the following

  1   2   >