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 works for me: dfr <-

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, print.level = 0) I see som

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 __ R-help@r-project.org mailing lis

Re: [R] Randomness tests

2015-09-23 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 problem is that the alternative

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

2015-07-28 Thread Johannes Huesing
Faris Zaidi Bin Mohd Nor [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 is something fascinating

Re: [R] problem understanding grid coordinate systems

2015-07-15 Thread Johannes Huesing
Paul Murrell [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 There is somethin

[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

[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

[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] 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) {

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

2014-10-18 Thread Johannes Huesing
Steven Yen [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 the > erro

[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 <- fu

Re: [R] (no subject)

2011-06-02 Thread Johannes Huesing
Ungku Akashah [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 from a title. M

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

2011-03-20 Thread Johannes Huesing
Russ Abbott [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 effects. And hardly anyone

Re: [R] Data manipulation

2011-02-12 Thread Johannes Huesing
mathijsdevaan [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 1 c 1999 >

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

2011-01-28 Thread Johannes Huesing
Johannes Huesing [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

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

2011-01-27 Thread Johannes Huesing
David Winsemius [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 I'll be checking Gm

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 [Mon, Jan 10, 2011 at 09:26:37PM CET]: > I am trying to model survival data with a Weibull distribut

[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 so

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

2010-12-29 Thread Johannes Huesing
Manoj Aravind [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 and > OpenEpi.

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

2010-12-29 Thread Johannes Huesing
Marc Schwartz [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 speaking, exac

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

2010-12-28 Thread Johannes Huesing
Marc Schwartz [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 context of yo

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

2010-12-28 Thread Johannes Huesing
Marc Schwartz [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 an exact

Re: [R] diverting output from nested loops

2010-11-24 Thread Johannes Huesing
Lara Poplarski [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] Sweave question

2010-11-13 Thread Johannes Huesing
Ralf B [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 http://www.statistik.lmu.de/~leisch/Sweave/Sweave-Rnews-2002-3.

Re: [R] setting attributes

2010-11-04 Thread Johannes Huesing
David Winsemius [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 it is meant to do. -

Re: [R] stacking effectively - without loops

2010-11-02 Thread Johannes Huesing
Dimitri Liakhovitski [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 > a 2 > a 3 >

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

2010-10-31 Thread Johannes Huesing
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. > First sect

[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 i

Re: [R] Mapping the coordinates!

2010-10-09 Thread Johannes Huesing
Mehdi Zarrei [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, longitude

Re: [R] Non-Parametric Adventures in R

2010-10-03 Thread Johannes Huesing
Jamesp [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 out now. summary(X

Re: [R] Sweave and LaTeX beamer class

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

[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 \uncover<2->{25}\uncover<3->{\Sexpr{5*5}} \uncover<4->{ <>= 5*5*101 @ } \end{frame} \

Re: [R] Combining two variables in text

2010-09-28 Thread Johannes Huesing
trb1 [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 separated

Re: [R] Statistical mailing list

2010-08-14 Thread Johannes Huesing
Tal Galili [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 There is so

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

2010-06-28 Thread Johannes Huesing
Johannes Huesing [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

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

2010-06-28 Thread Johannes Huesing
GL [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) dbs1$Payor.Group[i] = >

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

2010-06-22 Thread Johannes Huesing
Jim Hargreaves [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, 6, 7, 8,

Re: [R] Read code from character string

2010-06-20 Thread Johannes Huesing
Gabor Grothendieck [Sat, Jun 19, 2010 at 11:19:35AM CEST]: > On Sat, Jun 19, 2010 at 4:56 AM, Johannes Huesing > wrote: > > A rule would be a data structure containing the ID of the rule, the > > rule in human readable language, an expression evaluating variables > > w

Re: [R] Read code from character string

2010-06-19 Thread Johannes Huesing
Greg Snow [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. I could not pos

Re: [R] Read code from character string

2010-06-18 Thread Johannes Huesing
Greg Snow [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 frowned upon in R

Re: [R] Read code from character string

2010-06-17 Thread Johannes Huesing
Peter Langfelder [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 > wrote: > > Dear expRts, > > I have a character string, say a <- &q

[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 sci

Re: [R] Sampling with replacement

2010-06-16 Thread Johannes Huesing
Somnath Somnath [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? m^n -- Joh

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 [Fri, May 21, 2010 at 04:47:21PM CEST]: > Dear R users, > # > #There is table containing 1000 (lat, lon, time) GPS positions, wchich should > be recognized

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

2010-05-21 Thread Johannes Huesing
Jan van der Laan [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 (reading line by line) sound

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

2010-04-01 Thread Johannes Huesing
Rolf Turner [Tue, Mar 30, 2010 at 09:48:12PM CEST]: [...] > Who designed Excel? Among others, Joel Spolsky. Is that an appeal to authority? -- Johannes Hüsing There is something fascinating about science.

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

2010-03-28 Thread Johannes Huesing
Ravi Kulkarni [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" and "Low".

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

2010-03-18 Thread Johannes Huesing
Jim Lemon [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 (sometimes succeeded wi

Re: [R] two questions for R beginners

2010-02-27 Thread Johannes Huesing
Dieter Menne [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 S it was w

Re: [R] Novel (Maybe?) Visualizations

2010-02-10 Thread Johannes Huesing
D. Dashcle [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 approach in mind, just look

Re: [R] Learning R

2009-11-30 Thread Johannes Huesing
Esmail [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 judging by the author's

Re: [R] Best R text editors?

2009-09-10 Thread Johannes Huesing
Martin Maechler [Wed, Sep 02, 2009 at 09:17:42AM CEST]: > > "PaCo" == p connolly > > 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 > PaCo> |> has that annoying tr

Re: [R] Reading from Google Docs

2009-07-09 Thread Johannes Huesing
Gabor Grothendieck [Wed, Jul 08, 2009 at 10:31:30PM CEST]: [...] > >>> > >>>    [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.et

Re: [R] How to write loop

2009-06-13 Thread Johannes Huesing
Alex Roy [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 regression on remain

Re: [R] Gamma

2009-05-16 Thread Johannes Huesing
Kon Knafelman [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 into R, usin

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

2009-04-19 Thread Johannes Huesing
Rolf Turner [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 explicit

Re: [R] Binomial simulation

2009-04-19 Thread Johannes Huesing
beetle2 [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 the doc, this is

Re: [R] Binomial simulation

2009-04-17 Thread Johannes Huesing
beetle2 [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 the discrete dens

Re: [R] Separating variables in read.table

2009-04-17 Thread Johannes Huesing
Richardson, Patrick [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 0.0865 7.388 0.8976

[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 b

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

2009-03-28 Thread Johannes Huesing
Tal Galili [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 my mind, i

Re: [R] generalized Poisson regression

2009-03-27 Thread Johannes Huesing
Arafat Ud Zaman [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 suggestion side by si

Re: [R] Selecting closest values

2009-03-22 Thread Johannes Huesing
P_M [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 closest to ID 100

Re: [R] libRlapack.so not found

2009-03-22 Thread Johannes Huesing
Johannes Huesing [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://l

[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 R.Ver

Re: [R] standard error for median

2009-03-07 Thread Johannes Huesing
Ralph Scherer [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 fascinating about science.

Re: [R] R on netbooks et al?

2009-03-06 Thread Johannes Huesing
chaogai [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 am even think

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

2009-02-28 Thread Johannes Huesing
Erich Studerus, Psychiatrische Uni-Klinik [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 model by including timep

Re: [R] Problems in Recommending R

2009-02-12 Thread Johannes Huesing
Andrew Choens [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 for "inferno". --

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

2009-01-29 Thread Johannes Huesing
Bishwa S Koirala [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 is somethin

Re: [R] R in the NY Times

2009-01-10 Thread Johannes Huesing
Bert Gunter [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 manufacturers being secretiv

Re: [R] AT&T Researchers and the New York Times

2009-01-08 Thread Johannes Huesing
stephen sefick [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 wrote: > > And by the way, ARE YOU GUYS EVER GOING TO FIX your mailing list > > platform? It is extremely user-unfriendly and a technologic

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 othe

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 >

[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 m

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, bu

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 tha

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 is.na(

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 funct

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] 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 top

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 | X>x) in place of Pr(Y=1 | X=x)? Because regulatory bodies demand it? Being employed in

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:

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] 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: > >>>

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

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. > > [..

[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

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. Wh

[R] Another packaging question

2008-07-12 Thread Johannes Huesing
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 for files in the tests directory to access files in the data directory? -- Johannes Hüsing T

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] 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: Zugrif

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. -- Joh

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

[R] Interpreting messages when building packages

2008-07-06 Thread Johannes Huesing
Dear expRts, when I try to build a package by myself, the process fails writing the following messages to 00install.out: * Installing *source* package 'dynalc' ... ** libs WARNING: no source files found chmod: Zugriff auf „/home/hannes/texte/forschung/clot/programme/dynalc.Rcheck/dynalc/libs/*“ n

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 >

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 D

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: > > li

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 d

  1   2   >