[R] aggregate function with 'NA'

reciate any help. Yours Frank [[alternative text/enriched version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and pr

Re: [R] question on boot

boot will handle multivariate data as well as multistratum data, though it isn't clear from your question what is multivariate about your problem or what you are bootstrapping. Lisa Wang wrote: > Hello there: > > Dear Dr. Murdoch, > > I'm a statistician at Princess Margaret Hospital. Could you

[R] possible bug in merge with duplicate blank names in 'by' field.

h I expected substantial overlap in the labels, were taking so long to merge and ultimately generating 57000 labels. The culprit turned out to be a few hundred blanks buried in the middle. Why does the empty ("null") name merit special treatment? Perhaps I'm missing something. I

Re: [R] possible bug in merge with duplicate blank names in 'by' field.

I'll upgrade my Linux installation to the most recent version (1.9.1 is a year old, I guess). -Frank At 03:26 AM 6/17/2005, Prof Brian Ripley wrote: >What version of R is this (please do see the posting guide)? > >In both 2.1.0 and 2.1.1 beta I get > >>all > Promote

[R] Trouble building R2.1.0 from source on Linux: package VR

ith merge when there are empty labels), any assistance greatly appreciated. Thanks in advance. -Frank PhD, Computational Biologist, Harvard Medical School BCMP/SGM-322, 250 Longwood Ave, Boston MA 02115, USA. Tel: 617-432-3555 Fax: 617-432-3557 http://llama.med.harvard.edu/~fgibbo

Re: [R] Binary outcome with non-absorbing outcome state

Singer & Willett (2003) also cover this ground. Singer, JD & Willett, JB (2003). Applied longitudinal data analysis: Modeling change and event occurrence. New Yok: Oxford University Press. -Original Message- From: Frank E Harrell Jr <[EMAIL PROTECTED]> Sent: Jul 29, 2

[R] fitting a model with the nlme package

3,1,2,3,1,2,3) f<-c(1,2,1,2,1,2,1,2,1) p<-c(45,32,35,12,23,12,2,9,12) set<-data.frame(p,h,j,t) out<-lme(p~ -1 + f,data=set, random=~1|h/t) Thanks a alot, frank. -- unladen european swallow __ R-help@stat.math.ethz.ch mailing li

[R] optim error

attempt to fix it? Starting values? Something in the log-likelihood itself? Perhaps the data (which is generated)? Any thoughts would be greatly appreciated. Thanks, Frank __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listi

[R] scale function within a for-loop?

I have a data matrix containing around 180 variables and more than 470 observations for each and no missing values. Within a for-loop, the first step of calculations is to standardize each column, such that the mean of each column is zero and the sd is one. The for-loop starts with a subset of t

[R] R: extracting elements in a matrix

Dear R-users For a given matrix of dimension, say (n,p), I'd like to extract for every column those elements that are bigger than twice the interquartile range of the corresponding column. Can I get these elements without using a loop? Thank you for your help

[R] Save output

3 dimensional array and saved it just once after the loop. Or is there another way so that I can keep the resulting matrix of every step in the loop? Thank you very much for your help Frank Schmid [[alternative HTML version deleted]] __ R-he

[R] subset by multiple columns satisfying the same condition

Hi All, I have a very simple question. Suppose I had a data frame with 100 columns, now I wanted to select rows with the values of some columns satisfying the same condition, like all equal to "Tom". I know I can use the 'and' operator "&", but it's painful if there were many columns. Can anyone

[R] snow parLapply standard output

al where I invoked it? I intend to transmit some status information in advance to the final result of the function. I investigated the chain of functions called by parLapply and it seems that snow is designed to just retrieve the final result of the computation. Thanks in advance for any hints, Frank

[R] impossible to open SPSS file

eur dans read.spss("C:RdataESS1_ICT.SAV") : impossible d'ouvrir le fichier I guess the error is evident for everyone else, but I don't see it. Thanks Frank Thomas -- .. Dr. Frank Thomas FTR Internet

Re: [R] impossible to open SPSS file

possible to > Excel) and import it to R > SPSS does not externalize the data format :-( > Regards Knut > We don't have SPSS so unfortunately this work around is not possible. Thank you both, Frank -- ...... Dr. Frank Thomas FTR Internet R

[R] Splitting a character variable into a numeric one and a character one?

Hi All, I have a data with a variable like this: Column 1 "123abc" "12cd34" "1e23" ... Now I want to do an operation that can split it into two variables: Column 1Column 2 Column 3 "123abc" 123 "abc" "12cd34" 12"cd34" "1e23"

Re: [R] Splitting a character variable into a numeric one and a character one?

<- strapply(s, "^([[:digit:]]+)(.*)", c) > out <- do.call(rbind, out) # as a matrix > > data.frame(x = out[,1], num = as.numeric(out[,2]), char = out[,3]) # > as a data.frame > > > On 9/25/06, Frank Duan <[EMAIL PROTECTED]> wrote: > > Hi All, > > &

Re: [R] [R-pkgs] New versions of Matrix and lme4 packages for R-2.4.0

Thanks, Douglas. It works. Only one gripe. Every time I want to try a new version of your software, I have to get, compile, patch the latest bleeding edge version of R as well: ERROR: This R is version 2.3.1 package 'Matrix' needs R >= 2.4.0 Douglas Bates wrote: > Versions 0.9975-1 o

[R] extracting nested variances from lme4 model

the cases. How do I get the other variance estimate? Thanks for any help. -Frank My fake data set: --- data: $frame x rtr trth cs 1 -4.4964808 1 TRUE 1 2 -0.6297254 1 TRUE 2 31.8062857 1 TRUE 3 42.7273275 1 TRUE 4

[R] Install RMySQL with R 2.4.0

I just installed RMySQL 0.5-9 with R 2.4.0 on Windows XP and got the following error message when trying to run a script with RMySQL: Error in library(RMySQL) : 'RMySQL' is not a valid package -- installed < 2.0.0? Any ideas? Thanks, Frank ___

[R] Misunderstanding escape (backslash)

other\\\'s day" > term [1] "mother\\'s day" <--- should be "mother\'s day" Thanks, Frank __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] Misunderstanding escape (backslash)

> And check out the new FAQ 7.37 > > http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-does-backslash-behave-strangely-inside-strings_003f > Thanks, guys, for the help. Frank __ R-help@stat.math.ethz.ch mailing list https://stat.eth

[R] Turn off border on barplot

with only the bar's color visible. Thanks, Frank -- Frank McCown Old Dominion University http://www.cs.odu.edu/~fmccown/ __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting gui

Re: [R] Turn off border on barplot

> barplot(1:10, border=NA) That's it! Thanks a ton, guys. Frank __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide c

[R] Adding to user-contributed documentation

ted called "Producing Simple Graphs with R" http://www.cs.odu.edu/~fmccown/R/ I plan on adding additional graphs (like scatterplots) as time permits. Thanks, Frank __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-hel

Re: [R] Install RMySQL with R 2.4.0

to be no longer valid. Thanks, Frank -- Frank McCown Old Dominion University http://www.cs.odu.edu/~fmccown/ __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/p

Re: [R] Install RMySQL with R 2.4.0

for others to > download. > > Good Luck > Joe Joe, Thanks for telling us how you got RMySQL installed. Would you mind posting the dll files so the rest of us wouldn't have to recompile anything? Thanks, Frank __ R-help@stat.math.et

[R] lme4 package: Fitted values and residuals

w, I would like to calculate a fitted value, 'yhat', from only the parameter estimates for 'a' and 'b|c'. Can this be done in lme4. That is; can I pull out all of the requisite elements from the output to calculate 'yhat'? Thanks for your help, Frank. --

[R] loop problem

the "print" or "cat" command, but could not turn the print or cat results into an accessible object. Help would be appreciated. Thank you, Frank. -- or over the web __ R-help@stat.math.ethz.ch mailing lis

[R] help with simple 3d graph

try commands like "persp" it wants ordered values. The only command that approximates what I want is "scatterplot3d" (see below) library(scatterplot3d) scatterplot3d(a,b,c) Can anybody help me with this. Thanks very much. Frank. -- _

[R] tapply question

HI, Suppose I have the following data structure. LRT tp 1 1.50654010 522 2 0.51793929 522 3 0.90340299 522 4 1.20293325 522 5 1.05578774 523 6 0.01617942 523 7 0.68183543 523 8 0.43820244 523 9 1.14123995 524 10 0.05809550 524 11 0.93061597 524 12 1.39739700 524 13

[R] Looking for packages to do Feature Selection and Classification

100 useful features or 2) using all 200 features. Any suggestions or point me to a reference? Thanks in advance! Frank [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-hel

Re: [R] Looking for packages to do Feature Selection and Classification

lowing paper on svm (using libsvm) as well as random > forest in the context of feature selection. > > http://www.csie.ntu.edu.tw/~cjlin/papers/features.pdf > > HTH > > On 1/4/06, Diaz.Ramon <[EMAIL PROTECTED]> wrote: > > Dear Frank, > > I expect you'

[R] Newton-Raphson algorithm

, available in R? If not, does anybody have an idea how to best go about solving simultenous equations numerically in R? Thanks for your help. Frank. -- love email again __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch

Re: [R] problem to install R on linux

For the rpm install you may also need to install a package called info-x.x.x.x-x.rpm. For the compilation you'll need to be certain that you've installed readline-x.x.x.x.x.rpm and and perhaps readline-devel-x.x.x.x.rpm Those should be in your mandriva distn André Beló wrote: > Dear members, >

[R] tk mouse cursor icon widget tkwinfo tkfocus questions

messy. 2. How do I specify the force option with tkfocus? Or, more generally, any option in a tk function that doesn't take a value? I can see the force option is available from the errors at my failed attempts to specify it. MS keeps popping up my Tk window behind an

Re: [R] using kernel density estimates to infer mode of distribution

A density() fit calls the eval x and estimate y: fit<-density(data) plot(fit$x,fit$y) Adelchi Azzalini wrote: > On Wed, 15 Feb 2006 18:28:25 -0500, Dan Rabosky wrote: > > DR> > DR> Is it possible to use "density" or another kernel density > DR> estimator to identify the mode of a distributi

[R] exces return by mktcap decile for each year

I have a data frame, lets call it dat, with 3 columns ( mc, yr, ret) which represent market cap, year, and return. mc is a factor, mc, and ret are real numbers. I want to add a column to the data calculated as follows. For each year, I want to split the data by mc decile, then calculate the mean

Re: [R] Ordination of feature film data question

[[alternative HTML version deleted]] >> >> __ >> R-help@stat.math.ethz.ch mailing list >> https://stat.ethz.ch/mailman/listinfo/r-help >> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] How to compare areas under ROC curves calculated with ROCR package

d in performing hypothesis testing on the difference in AUC of two parameters, I would suggest a permutation test. Permuting the ranks of the data between parameters is simple and works well. -Frank ## AuROC<-function(neg,pos) { #empiri

[R] Branching, depending upon operating system

I cant get it to work. It really shouldn`t be difficult, but how can I do an if-branch using that (or a different) command? -- Regards, Frank Rossdeutscher [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https:/

[R] Feeding a sequence to a function

Dear all, I have written a function that takes two arguments, and I would like to feed it all pairs of values from -10 to 10. The following code works for all pairs of values between 1 and 10, but given R's indexing, I can't extend it back to cover the zeros and the negative values. I'd appre

[R] a clipboard problem while using R2HTML

Error: clipboard connection is open for reading only The only setting about the clipboard is in my .emacs is: '(x-select-enable-clipboard t)) Thank you in advance for any hint. I am using Linux Fedora 4. Kernel 2.6.15-1.1833_1. R version 2.3.0, ESS 5.3.0. -Frank _

Re: [R] Net courses for R?

gement tasks in R, including how to get data into R, how to recode variable and how to subset data. The R code and the output produced are shown, as well as tips on things to look out for. Cheers, Frank Thomas Scott Cunningham wrote: > Are there any online courses

Re: [R] R pour débutant

__ > Denis POIGNONEC > http://denis.poignonec.free.fr > > IRD CoRéUs > BP A5 > 98848 Nouméa Cedex > Nouvelle Calédonie > mob: +(687) 92.09.17 > > __ > R-help@stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman

[R] installing SciViews

an install process. I guess I miss something. Cheers, F. Thomas -- .. Dr. Frank Thomas FTR Internet Research 93110 Rosny-sous-Bois France -- __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/

Re: [R] How-Understand-Kmeans-Cluster!!

ustering > from the datset (DDset)? How can i understand the *similarities and > dis-similarities*? Will Kmeans-Clustering > helps on this problem? Where can i get *more tutorials* on > kmeans-clustering? > > Thanks for reading this email. > > Thanks in Advance > > JJ

Re: [R] all.equal() and which()

(abs(A-B)) and have it in a useful numeric form. -Frank Marc Schwartz (via MN) wrote: > On Fri, 2006-02-03 at 10:41 -0500, tom wright wrote: > >>Please excuse the lack of a complete dataset here, if its needed I'll be >>happy to provide it. >>Can anyone show me how t

Re: [R] executable file with R

Romain Lorrilliere wrote: > Hi, > > I made an R function, and I want make an executable applet with it. Do > you know how it is possible? > More details about what you want would be helpful. Here is what I do and it may be useful to you. Under *nix,OSX, etc use bash's "here document" feature.

[R] SPSS variable lables import

ance, F. Thomas -- .. Dr. Frank Thomas FTR Internet Research 93110 Rosny-sous-Bois France -- __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guid

[R] Drawing a 3-D plot for PCA?

Hi All, Can anyone point me a hint (package) how to draw a 3D plot using the first 3 components from PCA? Thanks a lot, FD [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-h

[R] Bayesian data mining

Hi, I need a help to know whether I can perform the following in R: I have a set of observations (Ns) and each observation is drawn from a poisson distribution with an unkown mean, lambda. The set of lambdas in their turn are drawn from a common prior distribution which is supposed to be a a mix

[R] Bayesian data mining (fwd)

Hi, I need a help to know whether I can perform the following in R: I have a set of observations (Ns) and each observation is drawn from a poisson distribution with an unkown mean, lambda. The set of lambdas in their turn are drawn from a common prior distribution which is supposed to be a a mixtur

[R] Row limit for read.table

2 columns I need. I just wonder why the above 2 points aren't implemented. Maybe they are and I'm totally missing it. Thanks, Frank -- Frank McCown Old Dominion University http://www.cs.odu.edu/~fmccown/ __ R-help@stat.math.ethz.ch ma

Re: [R] Row limit for read.table

er","NULL","NULL","NULL","NULL")) > > dim(f) > [1] 122271 1 > Did you read the help page? Of course I did. For me the definition of colClasses wasn't clear... "A vector of classes to be assumed for the columns&qu

Re: [R] Stacking Vectors/Dataframes

or something like new.frame<-rbind(x.frame,y.frame); # A frame of the right size. new.frame[seq(1,nrow(x.frame),by=2),] <- x.frame # Assign every other row new.frame[seq(2,nrow(x.frame),by=2),] <- y.frame # Assign every other row. -Original Message- From: Laura Quinn [mailto:[EMAIL PRO

RE: [R] Stacking Vectors/Dataframes

Oops there was a bug... new.frame<-rbind(x.frame,y.frame); # A frame of the right size. new.frame[seq(1,nrow(new.frame),by=2),] <- x.frame # Assign every other row new.frame[seq(2,nrow(new.frame),by=2),] <- y.frame # Assign every other row. -Original Message- From: Samuels

[R] precision problems in testing with Intel compilers

performance. If you are not sure whether your application needs this option, try compiling and running your program both with and without it to evaluate the effects on both performance and preci- sion. Has anyone else encountere

RE: [R] More precision problems in testing with Intel compilers

tx[tx <1])#-- Negative exponent -- 7 values [1] 2.002566e-308 2.222849e-308 2.225074e-308 2.225074e-308 2.225074e-308 2.227299e-308 2.447581e-308 Does anyone really care about being correct to 1 unit of machine precision? If you do, you have a bad algorithm. -Original Message---

RE: [R] R-devel and gcc 3.4

I encountered a similar problem with 1. Intel compilers. 2. gcc when I turn off optimizations (no -Ox) -Original Message- From: Robert Kruus [mailto:[EMAIL PROTECTED] Sent: Friday, August 20, 2004 3:14 PM To: [EMAIL PROTECTED] Subject: [R] R-devel and gcc 3.4 I tried compiling

RE: [R] More precision problems in testing with Intel compilers

the < 100*.Machine$double.eps comparisons in nafns.R appear to fail by a factor of about 2.5. There are some examples with extra prints. -Frank -Original Message- From: Martin Maechler [mailto:[EMAIL PROTECTED] Sent: Friday, August 20, 2004 6:39 PM To: Samuelson, Frank* Cc: &#x

[R] Spare some CPU cycles for testing lme?

and devel 2.0.0 (not all possible combinations actually tested.) I've communicated w/ Douglas Bates about this and he doesn't appear to have the problem. Thanks for any help. -Frank __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listin

[R] .Random.seed in R-devel

I'm running R-2.0.0 (yesterday's snapshot)in its own directory, and everything works great, except: >.Random.seed Error: Object ".Random.seed" not found Does 2.0.0 not use .Random.seed for saving, etc, like it says in the help page? Than

Re: [R] .Random.seed in R-devel

Oh. I guess I had a different definition of "when one is required" than the help page. :) Thanks. Prof Brian Ripley wrote: `Like it says on the help page'! Initially, there is no seed; a new one is created from the current time when one is required. Hence, different sessions will

Re: [R] R in Linux Redhat 9.0

rpm -ql R should show you your installed files. Did you try typing "R"? Bryan Tseng wrote: Hi, my name is Bryan and I just installed R on my Redhat 9.0. I imported the key from the net as directed, and tested the integreity of the rpm package- it was fine. I then went ahead and installed the i386.r

[R] Could anyone tell me how to extract pvalue from "lm" fitting?

Dear R people, I have a naive question: after fitting "lm" to a data, I can't extract the pvalue corresponding to a specific covariate in a direct way. Could anyone give me a hint? Thank you very much. Frank __ [EMAIL PROTECTED] mai

Re: [R] Help with Affymetrix data

use PM function in Affy package. Frank On Mon, 4 Oct 2004 16:48:12 -0400, kvyas <[EMAIL PROTECTED]> wrote: > I have CEL files from Affymetrix Mouse Array 430_2 and am trying to get the > the individual PM intensities (11 per gene) for each sample. I would like to > write out

[R] lattice, xyplot layout questions

at didn't work. Is there some other flag to try? Thanks for any help. -Frank __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] pasting indexes to variables within loops

as.numeric(apply(haplo, 1 ,function(x) paste(x,collapse=""))) perhaps? Kristin Kay Nicodemus wrote: Hi all, Hope someone can help me. I start out with a matrix called haplo with 600 rows and two columns. To start with, the elements of the matrix are character strings of 20 numbers. I then want

[R] Could anyone tell me how to plot a hierarchical clustering in a black-white color?

Dear R people, I am trying to do a hierarchical clustering to a matrix using heatmap() function. But I like to have the plot to display in a black-white manner? Could anyone tell me how to do that? Thank you. Frank __ [EMAIL PROTECTED] mailing list

Re: [R] Could anyone tell me how to plot a hierarchical clustering in a black-white color?

Thank you very much, Sundar. That's exactly what I want. Frank On Thu, 11 Nov 2004 15:19:19 -0800, Sundar Dorai-Raj <[EMAIL PROTECTED]> wrote: > > > > > Frank Duan wrote: > > > Dear R people, > > > > I am trying to do a hierarchical clusterin

[R] How to calculate the stratified means in a data frame?

ds of "apply" functions. Could anyone tell me how to do that? Thank you. Frank __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] How to calculate the stratified means in a data frame?

Got it. Many thanks. On Thu, 18 Nov 2004 13:10:53 -0800, Jeff Laake <[EMAIL PROTECTED]> wrote: > look at "by" > > > > Frank Duan wrote: > > > Dear R people, > > > > I have a simple question to ask. Suppose I have a data.frame with two >

[R] Is there a package in R that lets me fit a robust linear mixed model?

Dear R people, Happy Thanksgiving! I just wonder if there is a R package that can supply some kind of "robust" way to fit a linear mixed model. I mean assigning small weights to those observations with large residuals, like iteratively-reweighted-least-squares approach. Many than

[R] Could anyone help me reshape this "wide" data into "longitudinal" one? Thanks

1 2 b 1 3 b 1 4 b 2 1 b 2 2 b 2 3 b 2 4 ... ...... ... ...... ... ...... I always met errors when using "reshape" function.

[R] How to interpret and modify "plot.svm"?

? 2. I want to draw a white-gray colored plot and modify the different colored crosses or circles by different shaped points. Could anyone give me a hint? 3. Is it possible for me to draw a "hyperplane" on the plot? 4. What is the algorithm to plot the contour region? Thank you very mu

Re: [R] How to interpret and modify "plot.svm"?

quot; points into the shapes I want. Basically, I don't want to have any color in the plot. Could you give me a hint how to do that? The code I am using to plot the gray-contour figure is: m <- svm(, data = mydata) plot(m, mydata, col = c(gray(0.2), gray(0.8))) Thank you very much a

[R] help with Tukey Mean-Difference Plot

ver I think I don't get a Tukey Mean-Difference plot. It looks to me that only x=x+y/2 and y=x+y. transformation is performed (a=a+b; b=b+a/2). and the result is displayed as lattice graphic. It would be great if anyone could help me to find the error in my plot. May thanks Frank -- Fr

[R] LDA in R: how to extract full equation, especially constant term

mputing the constant term in its discriminant function. Thanks, -Frank Gibbons PhD, Computational Biologist, Harvard Medical School BCMP/SGM-322, 250 Longwood Ave, Boston MA 02115, USA. Tel: 617-432-3555 Fax: 617-432-3557 http://llama.med.harvard.edu/~fgibbons ___

[R] advise for modeling a linear mixed model

ur measures from the same sample and an effect of the variable "gated" on the frequency. I'm appreciate any help, and apologize in case I have not read carefully the available documentation. Frank -- Frank Mattese-mail: [EMAIL PROTECTED] Department

Re: [R] discriminant function

ust using the '$scaling' That's a start, but there's a long way to go from there -Frank At 09:24 AM 8/26/2003, you wrote: How can I extract the linear discriminant functions resulting from a LDA analysis? The coefficients are listed as a result from the analysis

[R] Cook-distance-type plot (vertical bars)

look in the places I could think of (plot.lm.R, termplot.R, plot.R, aov.R), and couldn't find anything. Maybe I overlooked it? Thanks -Frank Gibbons PhD, Computational Biologist, Harvard Medical School BCMP/SGM-322, 250 Longwood Ave, Boston MA 02115, USA. Tel: 617-432-3555

Re: [R] Cook-distance-type plot (vertical bars)

, this would be something with lots of figures, so that I could browse the figures to find what I want to do, and then look up how to do it. "R for Beginners" goes some way along this path, but perhaps there's something more comprehensive? Thanks again, -Frank PhD, Computation

Re: [R] lda source code

ps. There are certainly others out there who are more au fait with this than me. -Frank Gibbons At 05:47 PM 10/1/2003, you wrote: Wei Geng wrote: I am new to R. Trying to find out how lda() {in MASS R1.8.0 Windows} was implemented in R. Does anyone know where to find out lda source code ? Thanks. He

[R] I have a problem with the log2 function

Dear R users, according the help(log), the function log2(x) should give the natural logarithm of x. I expect in case of x=2 to to get 0.6931, however, R gives me 1 as a result. Similar, logb(2,2) gives 1 again. I'm wondering if I have missed something ? Yours Frank -- Frank Matte

[R] help installing Design package OS-X

kages(ui.pks, CRAN= getOption(where), lib = .libPath()[1] or install.fom.file() npackage installatiion failed in case someone has the binary package for OS-X , please let me know. frank mattes -- Frank Mattes, MDe-mail: [EMAIL PROT

[R] help drawing kaplan-meier plot starting from 0

tive". It would be great, if someone could point me in the right direction where I can find the answer. Many thanks in advance Frank Frank Mattes, MD Department of Virology Royal Free and University Medical School London __ [EMAIL PROTECTED] mailing list http://www.stat.math.ethz.ch/mailman/listinfo/r-help

[R] different plot symbols according to a factor

c I tried plot(a,b, pch=c(1,2)[c]) which didn't work. I appreciate if someone from the list could direct me in the right direction Yours Frank Frank Mattes, MD Department of Virology Royal Free and University Medical School London __ [EMAIL PROT

[R] type III Sum Sq in ANOVA table - Howto?

Hello, as far as I see, R reports type I sums of squares. I'd like to get R to print out type III sums of squares. e.g. I have the following model: vardep~factor1*factor2 to get the type III sum of squares for factor1 I've tried anova(lm(vardep~factor2+factor1:factor2),lm(vardep~factor1*factor2

[R] losing color in panal.xyplot (lattice)

mbols according to the grouping variable. Yours Frank -- Frank Mattes, MDe-mail: [EMAIL PROTECTED] Department of Virology fax 0044(0)207 8302854 Royal Free Hospital and tel 0044(0)207 8302997 University College Med

[R] how to get a line plot before/after treatment

"after". I would like and generate a dotplot with before and after as x-scale, the elements of "before" and "after" on the y-scale and the corresponding elements connected with a line. However, so far I couldn't figure out how to do this in R. If anyone as a

[R] ? building a database with a the great examples

orted in topics, like the help document "Rtips". I have no clue if this is possible to do nor how time consuming the maintaining would be. This is just my view how the help list could be improved Yours Frank -- Frank Mattes, MDe-mail: [EMAIL PROTECTED] Departme

Re: Rtips (was Re: [R] ? building a database with a the great/cookbook

(suggestion , might be a solution http://fawn.unibw-hamburg.de/cgi-bin/Rwiki.pl?RwikiHome. Frank Message: 67 Date: Tue, 3 Jun 2003 00:05:02 -0700 From: Claude Rubinson <[EMAIL PROTECTED]> Subject: Re: Rtips (was Re: [R] ? building a database with a the great I'm wondering if there's dema

Re: [R] avoiding loops

there existing routine that does this? Thanks. -Frank Liaw, Andy wrote: Doesn't this work? d1 <- data.frame(x=rnorm(100), y=rnorm(100), z=rnorm(100)) d2 <- data.frame(x=rnorm(100), y=rnorm(100), z=rnorm(100)) dist.12 <- sqrt(rowSums((d1 - d2)^2)) Andy From: dax42 Hi again, thanks a l

[R] netcdf

month t.unit="month" t.unit=MONTH t.unit="MONTH" all of these result in a return of "unused argument t.unit" I'm apparently not entering the correct value for t.unit, but I can't find a list of options anywhere Any help? Thanks Frank Aebly _

[R] How to create a vector with "one", "two", "three", ...?

Hi R people, I met a naive prolem. Could anyone give me a hint how to create such a vector with entries: "one", "two", "three", ...? I must have seen this problem somewhere else but I can't find that source now. Sorry to bother you with such a s

Re: [R] How to create a vector with "one", "two", "three", ...?

Sorry, I didn't get the question clear. What I meant is to create a character vector with length 200: "one", "two", "three", ..., "two hundred" On 4/15/05, Federico Calboli <[EMAIL PROTECTED]> wrote: > On Fri, 2005-04-15 at 14:30 -0400, Fr

Re: [R] How to create a vector with "one", "two", "three", ...?

Dear John, That's exactly what I want. Millions of thanks, Frank On 4/16/05, John Fox <[EMAIL PROTECTED]> wrote: > Dear Frank, > > This was an interesting exercise. Here's a solution: > > numbers2words <- function(x){ >helper <- function(x){ >

[R] looping over factors

factors? Thanks for any help. -Frank __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] moments, skewness, kurtosis

I checked the help and the mailing list archives, but I can find no mention of a routine that calculates higher moments like skewness and kurtosis. Of course, these are easy enough to write myself, but I was thinking that they MUST be in here. Am I wrong? Thanks. -Frank

  1   2   3   4   5   6   7   8   >