Re: [R] BUG: atan(1i) / 5 = NaN+Infi ?

On Fri, Sep 13, 2024 at 1:10 PM Duncan Murdoch wrote: > [You don't often get email from murdoch.dun...@gmail.com. Learn why this is > important at https://aka.ms/LearnAboutSenderIdentification ] > Caution: External email. > On 2024-09-13 8:53 a.m., Jonathan Dushoff wrote

Re: [R] BUG: atan(1i) / 5 = NaN+Infi ?

> Message: 4 > Date: Thu, 12 Sep 2024 11:21:02 -0400 > From: Duncan Murdoch > That's not the correct formula, is it? I think the result should be x * > Conj(y) / Mod(y)^2 . Correct, sorry. And thanks. > So that would involve * and > / , not just real arithmetic. Not an expert, but I don't s

[R] Subject: Re: BUG: atan(1i) / 5 = NaN+Infi ?

t might also make sense to calculate x / y using real arithmetic (as x*y / |y|²) Jonathan __ 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 https://www.R-proj

[R] R

Hi I found your email on a website Can I ask some questions about R please Many thanks Jonathan [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo

[R] Best settings for RStudio video recording?

Folks: I was wondering if you all would suggest some helpful RStudio configurations that make recording a session via e.g. zoom the most useful for students doing remote learning. Thoughts? --j -- Jonathan A. Greenberg, PhD Randall Endowed Professor and Associate Professor of Remote Sensing

[R] Checking for a proper "stop" statement...

27;m having a hard time figuring out how to do this test. --j -- -- Jonathan A. Greenberg, PhD Randall Endowed Professor and Associate Professor of Remote Sensing Global Environmental Analysis and Remote Sensing (GEARS) Laboratory Natural Resources & Environmental Science University of Nevada, R

[R] [R-pkgs] New package: rDotNet

rg/web/packages/rDotNet/index.html <https://cran.r-project.org/web/packages/rDotNet/index.html> The package is stable, as has been in use for some years, but only now packaged up for public use on CRAN. Feel free to contact with questions or suggestions on GitHub or by email. Regards -- Jo

Re: [R] Odd results from rpart classification tree

replicate(1000, test_split(0.5))) # 2-12, i.e. splits only sometimes... Adding a constant to y and getting different trees is a bit strange, particularly stochastically. Will see if I can track down a copy of the CART book. Jonathan From: Therneau, Terry

[R] Odd results from rpart classification tree

t <- gini(prop.table(table(y))) impurity_L <- gini(prop.table(table(obs_L))) impurity_R <- gini(prop.table(table(obs_R))) impurity <- impurity_root * n - (n_L*impurity_L + n_R*impurity_R) # 2.880952 Thus, an improvement of 2.88 should result in a split. It does not. Why? Jonathan __

[R] [R-pkgs] New package: ggghost 0.1.0 - Capture the spirit of your ggplot2 calls

me any feedback or suggestions you may have. Kind regards, - Jonathan Carroll. [[alternative HTML version deleted]] ___ R-packages mailing list r-packa...@r-project.org https://stat.ethz.ch/mailman/listinfo/r-pa

[R] 2x2x2 rm ANOVA, varying results

Hello, I ran a 2x2x2 repeated measures ANOVA which turned out fine: DfSum Sq Mean Sq F value Pr(>F) Attend1 0.5540 0.55402 7.03740.01079 *PercGrp 1

[R] Problems with pooling Multiply Imuputed datasets, of a multilevel logistic model, using (MICE)

I am having problems with the MICE package in R, particularity with pooling the imputed data sets. I am running a multilevel binomial logistic regression, with Level1 - topic (participant response to 10 questions on different topics, e.g. T_Darkness, T_Day) nested within Level2 - individuals. T

[R] R Licensing Question

I could not find an answer elsewhere. Thank you, Jonathan __ Jonathan Gellar Statistician Mathematica Policy Research 1100 First Street NE, 12th Floor Washington, DC 20002 __ R-help@r-project.org mailing list -- To

Re: [R] Replace NaN with value from the same row

- > Sent from my phone. Please excuse my brevity. > > On October 18, 2015 12:24:42 PM PDT, Jonathan Reardon > wrote: > >Hi, Sorry to be a pain. Would you be kind enough to briefly explain > >what the lines

[R] Replace NaN from 1 column with a value from the same row

Hi everyone, Ignore my previous post, i realised that the rows and columns i typed into the email were unreadable, sincere apologies for this. A simple question, but i cannot figure this out. I have a data-frame with 4 columns (onset, offset, outcome, mean): df<-data.frame(onset=c(72071,142598,293

[R] Replace NaN with value from the same row

eplace any mean row NaN's with the offset value from that very same row? I don't want to use any pasting etc as this needs to be used as part of a function working over a large dataset than the one shown here. Cheers Jonathan

Re: [R] kknn::predict and kknn$fitted.values

oses kknn$best.parameters, but want to verify this. Hopefully that clarifies the issue. I post here in case future users have a similar question. Thanks to any who took the time to think about this! Jonathan -- View this message in context: http://r.789695.n4.nabble.com/kknn-predi

[R] kknn::predict and kknn$fitted.values

2 3 5 6 7 1 70 0 0 0 0 0 2 0 76 0 0 0 0 3 0 0 17 0 0 0 5 0 0 0 13 0 0 6 0 0 0 0 9 0 7 0 0 0 0 0 29 Can anyone clarify what fitted.values and predict actually do? I would have expected they would give the same output. Thanks... Jonathan -- View t

[R] R command to open a file "browser" on Windows and Mac?

Folks: Is there an easy function to open a finder window (on mac) or windows explorer window (on windows) given an input folder? A lot of times I want to be able to see via a file browser my working directory. Is there a good R hack to do this? --j [[alternative HTML version deleted]]

Re: [R] Correlation question

a linear function of speed! Hence > cor(cars$dist,fitted.right)^2 and cor(x=cars$dist,y=fitted.wrong)^2 must be > the same. > > HTH > d > > Feladó: R-help [r-help-boun...@r-project.org] ; meghatalmazó: Jonathan > Thayn [jth...@il

[R] Correlation question

I recently compared two different approaches to calculating the correlation of two variables, and I cannot explain the different results: data(cars) model <- lm(dist~speed,data=cars) coef(model) fitted.right <- model$fitted fitted.wrong <- -17+5*cars$speed When using the OLS fitted values, the

[R] prediction intervals for robust regression

for robust regression, I haven't had much success in finding out how to create prediction intervals for the results. I was wondering if anyone would be able to provide some direction on how to create these prediction intervals in the robust regression setting. Thanks, Jonathan Bur

Re: [R] Using PCA to filter a series

This is exactly what I was looking for. Thank you. Jonathan Thayn On Oct 3, 2014, at 10:32 AM, David L Carlson wrote: > You can reconstruct the data from the first component. Here's an example > using singular value decomposition on the original data matrix: > >> d <

Re: [R] Using PCA to filter a series

p. It's been a long time since I've played with PCA. Jonathan Thayn On Oct 2, 2014, at 4:59 PM, David L Carlson wrote: > I think you want to convert your principal component to the same scale as d1, > d2, d3, and d4. But the "original space" is a 4-dimensional space

[R] Using PCA to filter a series

to have created the clean pattern I want, but I would like to project the first component back into the original axes? Is there a simple way to do that? Jonathan B. Thayn [[alternative HTML version deleted]] __ R-help@r-project.org mai

Re: [R] Building R for better performance

All, I’ve attached the actual benchmark TACC and I used. I’ve also attached a paper I wrote covering this in a little more detail. The paper specifies the hardware configuration I used. Let me know if you have any other questions. Regards, Jonathan Anspach Sr. Software Engineer Intel Corp

Re: [R] Building R for better performance

Yes, that's the original. Then TACC increased the matrix sizes for their tests. Jonathan Anspach Intel Corp. Sent from my mobile phone. On Sep 11, 2014, at 9:18 AM, "Henrik Bengtsson" mailto:h...@biostat.ucsf.edu>> wrote: You'll find R-benchmark-25.R, which I

Re: [R] Building R for better performance

I'm out of the office today, but will resend it tomorrow. Jonathan Anspach Intel Corp. Sent from my mobile phone. On Sep 11, 2014, at 3:49 AM, "arnaud gaboury" wrote: >>> I got the benchmark script, which I've attached, from Texas Advanced >>> Computing

[R] table over a matrix dimension...

ut requires two nested loops (apply and for() ), so # not efficient for very large datasets. ### Is there a more elegant solution to this? --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmental Analysis and Remote Sensing (GEARS) Laboratory Department of Geography and Ge

[R] StatET and R 3.1.0

arting the R engine failed. Please make sure that R package 'rj' (1.1 or compatible) is installed and that the R library paths are set correctly for the R environment configuration 'R'. --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmental Analysis

Re: [R] Ignore escape characters in a string...

gineer (Solar/BatteriesO.O#. #.O#. with > /Software/Embedded Controllers) .OO#. .OO#. rocks...1k > --- > Sent from my phone. Please excuse my brevity. > > On April 8, 2014 8:

[R] Ignore escape characters in a string...

k to using a copy/paste from Windows explorer I'm not aware of? --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmental Analysis and Remote Sensing (GEARS) Laboratory Department of Geography and Geographic Information Science University of Illinois at Urbana-Champaign 25

[R] numeric to factor via lookup table

ne the correct levels (ID matches the values), and code is the label? --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmental Analysis and Remote Sensing (GEARS) Laboratory Department of Geography and Geographic Information Science University of Illinois at Urbana-Champai

[R] Determine breaks based on a break type...

the break component: mybreaks <- hist(runif(100))$breaks But is there a faster way to do this, if this is all I want? --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmental Analysis and Remote Sensing (GEARS) Laboratory Department of Geography and Geographic Informat

Re: [R] Overriding predict based on newdata...

data' and the the definition of the new > function would persist through the rest of the session and could be > source()-ed in further sessions. > > -- > David. > > > On Mar 16, 2014, at 2:09 PM, Jonathan Greenberg wrote: > >> R-helpers: >> >>

[R] Overriding predict based on newdata...

g) is that predict is dispatched based on the object parameter, this is causing me confusion -- my object should still remain the model, I'm just allowing a new data type to be fed into the predict model(s). Cheers! --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmen

Re: [R] Building R for better performance

CE, Sorry for the delay. I haven't installed any additional packages, so I don't know the answer to your question. Let me look into it and get back to you. Regards, Jonathan Anspach Sr. Software Engineer Intel Corp. jonathan.p.ansp...@intel.com 713-751-9460 -Original Message

Re: [R] Building R for better performance

6.02 Total 274.93 21.01 Regards, Jonathan Anspach Sr. Software Engineer Intel Corp. jonathan.p.ansp...@intel.com 713-751-9460 -Original Message- From: Simon Zehnder [mailto:szehn...@uni-bonn.de] S

[R] Building R for better performance

ion instructions and FAQ? I can post my data if anyone is interested. Thanks, Jonathan Anspach Sr. Software Engineer Intel Corp. jonathan.p.ansp...@intel.com 713-751-9460 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-he

[R] install packages from R-forge SVN

eers! --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmental Analysis and Remote Sensing (GEARS) Laboratory Department of Geography and Geographic Information Science University of Illinois at Urbana-Champaign 259 Computing Applications Building, MC-150 605 East Springfield Avenue Cha

[R] Checking for and adding "..." arguments to a function...

ification should be done without invoking the function, e.g. I'm not asking how to test for this WITHIN the function, I'm asking how to test "myfunction" directly as an R object. Thanks! --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmental Analysis and

Re: [R] Official way to set/retrieve options in packages?

rvey.ultimate.cluster = FALSE) >> if (is.null(getOption("survey.want.obsolete"))) >> options(survey.want.obsolete = FALSE) >> if (is.null(getOption("survey.adjust.domain.lonely"))) >> options(survey.adjust.domain.lonely = FALSE) >> if

[R] Multimodal multidimensional optimization

e any suggestions, I'm quite a bit beyond my expertise here. Thanks. Jonathan B. Thayn, Ph.D. Ridgely Fellow of Geography Department of Geography – Geology Illinois State University Felmley Hall of Science, Rm 200A Normal, IL 61790 jth...@ilstu.edu my.ilstu.edu/~jthayn [[altern

Re: [R] Error: C stack usage is too close to the limit when using list.files()

detected; `dir/subdir/linkToUpperDir' is part of the > same file system loop as `dir'. > dir/IMPORTANT_1 > > Bill Dunlap > Spotfire, TIBCO Software > wdunlap tibco.com > > >> -----Original Message- >> From: r-help-boun...@r-project.org [mailto:r-h

Re: [R] Error: C stack usage is too close to the limit when using list.files()

ot;,...) That'll teach me not to copy/paste from my code... On Fri, Sep 27, 2013 at 2:36 PM, Ben Bolker wrote: > Jonathan Greenberg illinois.edu> writes: > >> >> R-helpers: >> >> I'm running a file search on my entire drive (Mac OS X) using: >>

[R] Error: C stack usage is too close to the limit when using list.files()

eturns ~ 3 results). I keep getting an error: Error: C stack usage is too close to the limit when running this command. Any ideas on 1) how to fix this or 2) if there is an alternative to using list.files() to accomplish this search without resorting to an external package? Cheers! --jonathan

Re: [R] Confusing behaviour in data.table: unexpectedly changing variable

Thanks for your help, and sorry for mis-posting. JD On Wed, Sep 25, 2013 at 3:18 AM, Matthew Dowle wrote: > Very sorry to hear this bit you. If you need a copy of names before > changing them by reference : > oldnames <- copy(names(DT)) > This will be documented and it's on the bug list

[R] Confusing behaviour in data.table: unexpectedly changing variable

I got bitten badly when a variable I created for the purpose of recording an old set of names changed when I didn't think I was going near it. I'm not sure if this is a desired behaviour, or documented, or warned about. I read the data.table intro and the FAQ, and also ?setnames. Ben Bolker crea

Re: [R] R-3.0.1 g77 errors

Sadly, I am limited to the Solaris 10 system. I wish that I could use Linux, the world uses it. -- Jonathan M. Prigot Partners Healthcare Systems On Wed, 2013-09-18 at 17:20 +0200, Simon Zehnder wrote: > On my systems Linux Scientific and Mac OS X I use as well for the F77 the > gf

[R] R-3.0.1 g77 errors

^ Invalid declaration of or reference to symbol `digits' at (^) [initially seen at (^)] dlamch.f:89: warning: INTRINSIC DIGITS, EPSILON, HUGE, MAXEXPONENT, ^ Reference to unimplemented intrinsic `EPSILON' at (^) (assumed EXTERNAL

[R] library() and install.packages() no longer working ("Access is denied" error)

o create a personal library 'D:/Users/[UID]/Documents/R/win-library/3.0' to install packages into? I clicked yes. The Mirror browser shows up, I select a mirror. A 3.0 directory is created, but I got the same error, and when examining the (new) 3.0 directory, nothing is created inside

[R] ifelse question (I'm not sure why this is working)...

ve never seen an ifelse statement used that way, and I was fully expecting that to NOT work, or to place the output of which(x%%d==0) in each location where the statement x%%d==0 was true. Any ideas on deconstructing this? --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmental Anal

[R] Vectorized version of colMeans/rowMeans for higher dimension arrays?

orized operation (so it doesn't execute as quickly). How would one vectorize this operation (if possible)? Is there an array equivalent of colMeans/rowMeans? --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmental Analysis and Remote Sensing (GEARS) Laboratory De

Re: [R] Parallel version of Map(rather, mapply)

Hi Saptarshi: There are quite a few parallel mapply's out there -- my recommendation is to use the foreach package, since it allows you to be flexible in the parallel backend, and you don't have to write two statements (a sequential and a parallel statement) -- if a parallel backend is running,

Re: [R] Determining the maximum memory usage of a function

happening. > > > On Thu, Jun 20, 2013 at 10:45 AM, Jonathan Greenberg > wrote: >> >> Folks: >> >> I apologize for the cross-posting between r-help and r-sig-hpc, but I >> figured this question was relevant to both lists. I'm writing a >> function to be a

[R] Determining the maximum memory usage of a function

for determining this? --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmental Analysis and Remote Sensing (GEARS) Laboratory Department of Geography and Geographic Information Science University of Illinois at Urbana-Champaign 607 South Mathews Avenue, MC 150 Urbana, IL 61801

Re: [R] Instant search for R documentation

n't think it's > quite ready yet. I think it might be wise to check for it later this year. > > > * The "sos" package with its vignette for searching CRAN > packages and getting the result sorted to place first the package with > the mo

[R] Instant search for R documentation

Hi, I just wanted to share with you that we made a website over the weekend that allows "instant search" of the R documentation on CRAN, see: www.Rdocumentation.org. It's a first version, so any feedback/comments/criticism most welcome. Best rega

Re: [R] Official way to set/retrieve options in packages?

ll(getOption("survey.adjust.domain.lonely"))) > > options(survey.adjust.domain.lonely = FALSE) > > if (is.null(getOption("survey.drop.replicates"))) > > options(survey.drop.replicates = TRUE) > > if (is.null(getOption("survey.multico

[R] Official way to set/retrieve options in packages?

quot; way to do this? I see there is an options() and getOptions() function, but I'm unclear how I would use this in my own package to create/save new options for my particular package. Cheers! --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmental Analysis and R

[R] RCurl: using ls instead of NLST

"ls" to list the directory contents -- is there any way to use "ls" remotely on this site? Thanks! --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmental Analysis and Remote Sensing (GEARS) Laboratory Department of Geography and Geographic Information Science

Re: [R] points overlay axis

> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On > Behalf Of John Kane > Sent: Tuesday, May 14, 2013 7:47 AM > To: Jonathan Phillips; r-help@r-project.org > Subject: Re: [R] points overlay axis > > Probably

[R] R help: Batch read files based on names in a list

* I am currently reading in a series of files, applying the same functions to them one at a time, and then merging the resulting data frames e.g.: >MyRows <- c("RowA", "RowB", "RowC")>>File1_DF <- >read.delim("DirectoryToFiles\\File1_Folder\\File1.txt", >stringsAsFactors=FALSE, check.names=

[R] points overlay axis

Hi, I'm trying to do quite a simple task, but I'm stuck. I've set xaxs = 'i' as I want the origin to be (0,0), but unfortunately I have points that are sat on the axis. R draws the axis over the points, which hides the points somewhat and looks unsightly. Is there any way of getting a point to be

[R] Stepwise regression for multivariate case in R?

iable lengths differ (found for 'grouping') What can be wrong here? I have checked and all variables in my.data is of the same length. //Jonathan [[alternative HTML version deleted]] __ R-hel

Re: [R] Singular design matrix in rq

26677.8 Max. :69.342 fit_spl <- rq(response ~ bs(predictor,df=15),tau=1,data=mydata) # Error in rq.fit.br(x, y, tau = tau, ...) : Singular design matrix --j On Fri, Apr 19, 2013 at 8:15 AM, Koenker, Roger W wrote: > Jonathan, > > This is not what we call a reproducible example..

[R] Singular design matrix in rq

) fit_spl <- rq(raw_data[,1] ~ bs(raw_data[,i],df=15),tau=1) > Error in rq.fit.br(x, y, tau = tau, ...) : Singular design matrix Any ideas what might be causing this or, more importantly, suggestions for how to solve this? I'm just trying to fit a smoothed hull to the top of the data clou

Re: [R] Check if a character vector can be coerced to numeric?

Yep, type.convert was exactly what I was looking for (with as.is=TRUE). Thanks! On Thu, Mar 21, 2013 at 1:31 PM, Prof Brian Ripley wrote: > On 21/03/2013 18:20, Jonathan Greenberg wrote: > >> Given an arbitrary set of character vectors: >> >> myvect1 <- c("abc&

[R] Check if a character vector can be coerced to numeric?

a numeric, but leaves character vectors which cannot be converted (myvect1) alone. Is there any simple way to do this (e.g. some function that tests if a vector is coercible to a numeric before doing so)? --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmental Analysis and Remote

[R] Issue with matrices within nested for-loops

bounds* Which I understand to mean that the bounds of the indicated vector/matrix have been violated. I am however at a loss as to how to resolve this. Any advice would be appreciated Cheers! JR -- Jonathan Richar Doctoral candidate UAF SFOS Fisheries Division 17101 Pt. Lena Loop Rd. Unive

[R] Mixed models with missing data

Hi, I am creating a mixed model based on a experiment where each subject has 2 repeats. In some instances though there is only data for one of a given subjects repeats for most there is data for both. Can I still justify having subject as a random effect? Thanks, Jonathan [X] [X] [X] [X] [X

[R] mixed effects regression with non-independent data

one subject from each game for the analysis (and so losing half the data). Is there a way to introduce this into the model instead, perhaps as a random effect?? Thanks, Jonathan [X] [X] [X] [X] [X] [X] [X] [X] [[alternative HTML version deleted

[R] Inserting rows of interpolated data

Dear help list - I have light data with 5-min time-stamps. I would like to insert four 1-min time-stamps between each row and interpolate the light data on each new row. To do this I have come up with the following code: lightdata <- read.table("Test_light_data.csv", header = TRUE, sep = ",") #

[R] Count of Histogram Bins using Shingles with lattice

I know that I can get a count of histogram bins in base R with plot=FALSE. However, I'd like to do the same thing with lattice. The problem is that I've set up shingles, and I'd like to get the count within each bin within each shingle. plot=FALSE doesn't seem to do it. [[alternativ

Re: [R] Loading a list into the environment

,b=20) > > myfun(myvariables) > a > b > > > Hope this helps, > > Rui Barradas > > Em 01-02-2013 22:24, Jonathan Greenberg escreveu: > > R-helpers: >> >> Say I have a list: >> >> myvariables <- list(a=1:10,b=20) >> >> I

[R] Loading a list into the environment

ent without having to explicitly define them. --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmental Analysis and Remote Sensing (GEARS) Laboratory Department of Geography and Geographic Information Science University of Illinois at Urbana-Champaign 607 South Mathews Avenue, MC

[R] scaling of nonbinROC penalties - accurate classification with random data?

x sum to some multiple of the number of categories), or (c) is one free to use arbitrarily-scaled penalty matrices? I apologise if I am wasting your by making an obvious mistake. I am a clinician, not a statistician. So, I do not understand the m

[R] scaling of nonbinROC penalties - accurate classification with random data?

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

Re: [R] Adding a line to barchart

t; > Petr > > > -Original Message- > > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > > project.org] On Behalf Of Jonathan Greenberg > > Sent: Tuesday, January 22, 2013 11:42 PM > > To: r-help > > Subject: [R] Adding a line to barchart &

[R] Adding a line to barchart

positioned at x = 2 that lays over the bars (say, using a dotted line). How do I go about doing this? --j -- Jonathan A. Greenberg, PhD Assistant Professor Global Environmental Analysis and Remote Sensing (GEARS) Laboratory Department of Geography and Geographic Information Science University

[R] scaling of nonbinROC penalties

or (c) is one free to use arbitrarily-scaled penalty matrices for estimates of the accuracy of an ordinal gold standard? Thanks, in advance, for your help, Jonathan Williams

[R] Best way to coerce numerical data to a predetermined histogram bin?

is example, the bins are not necessarily equal widths. I can, of course, cycle through each element of data, and then move through breaks, stopping when it finds the correct bin, but I feel like there is probably a faster (and more elegant) approach to this. Thoughts? --j -- Jonathan A.

[R] Anomalous outputs from rbeta when using two different random number seeds

Hi, in the code below, I am drawing 1000 samples from two beta distributions, each time using the same random number seed. Using set.seed(80) produces results I expect, in that the differences between the distributions are very small. Using set.seed(20) produces results I can't make sense of. Aro

[R] How to change smoothing constant selection procedure for Winters Exponential Smoothing models?

three smoothing constants. Can anyone help? Thank you, Jonathan [[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.or

Re: [R] Selecting the "non-attribute" part of an object

all.equal(m,mm,check.attributes=FALSE) > gives TRUE!!! I.e. sometimes attributes really are vital characteristics. > cheers, > Rolf Turner > On 16/11/12 08:52, Jonathan Dushoff wrote: >> I have two matrices, generated by R functions that I don't understa

[R] Selecting the "non-attribute" part of an object

I have two matrices, generated by R functions that I don't understand. I want to confirm that they're the same, but I know that they have different attributes. If I want to compare the dimnames, I can say > identical(attr(tm, "dimnames"), attr(tmm, "dimnames")) [1] FALSE or even: > identical(d

[R] Best R textbook for undergraduates

loop is, an if-then statement, etc). Suggestions? In particular, I'd like to hear from those of you who have TAUGHT classes using R. Thanks! --jonathan -- Jonathan A. Greenberg, PhD Assistant Professor Department of Geography and Geographic Information Science University of Illinois at U

[R] Can I please be taken off the mailing list

[[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 and provide commented, minimal, self-contained,

[R] Completely ignoring an error in a function...

creates two files, failing on the second of the two files -- but the first file it creates is what I wanted and there is no current way to create that single file on its own without a lot of additional coding). --j -- Jonathan A. Greenberg, PhD Assistant Professor Department of Geography and Geogr

Re: [R] Proposal: Package update log

Thanks to all for the responses and suggestions. I was primarily proposing a more detailed change log for packages on CRAN. To my mind, repositories like R-forge host packages more 'raw' than those on CRAN (i.e. CRAN seems to me to contain more 'finished' packages which occasionally are updated

[R] Proposal: Package update log

I'm relatively new to R and would first like to sincerely thank all those who contribute to its development. Thank you. I would humbly like to propose a rule which creates a standard (i.e., strongly encouraged, mandatory, etc.) for authors to include a `change log' documenting specific changes

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

ge: > In writeChar(character(1), fl, nchars = 1, useBytes = TRUE) : > writeChar: more characters requested than are in the string - will > zero-pad > > close(fl) > > > File "Test.txt" is now 1Kb in size. > > Hope this helps, > > Rui Barradas > Em 27-09-2012

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

y not larger than the filesystem can handle). However, length appears to be restricted by .Machine$integer.max (I'm on a 64-bit windows box): > .Machine$integer.max [1] 2147483647 Any suggestions on how to solve this problem for much larger file sizes? --j On Thu, May 3, 2012 at 10:44

Re: [R] list of funtions

t using some 'i' like variable in a loop. I don't think it's something obvious I've missed... On 13 September 2012 18:06, Uwe Ligges wrote: > > > On 13.09.2012 19:01, Jonathan Phillips wrote: >> >> Hi, >> I have a function called fitMicroProtei

[R] list of funtions

Hi, I have a function called fitMicroProtein which takes a value called form, this can be any integer from 0-38. In a larger function I'm making (it's called Newton), the first thing I want to do is construct a list of functions where form is already set. So in pseudocode fs[[1]](...) <- fitMicro

[R] Shading in prediction intervals

I have the following code for the minimum and maximum of my prediction interval > y.down=lines(x[x.order], set1.pred[,2][x.order], col=109) > y.up=lines(x[x.order], set1.pred[,3][x.order], col=109) domain=min(x):max(x) polygon(c(domain,rev(domain)),c(y.up,rev(y.down)),col=109) It doesnt seem t

[R] How do you learn/teach R?

us and fill out the survey on http://www.Rcademy.org. It takes less than 5 minutes. Thank you, Jonathan Cornelissen Doctoral researcher, KU Leuven [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch

[R] problem plotting in a grid

Hi all, I'm trying to generate a grid of four plots. The first 2 appear just fine, but the final 2 will not appear in the grid, instead overwriting the first two.Any ideas on how to get them all in the same window would be greatly appreciated. Cheers, Jonathan library(fields) par(

[R] fitting particular distributions

e alternative processes in R? Any help will be greatly appreciated. Jonathan [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do rea

[R] help with filled.contour() -

Dear all, I can't figure out a way to have more than one plot using filled.contour() in a single plate. I tried to use layout() or par(), but the way filled.contour() is written seems to override those commands. Any suggestions would be really appreciated. Jon

  1   2   3   4   5   6   >