Re: [R] Problem loading package 'JGR' using R-2.15.0 (Win32).

guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Robert W. Baer, Ph.D. Professor of Physiology Kirksville College of Osteopathic Medicine A. T. Still University of Health Sciences 800

Re: [R] Histogram from a table in R

oject.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Robert W. Baer, Ph.D. Professor of Physiology Kirksville College of Osteopathic Medicine A. T. Still University of Health Scien

[R] Applying a function to categorized data?

ck(v1, v1$ARI_MIT ~ v1$MONTH)) { write(mean(i), "") } This works, but how do I get the data into a data frame? With the month labels in a column? They are not avaliable inside the loop body. I know I need to get a book on R. Thanks, robert __

Re: [R] Applying a function to categorized data?

ow (record) oriented thinking and get my head into vectors and columns. I'm still nowhere near where I think I need to be in order to work mit my data. I'll get back to the list when I have pinpointed my problem a bit better, and I'll also supply some sample data. Have a nice weekend

[R] How to create a data.frame from several time series?

frame, but I don't know how to do that because I don't know beforehand how many series I create in the for() loop, how to distinguish them by (unknown beforehand) tool names, and how to supply them to ts.union. Thanks, robert CODE HERE library(zoo)

[R] Finding percentile of a value from an empirical distribution

a number z and vector x (i.e. the raw vector).It returns percentile of z wrt x. E.g. > f(71,x) Should return something around 0.708 or 0.709. I am wondering if there is any pre-packaged functions that do this in R? If not, how can I write such a function. Any suggestion would be ap

Re: [R] Finding percentile of a value from an empirical distribution

Thank you. That's easier than I thought. On Wed, Jan 11, 2012 at 12:06 AM, Jorge I Velez wrote: > Hi Robert, > > Try > > set.seed(123) > x <- seq(100) > x <- sample(x, 1000, replace = TRUE) > f <- ecdf(x) > f(10) > # [1] 0.099 > f(71) > # [1]

Re: [R] help to slip a file name using "strsplit" function

t;), "_", fixed = TRUE) Thanks in advance Gianni [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] x-axis label for boxplot

xlab = needs to be names = see ?boxplot boxplot(jdata[,"x"],jdata[,"y"],jdata[,"z"],names=c("x values","yvalues","zvalues")) -Original Message- From: John Sorkin Sent: Wednesday, February 08, 2012 6:48 PM Cc: r-help Subject: [R] x-axis label for boxplot I am trying to produce th

[R] Find exit status inside on.exit

these function with a tryCatch. Something like this: isError <- function() { parentError <- ??? if (!parentError) cat("No errors here!\n") } X <- function() { on.exit(isError()) 1+1 } > X() [1] 2 No errors here! Thanks, Robert Robert

[R] How would you program an Adverse Events statistical table using R code?

rent statistical procedures on the same page. Robert Wilkins __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] Cleaning up messy Excel data

s end, the choice might be context sensitive. Rob ---- -- Robert W. Baer, Ph.D. Professor of Physiology Kirksville College of Osteopathic Medicine A. T. Still University of Health Sciences 800 W. Jefferson St

Re: [R] ggplot2 dot chart-start at zero

___ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Robert W. Baer, Ph.D. P

Re: [R] How to create a data.frame from several time series?

ding how it works and why. Regards, robert __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] r-package RDCOMClient

These sorts of issues are often machine dependent and not R package or R software dependent, especially when it comes to processes like com. I a,m not a RDCOMClient user, but it's been on my list to look at. You don't provide much info on your various component versions (see posting guide) so

Re: [R] [R-sig-hpc] Quickest way to make a large "empty" file on disk?

machine: > ptm <- proc.time() > zz <- file("testbin.bin", "wb") > for(i in 10) writeBin(rep(0,1),zz, size=16) > close(zz) > proc.time() - ptm user system elapsed 2.416 1.728 16.705 Otherwise I would suggest writing a little piece of c cod

[R] Replacing tick labels in a plot

effort into actually giving it a shot. Any ideas? BTW, is it normal that each and every post to this list, even by list members, has to go through moderator approval? Thanks, robert __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listi

Re: [R] Replacing tick labels in a plot

sed when you signed up? > That's a frequent cause of moderation. > > Sarah > > On Fri, May 4, 2012 at 1:09 PM, Robert Latest wrote: >> Hello, >> >> is it possible to replace the text of tick marks in a plot? >> Specifically, I'd like to have a ppnor

Re: [R] How to read ANOVA output

with plenty of R code and examples. The only downside to the latter text is that it has numerous typos that seem to have escaped the editing process. If you haven't already, I'd check these two resources. Rob -- Robert W. Baer, Ph.D. Professor of

[R] How to deal with a dataframe within a dataframe?

hanks, robert > agg <- aggregate(formula=df$value ~ df$quarter + df$tool, + FUN=cp.cpk, lsl=1300, usl=1500) > head(agg) df$quarter df$tool df$value 1 09Q3VS1A 1.800534, 1.628483 2 10Q1VS1A 1.299652, 1.261302 3 10Q2VS1A 1.699018, 1.381570 4

Re: [R] How to deal with a dataframe within a dataframe?

es = c("df$quarter", "df$tool", "df$value"), row.names = c(NA, 6L), class = "data.frame") > I would like this in either the form of a "flat" data frame (i.e., the contents of "df$value" as two separate columns), or -- even preferable --

[R] How to apply functions across columns?

417 1393 1346 1328 I'd like to add a column to this that gives, for each row, the averages of the values in the columns s1 to s5. Really primitive. But I totally absolute don't understand how to do this. I don't need any "intelligence", I know my values are always in columns

Re: [R] How to apply functions across columns?

I need: tencor$mean <- apply(tencor[4:8], 1, FUN=mean) This way I'm also not tied to just mean() as aggregator but can use any homemade function (this would have been my followup question had I followed your advice ;-) Thanks! robert __ R-help@r-p

[R] ANOVA question

"on foot" by looping over the levels of my single factor, but I suppose this gets complicated and messy with more complex models. Any ideas? Thanks, robert set.seed(0) N <- 100 # sample size MEAN <- c(10, 20, 30, 40, 50) VAR <- c(20,20,1, 20, 20) LABELS <- c("A", &q

Re: [R] ANOVA question

E 49.78282 0.3937038 126.44738 0.00e+00 Thanks, robert __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-co

[R] How to interpret an ANOVA result?

out six explanatory variables, some continuous, some factorial). Already in the DOE phase I want to understand what I will be doing with the data later (this is for a Six Sigma project in an industrial production environment, in case anybody wants to know). Thanks, robert Here's the full d

Re: [R] How to interpret an ANOVA result?

> > Anyway, I thought your data set was interesting, so I took the approach that > comes to my mind.  Here it is.  It might be pretty much self-explanatory, if > not, try ?aov and ?TukeyHSD for details. TukeyHSD looks interesting. I'll look into it. Thanks, robert _

[R] How to include known errors in a regression?

d pitch all of this together an a linear model. Makes sense? My goal is to replace one (slow, expensive) measurement by another (fast, cheap) one, and I need to establish the correlation (and especially the expected error margin) between the two to see if it is feasible. Thanks, robert _

Re: [R] Regression Analysis or Anova?

orm. Good luck, robert On Tue, May 15, 2012 at 10:49 PM, Andrea Sica wrote: > Dear all, > > I hope to be the clearest I can. > Let's say I have a dataset with 10 variables, where 4 of them represent for > me a certain phenomenon that I call Y. > The other 6 represent for me anot

Re: [R] how to find outliers from the list of values

ead the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Robert W. Baer, Ph.D. Professor of Physiology Kirksville College of Osteopathic Medicine A. T. Still University of Health

Re: [R] How to create axis y axis for horizontal bar plot?

with scale 0 to 12 but i need y axis. How can i implement it? http://r.789695.n4.nabble.com/file/n4630478/barplot2.jpg Thanks Try: gears = c(3,4,5) b = barplot(gears, horiz=TRUE) axis(2, at=b, labels=c('three','four','five')) Rob -----

Re: [R] Names of Greek letters stored as character strings; plotmath.

eed each time by frustrating trial and error. Can I ask how you guys got a handle on the bigger (besides your usual brilliance )? Is there more comprehensive documentation in the developer literature or is there a user wiki that you would recommend for those who never quite get the big picture? If not, this would be a worthy t

Re: [R] how to check given number seq. is time series or not?

://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. -- Robert W. Baer, Ph.D. Professor of Physiology Kirksville College of Osteopathic

[R] inconsistent output when using variable substitution

nd correct behavior for R, or am I missing something? And if it is normal and correct, what do I need to do to get the behavior I am expecting. Thanks in advance for any advice. Rob -- Robert Felty http://robfelty.com "In this house we obey the laws of thermodynamics!" --

Re: [R] inconsistent output when using variable substitution

j[1] [1] "49" > thisSubj = subset(Trials,Trials$Subj=="49") > thisSubj$Ansr[1] [1] "able" > thisSubj = subset(Trials,Trials$Subj==Subj[1]) > thisSubj$Ansr[1] [1] "abacus" The problem is that I want to use this in a loop, so that I can get

Re: [R] inconsistent output when using variable substitution

e final comma, and had come to the conclusion that this method only worked when selecting a particular column of a data frame. Now I know better. Thanks so much for everyone's help. I am very impressed by the R community Rob -- Robert Felty http://robf

Re: [R] tolstoi.newcastle dead?

Hello useRs, tolstoy is back up and running again. Robert On Sunday 01 June 2008 8:07, Carl Witthoft wrote: > Anyone know what's up w/ tolstoi.newcast.edu.au (home of the threaded > R-help archives)? I've got a "can't connect" for several days now. > > Carl

[R] read.spss {foreign} doesn't work over network?

n" "C:/PROGRA~1/R/R-27~1.0/library" "0.1" multcomp "multcomp" "C:/PROGRA~1/R/R-27~1.0/library" "1.0-0" mvtnorm "mvtnorm" "C:/PROGRA~1/R/R-27~1.0/library" "0.9-0" nlme "nlme&

Re: [R] read.spss {foreign} doesn't work over network?

Sorry, the error was "OTOSOTK" {On The Other Side Of The Keyboard}. Someone renamed one of the subdirectories (60 to 06) and I couldn't see it. Sorry for the bandwidth Robert Farley Metro www.Metro.net -Original Message- From: Prof Brian Ripley [mailto:[EMAIL P

[R] Survey questions

entry level discussion of the statistical ramifications of this? Is it possible (proper) to re-survey a population to obtain new weights (for previously un-weighted responses) for an original survey? Thanks in advance, Robert Farley Metro 1 Gateway Plaza Mail Stop 99

[R] kriging problem(?)

ng kriging. Just to avoid misunderstandings: I do not just want to get a spatial grid of Sa and Sb, but I want to find on which Scale Sb affects Sa. cheers, Robert, Robert Ptacnik Norwegian Institute for Water Research (NIVA) Gaustadalléen 21 NO-0349 Oslo FON +47 982 2

Re: [R] Survey questions

the details of designing a "post survey survey" to generate additional weights? Thanks in advance, Robert Farley Metro www.Metro.net ________ From: Farley, Robert Sent: Thursday, June 26, 2008 15:48 To: 'r-help@r-project.org' Subje

Re: [R] A question about using function plot

How about: plot.default(factor(x), y, xaxt='n', xlab='x') axis(1, at=factor(x), labels=x) - Original Message - From: "Jason Liao" <[EMAIL PROTECTED]> To: Sent: Monday, July 14, 2008 10:00 AM Subject: [R] A question about using function plot Hello, everyone! I have spent two hou

[R] r^2 for SSasymp?

Dear R Help, At the risk of sounding naïve, is there a way for obtaining something resembling an r^2 value for the fit obtained using SSasymp? I have read ?SSasymp and looked through the archives but to no avail. Many thanks, Robert Terwilliger __ R

[R] ARIMA and sig. tests

Why in arima output there is no significance tests of estimated parameters as in regression models? I know that I have s.e. and I could build these test by myself but are there any reasonable reasons? Best, Rob __ R-help@r-project.org mailing list ht

[R] Nonlinear equation

Dear R Users, I'm trying to find function that allow me to solve one nonlinear equation. Functions that I found are good for optimization problems. Any suggestions are welcome, rob __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/list

[R] Fwd: Re: Nonlinear equation

> > >From: [EMAIL PROTECTED] > > >Date: 2008/04/02 Wed AM 09:58:28 CDT > > >To: r-help@r-project.org > > >Subject: [R] Nonlinear equation > > > > hi: you need to give an example and details or > > you won't get much response, if any. Equation e.g. (A, B are known constants): 3log(gamma(x))-log(ga

Re: [R] Fwd: Re: Nonlinear equation

> On 4/2/2008 11:24 AM, [EMAIL PROTECTED] wrote: > >> > >From: [EMAIL PROTECTED] > >> > >Date: 2008/04/02 Wed AM 09:58:28 CDT > >> > >To: r-help@r-project.org > >> > >Subject: [R] Nonlinear equation > >> > > >> > hi: you need to give an example and details or > >> > you won't get much response, if

[R] newbie subset question

03 1 62 48 42 227 school 00 0 00 0 00 shop 00 0 00 0 00 soc/rec00 0 00 0 00 university 00 0 00 0

[R] SPSS to R Data file conversion

I wrote a little routine to convert multiple spss data files (as data frames) to R data files. The code is as follows: # list=dir(pattern=".sav") library(foreign) for (i in 1:length(list)){ # The saved data frame will be dat dat=read.spss(list[i],to.data.frame=TRUE) name=substring(list[i],1,nchar

Re: [R] Histogram Label Font Size

Here's an example to make the text 18 point. You can set ps=6 to make small text. x=rnorm(100) opar=par(ps=18) # Make text 18 point hist(x) opar Rob Baer - Original Message - From: "Sue Lee" <[EMAIL PROTECTED]> To: Sent: Monday, April 14, 2008 2:35 PM Subject: [R] Histogram Label Fon

Re: [R] [R-sig-ME] Post hoc tests with lme

At the risk of annoying Doug Bates, I'll point out that "glht" in the multcomp package works with lmer objects. In fact, you can supply your own degrees of freedom value via the "df" argument (a fact which is not immediately obvious in the glht help page). If you don't supply a df value, it will

Re: [R] Finding a path using the Graph package

c cba cb bc. > > Thanks in advance, > > Dan > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal,

Re: [R] Mantel-Haenszel for 2x2

See if this is what you are looking for: library(epicalc) data(Oswego) use(Oswego) cc(ill, chocolate) mhor(ill, chocolate, sex) You might also wish to read the details of: mantelhaen.test HTH], Rob Baer - Original Message - From: "Dr. Jeff Miller" <[EMAIL PROTECTED]> To: <[EMAIL PROTEC

[R] logit GLM without intercept

, binomial) Does this code yield the correct model or is there another possibility? I’d appreciate it very much if you could help me out with this. I attached some example data. Thanks & all the best Robert ---- Robert Junker Department of Animal Eco

Re: [R] Documentation General Comments

>> I realize the R developers are probably overwhelmed and have little time >> for this, but the documentation really needs some serious reorganizaton. >> A good through description of basic variable types would help a lot, >> e.g. the difference between lists, arrays, matrices and frames. > > Agre

[R] time zone conversion

2007-04-01 12:00:00 EDT = 2007-03-01 17:00:00 BST Is any of the time zone/daylight savings information coded up somewhere? Thanks, Robert Robert McGehee, CFA Geode Capital Management, LLC One Post Office Square, 28th Floor | Boston, MA | 02109 Tel: 617/392-8396Fax:617/476-6389 mailto:[EMAIL

Re: [R] time zone conversion

43537 language R version.string R version 2.6.1 (2007-11-26) > -Original Message- From: McGehee, Robert Sent: Monday, April 28, 2008 4:27 PM To: r-help Subject: time zone conversion Hello, I'm trying to convert times in the E

[R] A Maze Generator

I had some fun this afternoon coding up a 'maze generator' in R. I thought I'd pass along the fruits of my labor for everyone's amusement. As written, every point is connected to every other point, so feel free to 'start' and 'finish' anywhere you like. Ha

Re: [R] A Maze Generator

main="Simple Maze") plotMaze(makeMaze(25), main="Hard Maze", lwd=1) -Original Message- From: Peter Dalgaard [mailto:[EMAIL PROTECTED] Sent: Tuesday, April 29, 2008 5:34 PM To: McGehee, Robert Cc: r-help Subject: Re: [R] A Maze Generator McGehee, Robert wrote: > I ha

Re: [R] Gravity (spatial interaction) models in R

#x27;m only a newbie, but from what I've seen on the list, this should be simple to do As to calibrating F-Factor functions: one often "chooses" a functional F-Factor form, then "adjusts" the coefficients, until the weighted trip length distribution of the resulting tri

Re: [R] Monty Hall simulation

then you can develop a winning strategy. If you were wrong about either the Host or the prize mix, that strategy may be harmful. Robert Farley Metro www.Metro.net -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Johannes Huesing Sent: Sunday, M

[R] Help creating a correlation matrix

= "complete.obs")} This is wrong, of course, but I don't know how to "tell" R to run through D2_1, D2_2,...,D2_14 giving me the correlations for each pair. Any help or hints would be appreciated. Sincerely, Robert O'Brien __ R-help@r-

Re: [R] Rgraphviz

e posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Robert Gentleman, PhD Program in Computational Biology Division of Public Health Sciences Fred Hutchinson Cancer Research Center 1100 Fairview Ave. N, M2-B87

[R] fit linear regression with multiple predictor and constrained intercept

intercept of each factor level. Robert -- NIVAs hovedkontor har flyttet til nye lokaler i CIENS - Forskningssenter for miljø og samfunn; Gaustadalléen 21, 0349 Oslo. Meld deg på vå

[R] relative importance of predictors

Regression Models). Now I want to compare the relative importance of my predictors. Can I rely on the output of 'relaimpo' here? (which takes 'lm', not 'rlm&#

[R] The R Book - great resource for R beginners

Hi all, I've recently discovered "The R Book" by Micheal J. Crawley. I am new to R and I am finding it incredibly useful. With my background in math and programming, I expected that I would be able to pick it up quickly, but even after reading several other intro books, I found it was still ha

Re: [R] The R Book - great resource for R beginners

It was not my intention to start any kind of war. I do want to thank Katharine Mullen for sharing to link to the review in R News. It helps me understand why there might be some controversy over my positive experience with The R Book. As I mentioned in an earlier post, I'm working my through th

Re: [R] Analyzing Publications from Pubmed via XML

InternalNodes = TRUE) > sapply(c("//author", "//category"), xpathApply, doc = doc, fun = xmlValue) > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide

[R] termplot reference

Dear Group, Is there a specific reference for how termplot partitionates residuals (part=T)? I want to use figures from termplot in a publication and wonder about an appropriate reference (nothing in the help-file) Thanks Robert

[R] OS X binary: 32 or 64-bit?

3321675776) failed (error code=3) R(2081,0xa000d000) malloc: *** error: can't allocate region R(2081,0xa000d000) malloc: *** set a breakpoint in szone_error to debug My Mac is 64-bit, however, I'm wondering if R is compiled as a 32-bit app? Many

[R] binning data

I would like to 'bin' a vector of of values (>10 values) in much the same way as 'hist' does for plotting, but I just want a result with the bins and the frequencies. Is there a 'quick and dirty' (built-in) way to do this, or do I need to write some code?

[R] match and incomparables

ed? Note from ?match: "incomparables: a vector of values that cannot be matched. Any value in x matching a value in this vector is assigned the nomatch value. For historical reasons, FALSE is equivalent to NULL." Thanks in advance! Robert platform i386-pc-mingw32

Re: [R] How to draw the graph of f(x,y) = x * y ?

You could always try F(x,y) = f(x) + 0*y That is "zero out" the degenerate dimensions. Of course you'll be plotting what is essentially a two dimensional object as if it were three dimensional. The degeneracy in y means a 2-D curve will be "extruded" along the Y di

[R] (nlme) Repeated measures with continuous covariate in lme

has the same value within each year. On the other hand it has different values within each plot. I hope someone can point me in the right direction. Thanks in advance, Robert Buitenwerf South African Environmental Observation Network Phalaborwa South Africa _

[R] GPL and adds-on

I'm wondering is it in accordance with the law, taking into account GNU GPL on which R is based, that SPSS have adds-on module that allow you to do a analysis in R. I mean commercial software use R. I guess that nobody can create R-GUI (e.g. in JAVA) on commercial rules (licensing, pricing, sell

[R] newbie Syntax questions: Re-Calculating expansion factors

1 1 0 0 2 9013240 0 0 1 0 3 0 0 20 0 0 0 0 ... and so on Robert Farley Metro 1 Gateway Plaza Mail Stop 99-23-7 Los Angeles, CA

[R] pcnm - scaling

to using pcnm Thanks ! Robert Robert Ptacnik Norwegian Institute for Water Research (NIVA) Gaustadalléen 21 NO-0349 Oslo FON +47 982 277 81 FAX +47 221 852 00 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do re

[R] Ashlee Vance's article on R in the New York Times

Ashlee Vance's article on R in the New York Times. This is typical of the New York Times. Because they get to coast on the prestige and reputation of their brand , they have a history of just this sort of journalistic sloppiness. Whether it's the author or the editor at fault doesn't really matter

[R] The AT&T researchers and the New York Times

Is anyone in the leadership of the R-project going to contact the New York Times and clarify that the article gave remarkably short shrift to the people who designed the user interface for R, to a large extent AT&T researchers from an earlier generation? It would be the appropriate thing to do. Th

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

Is anyone in the leadership of the R-project going to contact the New York Times and clarify that the article gave remarkably short shrift to the people who designed the user interface for R, to a large extent AT&T researchers from an earlier generation? It would be the appropriate thing to do. Th

[R] survey statistics, rate/proportions with standard errors

what does R have to compare with , say , proc surveymeans, estimate survey means/proportions with standard errors, using Taylor methods? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listi

[R] What does R have for age-adjusted survey analysis?

A procedure that , after adjusting for sampling weights, also explicitly does an age adjustment to conform with an age distribution of an older census? __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the p

[R] Extraction from an output

08 How to use ca.test and extract p-value and chi-squared. Robert __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

[R] Useful books for learning the R software and the S programming language

any useful books for learning the R statistical software? are they pricey? and if the books recommended focus on S, how compatible will they be for someone learning R? thank you in advance for your help. P.S. specialized survey statistical procedures? Is R good at that? [[alternative H

[R] Significance levels of variance terms

I would like to be able to extract the significance levels (or standard errors) of the variance-covariance parameters in nlme objects (i.e. from the varClasses and/or corClasses). Is there a way to do it, or do I have to reverse the intervals function calculations ? Thanks. Rob

[R] Gentleman and Ihaka's integrity in question

It does look like Gentleman and Ihaka not only lied to the New York Times, but also to the New Zealand Herald and who knows who else. This is disgusting. The R programming language is the S programming language, and Gentleman and Ihaka are not the ones who designed it. http://thenewyorktimesisslop

[R] Comparing weighted histograms?

5 1.15 1.15 1.15 1.15 1.15 1.15 1.15 1.15 1.15 1.15 1.15 1.15 1.15 1.15 mode = >1 mode Valid n = 87 87 categories - only first 18 shown > range(Survey

[R] Problems with svyhist

der=TRUE, sep=",", na.strings="NA", dec=".", strip.white=TRUE) MyBreaks <- c(0,5,10,15,20,30,40,50,60,80,100,120) hist(SurveyData$Miles2LAX, breaks = MyBreaks, main = paste("Histogram of unweighted distances" ) ) ##

[R] How do I get xtabs to show missing data?

8576 0.00 ..Refused 0.00 0.00 > Robert Farley Metro 1 Gateway Plaza Mail Stop 99-23-7 Los Angeles, CA 90012-2952 Voice: (213)922-2532 Fax:(213)922-2868 www.Metro.net [[alternative HTML version deleted]] __

[R] R Style Guides?

rably with examples. Suggestions, please? Thanks Robert Biddle __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self

Re: [R] KM estimation for interval censoring?

HBERGER, M. Survival Analysis. New York: > Springer Verlag, 1997, for the algorithm. > > This MIT-CSAIL report says they have an available R package that handles > interval censoring and time-dependent covariates: > > http://dspace.mit.edu/bitstream/1721.1/33957/1/MIT-CSAIL-TR-2

[R] Cumulative multinomial regression using VGAM

es either. Am I missing something ? Thanks. Robert _ [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLE

[R] j and jcross queries

liklihood that Y and Z come from different processes? I would appreciate advice. Cheers Robert Biddle __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

[R] compiling 2.6.2 using icc

ow shows that it is using gcc rather than icc when compiling the package, and the CFLAGS wd188 and ip . How do I get it to use icc for both building R and building the packages? Robert .. make[2]: Entering directory `/opt/freewaresrc/R/R-2.6.2/src/library/Recommended' begin installing rec

Re: [R] compiling 2.6.2 using icc

-lipgo -lirc -lgcc_s -lirc_s -ldl FCPICFLAGS = -fpic FPICFLAGS = -fpic The only reference to gcc I could find was OBJC = gcc Could this be a problem? I tried changing this to icc -c99 but this made no difference. Since it seems to work for other users, I assume there is something in my s

Re: [R] compiling 2.6.2 using icc

Just in case others make my mistake, the problem was that I had a .R directory in my home directory, which had a Makevars file in it with CC=gcc. I don't remember why I had this, but it was being read for package building and throwing everything out. Robert -Original Message-

[R] using subset() in data frame

ect=macro)) None of th above, plus many permutations thereof, works. Can anyone educate me? Thanks, Robert Walters __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] using subset() in data frame

Erin Hodgess wrote: > Hi Robert! > > Could you please check > str(data.b) > > and see what you have for porosity? It needs to be a factor > > Thanks, > Erin > > Erin, Yes, porosity is a factor. See output below. > str(data.b) 'data.frame': 9

Re: [R] using subset() in data frame

Chuck Cleland wrote: > On 2/23/2008 6:09 AM, Chuck Cleland wrote: >> On 2/22/2008 8:01 PM, Robert Walters wrote: Chuck, Thanks for the pointers on subset(). When I submit the two variants you suggested, below: > fit1 <- lm(y ~ x, subset(data.b, porosity == "macro"))

<    1   2   3   4   5   6   7   8   9   >