[R] Incomplete time series

te number of rows starting at the appropriate number of days since epoch, and then using match() to combine columns and add the desired rows. I'm new to time series in R, but it seems like there should be an easier way. I've gotten as far as Google and rseek can take me so any help would b

[R] Triangular distribution for kernel regression

options that I can't seem to use to produce estimated values (e.g. density). Any help is appreciated. Bryan [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do rea

Re: [R] Triangular distribution for kernel regression

ittle more direction. I've understand that it is connected to local polynomial regression but I can't seem to have any success from that direction either. At this point the only package that is giving smoothed estimates as I would expect is ksmooth - which doesn't include the appropriate d

Re: [R] Triangular distribution for kernel regression

atson estimator, which is where I need to be, but I'm not sure how to manipulate the bandwidth, as would be possible in other methods. Can you clarify this at all? Bryan On Sat, Sep 12, 2009 at 3:46 PM, Gabor Grothendieck wrote: > What about kknn -- that was listed as having t

Re: [R] Triangular distribution for kernel regression

t; You can specify the number of nearest neighbors. > > On Sat, Sep 12, 2009 at 5:16 PM, Bryan wrote: > > I originally looked over kknn because I need to be able to specify a > > bandwidth parameter. I am trying to replicate some previous non-R work > in > > R, so I can

Re: [R] predict and arima

Indeed, there was a bug ... my current play code looks like this ... get.best.arima <- function(x.ts, maxord=c(3,3,3,3)) { # function based on 'Introductory Time Series with R' # ... try and fit the best ARIMA(p,d,q,P,D,Q) model # using all permutations from 0 to max

[R] Rejection metho..d...

Generate random variables which flow f(x)=20*x*(1-x)^3 , 0http://r.789695.n4.nabble.com/Rejection-metho-d-tp4208138p4208138.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailm

Re: [R] Inefficiency of SAS Programming

of other contexts in which their pricing structure would work. Bryan On Fri, Feb 27, 2009 at 12:48 PM, Frank E Harrell Jr < f.harr...@vanderbilt.edu> wrote: > spam me wrote: > >> I've actually used AHRQ's software to create Inpatient Quality Indicator >> reports. I can

Re: [R] Bootstrapping in R

each time the code is run, correct? result <- boot(n_data, statistic = DataSummary, R = 100). Best, Bryan Mac bryanmac...@gmail.com > On Sep 29, 2016, at 12:16 PM, ruipbarra...@sapo.pt wrote: > > Hello, > > Read the help page ?boot::boot. > For instance, try the follo

Re: [R] Bootstrapping in R

orm(100) > y <- x + rnorm(100) > dat <- data.frame(x, y) > stat2 <- function(DF, f){ > model <- lm(y ~ x, data = DF[f,]) > coef(model) > } > > boot(dat, stat1, R = 100) > boot(dat, stat2, R = 100) Bryan Mac bryanmac...@gmail.com > On Oct 2, 20

[R] Hierarchical Cluster Analysis

Bryan Mac Data Scientist Research Analytics Ipsos Insight LLC [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] printing a data.frame that contains a list-column of S4 objects

Is there a general problem with printing a data.frame when it has a list-column of S4 objects? Or am I just unlucky in my life choices? I ran across this with objects from the git2r package but maintainer Stefan Widgren points out this example below from Matrix as well. I note that the offending o

[R] Least Median Square Regression

squared_lms_sqrtnar_sqrtnic) return(out) } Also, which value should be looked at decide whether this is best regression model to use? Bryan Mac bryanmac...@gmail.com [[alternative HTML version deleted]] __ R-help@r-project.org mailing l

Re: [R] Least Median Square Regression

I am confused reading the document. I have installed and added the package (MASS). What is the function for LMS Regression? Bryan Mac bryanmac...@gmail.com > On Oct 8, 2016, at 6:17 AM, Enrico Schumann wrote: > > On Sat, 08 Oct 2016, Bryan Mac writes: > >> Hi R-help,

Re: [R] Bootstrapping in R

Hi, After running the bootstrapping, I would like to the output of the bootstrapped samples. How can I view the bootstrapped samples of each variable? Bryan Mac bryanmac...@gmail.com > On Oct 18, 2016, at 3:57 AM, Rui Barradas wrote: > > It means that the sd of the bootstrap s

[R] Filtering out rows

Hi, I am looking to filter out rows in my data set, run my analysis and then removing the filter. Here is my regression. I would like to filter out a row. For example, I would like to filter out a row, "Case"=1511. npi=lm(npi_mvmt~cnavgpi,data=df) summary(npi) resid.vs.fitt.LS%https://stat.ethz.c

[R] Export R output in Excel

Hi, How do I export results from R to Excel in a format-friendly way? For example, when I copy and paste my results into excel, the formatting is messed up. Thanks. Bryan Mac bryanmac...@gmail.com __ R-help@r-project.org mailing list -- To

Re: [R] [FORGED] Export R output in Excel

not find most of the hints “useful”. if anything, it got me more confused. Thanks. Bryan Mac bryanmac...@gmail.com > On Dec 28, 2016, at 3:15 PM, Rolf Turner wrote: > > On 29/12/16 10:45, Bryan Mac wrote: >> Hi, >> >> How do I export results from R to Excel in

Re: [R] Export R output in Excel

a .CSV file with the good formatting ? I am looking to export the whole output if possible. I found this code, but it doesn’t cover the whole output of the console. write.csv(coef(summary(test)), file=“test.csv”) My whole output consists of descriptives and regressions. Best, Bryan Mac

[R] Conditionally Updating Lattice Plots

, from which I would infer that you can't update the update or I'm not calling it correctly. I have a nagging sense too that the "real" way to do this is with a non-standard use of panel.superpose but I don't quite see how to do that from available examples. TIF f

[R] Lattice Version of grconvertX or variant on panel.text?

to be what I need. Is there a function like grconvertX in Lattice, or is there a flag or some other method of making panel.text use plot coordinates? Thanks, Bryan __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] Lattice Version of grconvertX or variant on panel.text?

Never mind, I just hard-coded it using ratios. Simpler than I thought. Thanks, Bryan On 7/20/08 9:03 PM, "Bryan Hanson" <[EMAIL PROTECTED]> wrote: > Still playing with Lattice... > > I want to use panel.text(x, y etc) but with x and y in plot coordinates > (0,1),

[R] Properly Parsing Pre-Superscripts & Displaying Them With grid.text

n(^37*Cl), ": 25%", sep = "") # This doesnt' produce an error, but doesn't produce what is wanted either, # as the expression is taken (almost) literally: leg.txt3 <- paste("Based upon:\n", expression(""^35*Cl), ": 75%\n", expression(&quo

Re: [R] Properly Parsing Pre-Superscripts & Displaying Them With grid.text

, you have to quote it in the expression: "%" which I suppose is a general feature. I may have missed it, but that behavior doesn't seem to be mentioned in the plotmath help page - perhaps it's too obvious? Thanks, Bryan On 8/2/08 2:19 PM, "Gavin Simpson" <[EMAIL

Re: [R] Are there any guis out there, which will allow editing of the graph?

this was possible on some R code that went with a Wikipedia entry, and he tried it, and it works. YMMV. Bryan * Bryan Hanson Professor of Chemistry & Biochemistry DePauw University, Greencastle IN USA On 8/4/08 1:09 PM, "Bert Gunter" <[EMAIL PROTECTED]> wr

[R] Need Advice: Considering Converting a Package from S3 to S4

t; as they say. If methods are not critical, does it make sense to spend the time making the change? Any perspective and advice would be welcomed. Thanks in advance, Bryan * Bryan Hanson Professor of Chemistry & Biochemistry DePauw University, Greencastle IN USA

[R] Selecting/Accessing the last vector in a list of a list of data.frames

simple approach, but I may be too tired right now to see it! Thanks, Bryan * Bryan Hanson Professor of Chemistry & Biochemistry DePauw University, Greencastle IN USA __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/l

Re: [R] Selecting/Accessing the last vector in a list of a list of data.frames

think it should do, much less what it does do! I think an easier approach for me will be to re-write the function that generates "test" so it is simpler to extract what I need. I will think on it. Thanks, Bryan > temp <- lapply(test, lapply, '[', 'V1') >

Re: [R] Lattice in a loop does not produce output

Lattice objects must be assigned and deliberately printed: > png("test.png") > p <- xyplot(y~x|z) > plot(p) > dev.off() Should fix both problems. Bryan * Bryan Hanson Professor of Chemistry & Biochemistry DePauw University, Greencastle IN USA On 8

Re: [R] Google's R Style Guide (has become S3 vs S4, in part)

d made Myobject an S4 object, then I would have to go back, redefine the object, update validObject, and possibly write some new accessor and definitely constructor functions. At least, that's how I understand the way one uses S4 classes. Back to trying to get something done! Bryan **

[R] Recursion is slow

lt <- length(T) if (lt == 1) { if (0 <= m & m <= T[1]) { return(1) } else { return(0) } } R <- 0 for (u in 0:T[lt]) { R <- (R+(A(T[1:(lt-1)],(m-u } return(R) } - Bryan Keller, Doctoral S

Re: [R] Recursion is slow

. Bryan - Bryan Keller, Doctoral Student/Project Assistant Educational Psychology - Quantitative Methods The University of Wisconsin - Madison - Original Message - From: William Dunlap Date: Thursday, September 3, 2009 6:41 pm Subject: RE: [R] Recursion is slow To: Bryan Keller

[R] Matrix as input to xyplot {lattice} - proper extended formula syntax

his, or at least I don¹t recognize it if I have! Is there a more elegant way to tell xyplot I want to use each row of y repeatedly with the same x, in a loop-like fashion? TIA. Bryan ***** Bryan Hanson Professor of Chemistry & Biochemistry DePauw University, Greencastle IN USA

Re: [R] Matrix as input to xyplot {lattice} - proper extended formula syntax

of y columns, from y.1 to y.n let's say. Is there an easy way of creating the phrase y.1+y.2+...y.n to pass to xyplot, or even better, some sort of syntax that says "take all y.n" and plot them against x? Thanks, Bryan On 9/6/09 12:51 AM, "David Winsemius" wrote: > I

[R] xyplot {lattice} are different types possible for each panel?

, Bryan * Bryan Hanson Professor of Chemistry & Biochemistry DePauw University, Greencastle IN USA y <- rnorm(100) x <- rnorm(100) names <- rep(c("Set 1", "Set 2", "Set 3"), 4) df <- data.frame(y = y, x = y, names = as.factor(names))

Re: [R] xyplot {lattice} are different types possible for each panel?

Thanks Baptiste, your suggestion works wonderfully. Bryan For anyone following along, the following line needs to replace the similar one in my original example: names <- rep(c("Set 1", "Set 2", "Set 3", "Set 4"), 25) Or the data lengths will be wr

Re: [R] Recursion is slow

Thanks Martin! It seems to be right on and it is blazing fast! I greatly appreciate the responses from you and Bill both. As a beginning user of R it is really helpful to be able to compare my code with yours and try to figure out your tricks. Bryan - Bryan Keller, Doctoral

[R] Slow enumeration of vectors with restrictions

u] } } R } T <- rev(sort(T)) m <- as.integer(m) offset <- sum(T[-1]) - m + 1L nrow <- length(T) - 1L memo <- matrix(rep(NA_real_, nrow * (offset + m)), nrow=nrow) C(length(T), m) } #end of function A - Bryan Kell

[R] How to quickly find the position of the first non-zero element of a vector

31 Or, might the fastest way be to choose between the two methods on-the-fly based on length of the vector, etc.? Bryan - Bryan Keller, Doctoral Student/Project Assistant Educational Psychology - Quantitative Methods The University of Wisconsin - M

[R] Teasing out logrank differences *between* groups using survdiff or something else?

don't quite see it and wouldn't trust ourselves anyway, given the special nature of survival analysis. Manual instructions or a function suggestion would be great. Thanks in Advance, Bryan * Bryan Hanson Professor of Chemistry & Biochemistry

Re: [R] Teasing out logrank differences *between* groups using survdiff or something else?

ore hazardous than T1 vs WISO If in fact it is true? Maybe the answer is already in our output, in the sense that the CI's don't overlap much? Maybe we are wrong to seek a p value as well? Thanks again, Bryan * Bryan Hanson Professor of Chemistry & Biochemistry DePauw U

Re: [R] wilcox.test p-value = 0

exact p-values for sample sizes larger than 100 in each group! You can find it at http://pages.cs.wisc.edu/~klotz/ Bryan > Hi Murat, > I am not an expert in either statistics nor R, but I can imagine that since > the > default is exact=TRUE, It numerically computes the probab

Re: [R] generating unordered combinations

s,dis,dispr,numeric(lws),numeric(lws)) mat3[,4] <- cumsum(mat3[,3]) mat3[,5] <- cumsum(mat3[,3][lws:1])[lws:1] colnames(mat3) <- c("W","Freq","Probability","Sum up","Sum down") print(mat3) } > system.time(f(20,20)) user system

Re: [R] Combine vectors in order to form matrixes with combn

Moreno, I don't understand exactly what it is you are trying to do. Can you explain what you want your matrices to look like? Perhaps give an example. Bryan - Bryan Keller, Doctoral Student/Project Assistant Educational Psychology - Quantitative Methods The Universi

[R] Why doesn't vec[-real.number] give an error or warning? Kids do the darndest things!

behavior be valuable? I can see that it avoids doing floor or as.integer in computations where an index is needed but real numbers are generated for whatever reason, but is that common? Thanks for any insight. Bryan ***** Bryan Hanson Acting Chair Professor of Chemistry & Bioc

Re: [R] Why doesn't vec[-real.number] give an error or warning? Kids do the darndest things!

take in my thinking! Bryan On 3/13/10 11:39 PM, "jim holtman" wrote: > That notation  'v[-1.18]' is not going to remove the 12th entry if it was > 1.18.  They should have used: >   > v[-which(v == 1.18)] >   > to remove it, but it probably would not have wor

[R] How do I define the method for gcheckboxgroup in gWidgets?

. For the sliders, it's clear to me how the actions and drawing of the widgets differ, but not so for gcheckboxgroup. A big TIA, Bryan * Bryan Hanson Professor of Chemistry & Biochemistry DePauw University, Greencastle IN USA Full Script: x <- 1:10 y1 <- x y2 <- x^2 y3 &

Re: [R] How do I define the method for gcheckboxgroup in gWidgets?

Thanks Michael... I was working by analogy to the gbuttons, so I was trying to ³add² the gcheckboxgroup, which is apparently not necessary (due to, I guess, the intrinsic differences between the widgets). The index thing I was just screwed up on! I have it working now. Nice package. Bryan

[R] Managing Packages: Which functions call other functions in package?

might provide; maybe that¹s where I should be looking? Does such a thing exist? Thanks, Bryan * Bryan Hanson Professor of Chemistry & Biochemistry DePauw University, Greencastle IN USA [[alternative HTML version deleted]] __ R-

[R] panel.lmline - are m, b, and r^2 accessible somehow?

done manually within panel.groups ( http://www.nabble.com/add-trend-line-to-each-group-of-data-in%3A-xyplot(y1%2 By2-~-x-|-grp...-td3344023.html#a3382909) but that was a few versions back. Other suggestions? Thanks, Bryan * Bryan Hanson Professor of Chemistry & Biochemistry De

[R] str(data.frame) after subsetting reflects original structure, not subsetted structure?

() a lot to check my work, and in this case, it seems to have misled me. Thanks as always, Bryan * Bryan Hanson Professor of Chemistry & Biochemistry DePauw University, Greencastle IN USA __ R-help@r-project.org mailing list https:

Re: [R] str(data.frame) after subsetting reflects original structure, not subsetted structure?

Thanks Marc and Ben... Your answers were most helpful. I suspected something had been written about it, but was having trouble formulating a reasonable search query. I was looking in the help page for str(), which was sort of a dead end. Bryan * Bryan Hanson Professor of Chemistry

[R] What is mclust up to? Different clusters found if x and y interchanged

ers (not given here) and interchanged the x and y columns and got the same clusters. I'm at a loss. Hopefully someone can set my understanding on the right path. Thanks for any insight! Bryan > sessionInfo() R version 2.10.1 (2009-12-14) x86_64-apple-darwin9.8.0 locale: [1] en_US.UTF-8/

Re: [R] What is mclust up to? Different clusters found if x and y interchanged

much. Bryan On 4/20/10 11:41 AM, "Chris Fraley" wrote: > > > This is a really interesting case, but it is probably not that unusual. > > mod23 <- Mclust(tur[,2:3]) > mod32 <- Mclust(tur[,3:2]) > > mod23[c("modelName","G")] > mod32[

[R] Permutation tests using apply function with package coin

e matrix coin ends up taking longer. My guess is that turning each column into a data frame is what slows coin down. Any ideas on how to tune this up would be most appreciated. Bryan --- Bryan Keller, Doctoral Student/Project Assistant Educational Psychology - Quantitat

[R] ggplot2: proper use of facet_grid inside a function

nnet_7.2-48mclust_3.2 MASS_7.2-48 [21] lars_0.9-7 e1071_1.5-19 class_7.2-48 loaded via a namespace (and not attached): [1] cluster_1.12.0 Thanks for any help! Bryan * Bryan Hanson Professor of Chemistry & Biochemistry DePauw University,

Re: [R] ggplot2: proper use of facet_grid inside a function

;t work as I thought it would (like a hypothetical facet_grid_string()). Thanks so much. Bryan On 10/5/09 4:12 AM, "ONKELINX, Thierry" wrote: > Dear Bryan, > > In the ggplot() function you can choose between aes() and aes_string(). > In the first you need to hardwire t

[R] ggplot2: mapping categorical variable to color aesthetic with faceting

"set" the colors: Do I compute the length of all possible combos of A, B with lrg, sm, and then create one long vector of colors for the entire plot? I tried something like this, and was not successful, but perhaps could be with more work. All advice appreciated, Bryan (session info

Re: [R] ggplot2: mapping categorical variable to color aesthetic with faceting

orking up replaces a much much longer code in base graphics, so I am really liking the thought put into ggplot2 and the leanness of it - Thanks Hadley! Thanks again, Bryan On 10/6/09 12:36 PM, "baptiste auguie" wrote: > Hi, > > I may be missing an important design decision, but

Re: [R] ggplot2: mapping categorical variable to color aesthetic with faceting

s = theme_blank()) And then depending up on the method specified by the user, additional geoms are added and the plot created. This gets the job done, but if there are further suggestions, I'd love to learn other solutions. Bryan On 10/6/09 1:08 PM, "baptiste auguie" wrote: >

[R] General means of matching a color specification to an official R color name

search for the n closest matches of hex in df.rgb$hex.code # perhaps hex should be converted into a different color space 1st # eventually would like to display matches side by side w/hex } Thanks as always. Bryan * Bryan Hanson Acting Chair Professor of Chemistry & Bi

Re: [R] General means of matching a color specification to an official R color name

Works perfectly! Thanks Barry. I had actually seen some suggestions on using a distance, but by then I was thinking about hcl spaces and distance isn't so as simple there. I'm too tired I think. Anyway, you've got me running again! Thanks, Bryan On 10/13/09 5:12 PM, &q

[R] Subsetting/modifying a symbolic formula

rms.formula, ?formula etc all discuss these attributes but as I said, it seems a bit impenetrable. RSiteSearch("symbolic formula") returns too many answers. No doubt I am missing the obvious, as this is my first foray into using formulas in my own functions. TIA, Bryan * Bryan H

Re: [R] Subsetting/modifying a symbolic formula

Thanks Gabor, you taught me two useful things: all.vars, and the fact that a formula object always has length 3. Problem solved. Bryan On 10/21/09 11:57 AM, "Gabor Grothendieck" wrote: > Here is one way: > >> fo <- y ~ f1 * f2 >> one.x <- lapply(all.

Re: [R] Building an R package on MAC

Gina, at the terminal, make sure you are in directory Users/apple/Docments/R When you type >R CMD build TEST Sounds like you might be in the TEST directory, not the R directory. The build system looks for a subdirectory with your package name when it starts. Bryan ***** Bryan Han

[R] ggplot2: stat_bin ..count.. with geom_text when NA is present

is carried out as the last step. So, any suggestions about why the counts don't appear on my plot? I suppose I can always clean the data first, but it would be much more practical to do that in the background during the preparation of the plot. Thanks as always, Bryan * Bryan

Re: [R] ggplot2: stat_bin ..count.. with geom_text when NA is present

.. Adding to my original post... OK, here's a little function which demonstrates the behavior I described. Try it with rem = FALSE to see the annotation, then TRUE to see the annotations disappear. What's going on here? Thanks, Bryan res = runif(50, 0, 100) fac = rep(c("A

Re: [R] PDF Corrupted?

Looks like you need to do dev.off() after the plot to properly close the file. Bryan * Bryan Hanson Acting Chair Professor of Chemistry & Biochemistry DePauw University, Greencastle IN USA On 10/27/09 8:42 PM, "rkevinbur...@charter.net" wrote: > I am running R 2

Re: [R] ggplot2: stat_bin ..count.. with geom_text when NA is present

milar. Revised function below; n = 50 appears at the bottom of the plot regardless of how many values there are... Hmmm JUST FIXED IT: added the df <- na.omit(df) and the counts are correct! Thanks, Bryan res = runif(50, 0, 100) fac = rep(c("A", "B"), 50) df <- data.frame

Re: [R] plot.window arguments being ignored?

s automatically resize as necessary. Both lattice and ggplot2 have books and very nice web sites where you can find an example like yours. Bryan * Bryan Hanson Acting Chair Professor of Chemistry & Biochemistry DePauw University, Greencastle IN USA On 11/7/09 6:28 PM, "A

[R] Where are usages like "== 2L" documented?

Gurus: I keep seeing other people¹s code that contain ideas like If (x == 2L) X[-1L] X - 1L I have some idea of what¹s going on, but where is the use of concepts like ³2L² documented? Thanks, Bryan * Bryan Hanson Acting Chair Professor of Chemistry & Biochemistry DePauw Univer

Re: [R] Where are usages like "== 2L" documented?

at least not in the first few pages. Thanks for the tips. Bryan * Bryan Hanson Acting Chair Professor of Chemistry & Biochemistry DePauw University, Greencastle IN USA On 11/16/09 7:52 PM, "Gabor Grothendieck" wrote: > On Mon, Nov 16, 2009 at 7:25 PM, Duncan Murdo

Re: [R] Where are usages like "== 2L" documented?

explictly an integer for whatever reason? Are there other reasons, for instance, ways it saves lines of code? Bryan * Bryan Hanson Acting Chair Professor of Chemistry & Biochemistry DePauw University, Greencastle IN USA On 11/17/09 4:20 AM, "Patrick Burns&qu

[R] rgl keyboard shortcut for translation on Mac?

seem to modify anything. So, what I can't find is a key combination that translates the object. From reading various documents about openGL, it seems there might be a toggle somewhere, but I need a hint! Or do I need to install/invoke some additional command? Thanks, Bryan * Br

Re: [R] rgl keyboard shortcut for translation on Mac?

Thanks Duncan... I see how to do it. It appears that all 3 of my mouse buttons are actually set to the trackpad button on my MacBook Pro. So if I assign any button via pan3d I get the same behavior, and can't get back to the default. I'll play with it! Bryan On 12/6/09 9:14 P

Re: [R] Questíon regarding the use of write.csv2 , write.table ...

Write.table will give you all the control you need to get exactly what you want. Bryan On 6/18/09 7:50 AM, "xavier.char...@free.fr" wrote: > Hi, It sounds like the first column that is "added" is actually the row > names. That's why a previous answer point

[R] recursive function, how to avoid list structure in return value

ion could be fixed. Can anyone show me show? Thanks, Jenny Bryan Demo of my problem -- (Note: although my question has nothing to do with hierarchical clustering per se, my example does assume knowledge of the 'merge' object.) ## faking the key aspects of an hclust object myClus

[R] R-Squared with biglm?

I've been working with a rather large data set (~10M rows), and while biglm works beautifully for generating coefficients, it does not report an r-squared. It does report RSS. Any idea on how one could coax an R-squared out of biglm? Thanks in advance for any help with this! Brya

[R] ggplot2, building a simple formula interface

f1 <- rep(c("C", "D", "D", "C"), 25) f2 <- rep(c("A", "B"), 50) mydata <- data.frame(res, f1, f2) df <- simple(~ res + f1, mydata) p <- ggplot(df, aes(f1, res)) + geom_boxplot() plot(p) Thanks, Bryan * Br

[R] Comparing percentages

difference of percentages test, but I am not sure how to account for the third option. Any help is appreciated. Bryan [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] graph help

o the lattice or ggplot2 graphics systems, which handle groups of data related by some categorical variable much better. Look at http://had.co.nz/ggplot2/ for some ideas of what's possible, or http://lmdvr.r-forge.r-project.org/figures/figures.html?chapter=01;figure=;t heme=stdColor;code=ri

[R] Wordcloud warnings

ewer.pal(9,"Reds")[0:-3]) } Best, Bryan __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

[R] Wordcloud errors

s there a way to fix this script so that there are no warnings?? Thanks, Bryan [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do re

Re: [R] (no subject)

Look at the package chemometrics, it can certainly handle your number of variables (p > n is what that's called and it requires special considerations). I don't recall about missing values. The authors of that package also have a very helpful text. Good Luck. Bryan On Jul 14, 2

[R] Breaking out of multiple loops

tatement in the outer loop, but again that's inefficient, as it's checking that conditional hundreds of times.) So is there a way to "cleanly" break out of multiple loops? Thanks, -bryan __ R-help@r-project.org mailing list https:

Re: [R] Breaking out of multiple loops

head, or if functions are innately faster somehow. Still seems like there should be a way to break out of nested loops, though... -b On Tue, Dec 18, 2012 at 4:50 PM, Duncan Murdoch wrote: > On 12-12-18 1:02 PM, McCloskey, Bryan wrote: >> >> Hey all, >> >> I'm curren

Re: [R] Fitting this data with a gaussian would be great

Fortune candidate? I hear the landlord is hell, but the company good. Bryan > I've > already got an apartment reserved for me in one of Pat Burns's "R > Inferno" levels, and I don't want to descend even further. > > Best, > Bert

Re: [R] plotting CIE chromaticity diagram?

I can send you my not quite perfect attempt a little later today. Bryan **** Prof. Bryan Hanson Dept of Chemistry & Biochemistry DePauw University Greencastle IN 46135 USA academic.depauw.edu/~hanson/deadpezsociety.html github.com/bryanhanson academic.depauw.edu/~hanson/UMP/In

Re: [R] plotting CIE chromaticity diagram?

. Bryan On Mar 18, 2013, at 7:08 AM, Ken Knoblauch wrote: > ishi soichi gmail.com> writes: > >> Has anyone plotted or is it possible to plot >> >> CIE *xy* chromaticity diagram >> >> http://en.wikipedia.org/wiki/File:CIE1931xy_blank.svg >> >> I

Re: [R] OrgMassSpecR peak area issue

sity[k] - peakTime[k] * peakIntensity[j] } peakArea[i] <- abs(sum(x)/2) which looks pretty standard to me, though I'm not clear right off the top of my head why they are dividing by 2. You can always contact the maintainer. Bryan On Mar 18, 2013, at 1:34 P

Re: [R] Image segmentation

I've recently had a reason to work a little with image segmentation too, and in addition to EBImage, you should look at biOps. You can learn a lot by studying these packages. Bryan On Apr 6, 2013, at 10:04 AM, Eder Paulo wrote: > Hello, > > Thanks for replying me! > >

Re: [R] Linebreaks in cat() functions that call other variables?

The only permutation you likely didn't try: >> cat("df(between) is", a[1,1], "\ndf(within) is", a[2,1]) \n has to be inside the quotes. HTH. Bryan On May 22, 2013, at 4:34 PM, jordanbrace wrote: > Hi everyone, > > I'm having some difficul

[R] rgl crashes after one successful draw

cessfully. And the same thing happens if I just do the commands inside the function individually. SessionInfo() below. Perhaps something is corrupt with my X11 window system? Thanks, Bryan R version 3.0.1 (2013-05-16) Platform: x86_64-apple-darwin10.8.0 (64-bit) locale: [1] en_US.UTF-8

Re: [R] rgl crashes after one successful draw

Thank you Brian. Yes, problem is on R.app I will send to R-sig-mac. Thanks, Bryan On May 29, 2013, at 11:07 AM, Prof Brian Ripley wrote: > This is most relevant to R-sig-mac. There are two different rgl devices on > OS X, depending how you are running this. One based on X11 and

[R] Proper way to implement package internal functions

e seen this somewhere, but I can't find any resources just now. Perhaps a suggestion of a package that does things this way which I could study would be sufficient. Thanks, Bryan __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/li

Re: [R] Proper way to implement package internal functions

hould have documentation entries. if I understand correctly. I guess the reason I didn't find any documentation is the wide lattitude which is possible. Thank you. Bryan On Jun 12, 2013, at 10:57 AM, Duncan Murdoch wrote: > On 12/06/2013 10:44 AM, Bryan Hanson wrote: >> [p

[R] grDevices::convertColor XYZ space is it really xyY?

e CIE tristimulus values? It looks to me like the first case is true, but I would appreciate hearing from one of the people in the know. Thanks, Bryan __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] grDevices::convertColor XYZ space is it really xyY?

, from = "XYZ", to = "sRGB") # none of these are 1,1,1, namely white (they are ~ 0.6, 0.6, 0.6) So it looks like D65, a white standard, does not come back to something near white in the sRGB space. What am I doing wrong here, or what do I misunderstand? Please don't s

Re: [R] grDevices::convertColor XYZ space is it really xyY?

y are tristimulus values - is that correct? Thanks again. This solves several problems in a package I am developing. Bryan On Jun 12, 2013, at 4:36 PM, Ken Knoblauch wrote: > You seem to treating the input values as xyY when they should be XYZ > (case matters). > So, I would do somet

Re: [R] Proper way to implement package internal functions

hide, that the warning goes away. Since it is less work to change the export statement compared to even a minimal Rd, that's the way I went. It's interesting that there is not more info about these options available. Thanks, Bryan On Jun 12, 2013, at 6:46 PM, Rolf Turner wrote: >

  1   2   3   >