[R] New User Having Trouble Loading R Commander on Mac OS Yosemite

2015-04-30 Thread Aaron
I keep getting the same error message when trying to install R Commander. My operating system is Mac OS Yosemite 10.10 I have installed R 3.2, Rstudio, XQuartz (X11), and tcltk-8.x.x-x11.dmg. But I keep getting the following error:Loading required package: splinesLoading required package: RcmdrMis

Re: [R] Help with big data and parallel computing: 500, 000 x 4 linear models

2016-08-08 Thread Aaron Mackey
Don't run 500K separate models. Use the limma package to fit one model that can learn the variance parameters jointly. Run it on your laptop. And don't use %methylation as your Y variable, use logit(percent), i.e. the Beta value. -Aaron On Mon, Aug 8, 2016 at 2:49 PM, Ellis, Alici

[R] Unexpected errors in sparse Matrix arithmetic with zero-length dimensions

2019-02-10 Thread Aaron Lun
ero- dimension inputs. This ensures that my functions return sensible results or fail gracefully in edge cases that might be encountered by users. The current behaviour of sparse Matrix arithmetic causes my unit tests to fail for no (obvious) good reason. Best, Aaron Lun Research Associ

[R] lmer syntax, matrix of (grouped) covariates?

2008-08-18 Thread Aaron Mackey
(cov1 + cov2 + ... | class1) + (...) + (... + cov210 + cov211 | class8) But I'd like keep things syntactically simpler, and use the covariates and cov2class variables directly. I haven't been able to find the right syntactic sugar to get this done. Thanks for any help, -Aaron

[R] hex2RGB back to hex not the same?

2008-08-28 Thread Aaron Mackey
a call to rgb() where I can also set alpha transparency levels ... -Aaron [[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-proje

[R] hmm.discnp or other?

2009-08-12 Thread Aaron Mackey
the trained HMM, I then wish to label new sequences via viterbi(). Am I missing the purpose of this package? I also read through the msm package docs, but my data doesn't really have a time coordinate on which the data should be "aligned". Thanks for any pointers, -Aaron

[R] importing S4 methods using a namespace

2010-02-23 Thread Aaron Rendahl
(lme4, "summary") export("ss") --- But when I call my new function, I get the summary.default method instead of the mer method. > m <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy) > ss(m) Length Class Mode 1mer S4 Thanks, -- Aaron Rendahl,

Re: [R] importing S4 methods using a namespace

2010-02-24 Thread Aaron Rendahl
Thanks very much! Importing from Matrix as you suggest fixes it. -- Aaron Rendahl, Ph.D. Statistical Consulting Manager School of Statistics, University of Minnesota NEW OFFICE (as of June 2009): 48C McNeal Hall, St. Paul Campus 612-625-1062 www.stat.umn.edu/consulting

Re: [R] How to read a list into R??

2009-06-30 Thread Aaron McDaid
x27;s simply a convenient way to log what you see in the terminal. Aaron On Mon, Jun 29, 2009 at 23:58, Li,Hua wrote: > Dear R helpers: >      I have tried many times to find some way to read a list into R. But I > faid. Here is an example: >      I have a file 'List.txt'

[R] Specify CRAN repository from command line

2009-07-26 Thread Aaron Hicks
mand line option for this? Currently I'm doing it using an R script, but that's causing issues because it's not 'visible' to the installer. This would greatly streamline R installation with a standard package set. Regards, Aaron Hicks Please consider the environment befo

[R] Error using getBM() to query BioMart archives

2009-06-16 Thread Aaron Wolen
: type Valid are: name, description, options, fullDescription The same error is returned if I check to see what value type is required for a particular filter: > filterType("chromosome_name", mart=mart) Error in listFilters(mart, what = "type") : The function argument 'wh

[R] deleting/removing previous warning message in loop

2009-03-27 Thread aaron wells
ing("Resetting warning message") This problem has been addressed in a previous listserve string, however I do not follow the advice given. See the below web link. Any help would be greatly appreciated. Thanks! Aaron Wells https://stat.ethz.ch/pipermail/r-help/2008-Octobe

Re: [R] deleting/removing previous warning message in loop

2009-03-27 Thread aaron wells
William, The function keepWarnings that you wrote did the trick. Thanks for the help! Aaron > Subject: Re: [R] deleting/removing previous warning message in loop > Date: Fri, 27 Mar 2009 13:33:51 -0700 > From: wdun...@tibco.com > To: awell...@hotmail.com &g

[R] Public R servers?

2009-04-01 Thread Aaron Barzilai
identifed Wessa.net (http://www.wessa.net/mrc.wasp?outtype=Browser%20Blue%20-%20Charts%20White), but it's not clear it will meet my needs. Can anyone suggest any other resources? Thanks in advance, Aaron Barzilai [[alternative HTML

[R] function output with for loop and if statement

2009-04-22 Thread aaron wells
nd example output. Please help! Thanks, Aaron concov.test<-function(vegetation,specieslist) { test.veg<-vegetation names(test.veg)<-specieslist$LifeForm tmp<-matrix(nrow=nrow(test.veg),ncol=length(unique(names(test.veg for (i in unique(names(

Re: [R] function output with for loop and if statement

2009-04-22 Thread aaron wells
Mark, thanks for the suggestions. Unfortunately that did not fix the problem. I have experimented (with no success) with placing braces in different locations around the if/else statements and removing them all together. Thanks again, Aaron Date: Wed, 22 Apr

Re: [R] function output with for loop and if statement

2009-04-23 Thread aaron wells
, the if statement is an attempt to include the column of data from the dummy.vegdata into the output when there is only one species in a given lifeform. Examples of this condition in the dummy.vegdata include water (Bare_Ground) and popbal (Deciduous_Tree). Aaron > dummy.vegd

[R] rmysql query help

2009-05-06 Thread Aaron Sims
ble in the database and not the value of the variable. If I put it outside the quotes, then it complains about the space. Are there special escape characters or something else Im missing? This date format in a mysql table is pretty standard Any ideas? Thanks, Aaron require(RMySQL) startdatetime

Re: [R] optimization challenge

2010-01-13 Thread Aaron Mackey
n, and references therein: www.pnas.org/content/105/9/3209.abstract -Aaron P.S. Good to see you here Albyn -- I enjoyed your stat. methods course at Reed back in 1993, which started me down a somewhat windy road to statistical genomics! -- Aaron J. Mackey, PhD Assistant Professor Center f

[R] prop.test CI depends on null hypothesis?

2019-10-21 Thread Aaron Rendahl
t;, which was answered only with the quote from the help page: https://stat.ethz.ch/pipermail/r-help/2013-March/350386.html I also found several questions asking which Newcombe method is implemented, which didn't elicit specific answers; here's one from 2011: https://stat.ethz.ch/piper

Re: [R] prop.test CI depends on null hypothesis?

2019-10-22 Thread Aaron Rendahl
avior out there somewhere that I've missed? * if this really is a "new" discovery, how best to bring it to the attention of those who can decide what to do about it? Thanks! On Mon, Oct 21, 2019 at 11:33 AM Aaron Rendahl wrote: > Why does prop.test use continuity correction "o

Re: [R] XYZ data

2013-06-26 Thread Aaron Mackey
for plotting purposes, I typically jitter() the x's and y's to see the otherwise overlapping data points -Aaron On Wed, Jun 26, 2013 at 12:29 PM, Shane Carey wrote: > Nope, neither work. :-( > > > On Wed, Jun 26, 2013 at 5:16 PM, Clint Bowman wrote: > > > Jo

[R] points3d and ordirgl

2012-12-07 Thread Aaron Wells
er 1 to the nmds plot in the color red and size 3. Anyone have an ideas? Thanks, Aaron [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the postin

Re: [R] Renaming variables

2013-09-20 Thread Aaron Mackey
On Fri, Sep 20, 2013 at 10:10 AM, Preetam Pal wrote: > I have 25 variables in the data file (name: score), i.e. X1,X2,.,X25. > > I dont want to use score$X1, score$X2 everytime I use these variables. > attach(score) plot(X1, X2) # etc. etc. -Aaron [[alternative

Re: [R] combine glmnet and coxph (and survfit) with strata()

2013-12-09 Thread Aaron Mackey
I'm also curious how to use glmnet with survfit -- specifically, for use with interval regression (which, under the hood, is implemented using survfit). Can you show how you converted your Surv object formula to a design matrix for use with glmnet? Thanks, -Aaron On Sun, Dec 8, 2013 at

Re: [R] Generate random percentages and placing vectors

2010-10-28 Thread Aaron Lee
ill a room when its empty? Does the room, > the thing itself have purpose? Or do we, what's the word... imbue it." > - Jubal Early, Firefly > > > From: Aaron Lee To: r-help@r-project.org Date: > 10/27/2010 > 11:06 AM Subject: [R] Generate random percentages a

[R] Arrange elements on a matrix according to rowSums + short 'apply' Q

2010-12-02 Thread Aaron Polhamus
X/sum(X)) But my the custom function is still getting applied over the columns of the array, rather than the rows. For a check try: rowSums(TMAT) colSums(TMAT) Row sums here should equal 1... Many thanks in advance, Aaron [[alternative HTML version deleted]] _

Re: [R] Arrange elements on a matrix according to rowSums + short 'apply' Q

2010-12-02 Thread Aaron Polhamus
Ivan and Michael, Many thanks for the tips, those solved my queries. Still interested in how to force custom functions to work over rows rather than columns when using apply, but the MAT/rowSums(MAT) technique is definitely the most efficient way to go for this application. Cheers, Aaron 2010

[R] Writing out data from a list

2011-01-13 Thread Aaron Lee
t; "Random Text" 0 0 0.03 0.01 "Event 2" "Random Text" 0 0 0 0 0.01 etc. Is there some way to continually write text out using a loop and also attaching a string before each data segment? Thank you in advance! -Aaron

[R] Difficult with round() function

2011-01-17 Thread Aaron Polhamus
Even though my benchmark calculation seems to be returning a clean integers to serve as inputs for the creation of the 'positions' variable, for whatever reason R doesn't read it that way. I would be very grateful for any advice on how I can either alter my approach entirely

Re: [R] Shrink file size of pdf graphics

2011-05-19 Thread Aaron Mackey
You can try something like this, at the command line: gs -sDEVICE=pdfwrite -dCompatibilityLevel=1.5 -dPDFSETTINGS=/screen -dNOPAUSE -dQUIET -dBATCH -sOutputFile=output.pdf input.pdf evidently, the new compactPDF() function in R 2.13 does something very similar. -Aaron On Thu, May 19, 2011 at

[R] Variable in file name png

2011-06-08 Thread Aaron Coutino
Hi, I'm having trouble with getting the png function to properly produce multiple graphs. RIght now I have: for (z in data) { png(file=z,bg="white") thisdf<-data[[z]] plot(thisdf$rc,thisdf$psi) dev.off() } Which should take the "data" object, a list of data sets and produce a graph of eac

Re: [R] Running gplots package with Windows 7

2010-05-17 Thread Aaron Gusdon
Thank-you for the replies. I believe I figured out what the problem was. When I installed the package on linux it ran smoothly, but I just need to install a lot of accessory packages to make gplots work with Windows. Thanks again, Aaron 2010/5/17 Uwe Ligges > Additionally, please g

Re: [R] Complex sampling?

2011-03-09 Thread Aaron Mackey
What I think you need is something along the lines of: matrix(c(sample(3:7), sample(3:7), sample(3:7), sample(3:7), ...), nrow=2) now, each column are your random pairs. -Aaron On Wed, Mar 9, 2011 at 1:01 PM, Hosack, Michael wrote: > > -Original Message- > > From: r-help-

[R] Easy 'apply' question

2011-03-10 Thread Aaron Polhamus
what type of data that row/column contains. Here I want a simple mean if the column is numeric, and the first unique value if the column is a character column. As you can see, 'apply' treats all columns as characters the way I've written his function. Any thoughts? Many th

Re: [R] Easy 'apply' question

2011-03-11 Thread Aaron Polhamus
Perfect, thanks Josh! Cheers, A 2011/3/10 Joshua Wiley > Dear Aaron, > > The problem is not with your function, but using apply(). Look at the > "Details" section of ?apply You will see that if the data is not an > array or matrix, apply will coerce it to one (or tr

Re: [R] Reporting odds ratios or risk ratios from GLM

2011-03-15 Thread Aaron Mackey
OR <- exp(coef(GLM.2)[-1]) OR.ci <- exp(confint(GLM.2)[-1,]) -Aaron On Tue, Mar 15, 2011 at 1:25 PM, lafadnes wrote: > I am a new R user (am using it through the Rcmdr package) and have > struggled > to find out how to report OR and RR directly when running GLM models (not &

[R] Determining frequency and period of a wave

2012-06-01 Thread Aaron Patterson
m/tenderlove/rsausage Any help would be greatly appreciated! -- Aaron Patterson http://tenderlovemaking.com/ pgp8lz0FZdvoQ.pgp Description: PGP signature __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

Re: [R] Hardy Weinberg

2011-06-22 Thread Aaron Mackey
H-W only gives you the expected frequency of AA, AB, and BB genotypes (i.e. a 1x3 table): minor <- runif(1, 0.05, 0.25) major <- 1-minor AA <- minor^2 AB <- 2*minor*major BB <- major^2 df <- cbind(AA, AB, BB) -Aaron On Tue, Jun 21, 2011 at 9:30 PM, Jim Silverton wrote: &

Re: [R] Very slow optim()

2011-07-13 Thread aaron Greenberg
Why use a hammer when you need a wrench? Admb seems to be the best tool for the job. It has several slick interfaces with R. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-

Re: [R] HWEBayes, swapping the homozygotes genotype frequencies

2011-10-09 Thread Aaron Mackey
e (this is a common assumption: that AA is the reference, and thus the major, more frequent, allele). -Aaron On Sat, Oct 8, 2011 at 7:52 PM, stat999 wrote: > I evaluated the Bayes factor in the k=2 allele case with a "triangular" > prior under the null as in the example in th

Re: [R] How to selectively sum rows [Beginner question]

2011-10-24 Thread Aaron Siirila
778 2010 10 12419 5995ICN Seoul, South Korea778 LAXLos Angeles, CA 91 2010 Thanks, Aaron -Original Message- From: jim holtman [mailto:jholt...@gmail.com] Sent: Monday, October 24, 2011 11:58 AM To: asindc Cc: r-help@r-project.org Subject: Re: [R] How to sel

Re: [R] How to selectively sum rows [Beginner question]

2011-10-24 Thread Aaron Siirila
The count() function in the plyr package works beautifully. Thanks to Jim, Rainer and Dennis for your help. Best. -Original Message- From: Dennis Murphy [mailto:djmu...@gmail.com] Sent: Monday, October 24, 2011 12:05 PM To: asindc Cc: r-help@r-project.org Subject: Re: [R] How to selecti

Re: [R] Dependency-aware scripting tools for R

2012-04-19 Thread Aaron Mackey
execute non-PBS jobs. -- Aaron J. Mackey, PhD Assistant Professor Center for Public Health Genomics University of Virginia amac...@virginia.edu http://www.cphg.virginia.edu/mackey On Thu, Apr 19, 2012 at 3:27 PM, Sean Davis wrote: > There are numerous tools like scons, make, ruffus, ant, rake,

[R] MLE Estimation of Gamma Distribution Parameters for data with 'zeros'

2008-06-11 Thread Fox, Aaron
soring' algorithm can allow one to use MLE methods to estimate the gamma distribution for data with zero values (Wilkes, 1990, Journal of Climate). I have not, however, found R code to implement this, and, frankly, am not smart enough to do it myself... :-) Any suggestions? Has anyone else run up

[R] Types in grouped multi-panel (lattice) xyplot

2008-04-10 Thread Aaron Arvey
are not right for what I want to do. If anyone knows how to get points and lines in the same panel for the two different groups (probes/segments), I would love to hear about it. If you further know how to use the 'segment' plot in panels for the segments, I would really love to hear ab

Re: [R] Types in grouped multi-panel (lattice) xyplot

2008-04-10 Thread Aaron Arvey
On Thu, 10 Apr 2008, Deepayan Sarkar wrote: > On 4/10/08, Deepayan Sarkar <[EMAIL PROTECTED]> wrote: >> On 4/10/08, Aaron Arvey <[EMAIL PROTECTED]> wrote: >> > Apologetic prologue: I've looked through the mailing list for an answer to >> > this (sinc

[R] variable names when using S3 methods

2008-04-28 Thread Aaron Rendahl
the proper documentation, if it does exist, would be appreciated. Thanks, Aaron Rendahl University of Minnesota School of Statistics # set up two function that both use method "foo" but with different variable names fooA<-function(model,...) UseMethod("

[R] R crash using rpanel on mac os x

2007-11-27 Thread Aaron Robotham
up to date, and I can find no evidence of similar issues from searching the mailing lists. Any suggestions would be appreciated. Aaron __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http:/

[R] Help with lm and multiple linear regression?

2007-12-27 Thread Aaron Barzilai
ix of coefficients, I've only included two for illustration. Additionally, I'd actually like to weight the observations. How would I go about doing that? I also have that as a separate column vector. Thanks, Aaron Here's my session: > margin margin 166.67 2 -58.

[R] Help with lm and multiple linear regression? (Plain Text version)

2007-12-27 Thread Aaron Barzilai
rames. My actual problem will use a much wider matrix of coefficients, I've only included two for illustration. Additionally, I'd actually like to weight the observations. How would I go about doing that? I also have that as a separate column vector. Thanks, Aaron Here's

Re: [R] Help with lm and multiple linear regression? (Plain Text version)

2007-12-27 Thread Aaron Barzilai
se lm. Thanks, Aaron - Original Message From: Tim Calkins <[EMAIL PROTECTED]> To: Aaron Barzilai <[EMAIL PROTECTED]> Cc: r-help@r-project.org Sent: Thursday, December 27, 2007 6:55:57 PM Subject: Re: [R] Help with lm and multiple linear regression? (Plain Text version) consider m

Re: [R] Gumbell distribution - minimum case

2008-09-09 Thread Aaron Mackey
If you mean you want an EVD with a fat left tail (instead of a fat right tail), then can;t you just multiply all the values by -1 to "reverse" the distribution? A new location parameter could then shift the distribution wherever you want along the number line ... -Aaron On Mon, Sep 8,

Re: [R] database table merging tips with R

2008-09-11 Thread Aaron Mackey
I would load your set of userid's into a temporary table in oracle, then join that table with the rest of your SQL query to get only the matching rows out. -Aaron On Thu, Sep 11, 2008 at 2:33 PM, Avram Aelony <[EMAIL PROTECTED]> wrote: > > Dear R list, > > What is the

Re: [R] database table merging tips with R

2008-09-11 Thread Aaron Mackey
Sorry, I see now you want to avoid this, but you did ask what was the "best way to efficiently ...", and the temp. table solution certainly matches your description. What's wrong with using a temporary table? -Aaron On Thu, Sep 11, 2008 at 3:05 PM, Aaron Mackey <[EMAIL PROT

Re: [R] database table merging tips with R

2008-09-11 Thread Aaron Mackey
I guess I'd do it something like this: dbGetQuery(con, "CREATE TEMPORARY TABLE foo ( etc etc)") sapply(@userids, function (x) { dbGetQuery(con, paste("INSERT INTO foo (userid) VALUES (", x, ")")) }) then later: dbGetQuery(con, "DROP TABLE foo"); -

[R] XML package help

2009-01-23 Thread Skewes,Aaron
cteristics/[...@ename='one']", xmlAttrs) and it does not work. Might somebody help me with the syntax here? Thanks a lot!! Aaron [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/l

Re: [R] XML package help

2009-01-26 Thread Skewes,Aaron
Thanks! Works like a charm. -Aaron From: Duncan Temple Lang [dun...@wald.ucdavis.edu] Sent: Friday, January 23, 2009 6:48 PM To: Skewes,Aaron Cc: r-help@r-project.org Subject: Re: [R] XML package help Skewes,Aaron wrote: > Please consider this: >

[R] XML package- accessing nodes based on attributes

2009-02-09 Thread Skewes,Aaron
File> My requirement is to access eValues at each node based on FileTypeId. For example: How can I get the eValue of eName="PatientReference" for all Type="Patient" ,where the ? i.e. "TCGA-06-0875-01A" and "TCGA-06-0875-02A" For the life

[R] an S idiom for ordering matrix by columns?

2009-02-19 Thread Aaron Mackey
11122 2 1 2 639221 2 1 2 756221 2 1 2 Thanks for any suggestions! -Aaron [[alternative HTML version deleted]] __ R-help@r-

Re: [R] an S idiom for ordering matrix by columns?

2009-02-19 Thread Aaron Mackey
Thanks to all, "do.call(order, as.data.frame(y))" was the idiom I was missing! -Aaron On Thu, Feb 19, 2009 at 11:52 AM, Gustaf Rydevik wrote: > On Thu, Feb 19, 2009 at 5:40 PM, Aaron Mackey wrote: > > There's got to be a better way to use order() on a matrix than this:

[R] Help with R and MySQL

2009-02-23 Thread Aaron Barzilai
upgrading the packages, and upgrading R from 2.5.0 to 2.8.1. However, I still get the same errors. Has anyone run into this problem before? Any suggestions on how to solve it? Thanks in advance, Aaron [[alternative HTML version deleted]] _

Re: [R] Help with R and MySQL

2009-02-23 Thread Aaron Barzilai
Thanks Jeff, that was exactly the problem. When I unzipped the version at the page below for my version of MySQL (5.1), it worked fine. The version I downloaded through install.packages() must have been for 5.0. Thanks so much for the help and quick response, Aaron

[R] ordering

2009-03-10 Thread aaron wells
numeric(test[,2]) [1] FALSE > is.double(test[,2]) [1] FALSE > is.integer(test[,2]) [1] FALSE > is.real(test[,2]) [1] FALSE My questions are: Why is this happening? and How do I fix it? Thanks in advance!

Re: [R] ordering

2009-03-10 Thread aaron wells
Thanks Peter, that did the trick. I'll modify my function so that the numeric conversion is done automatically thus saving me the extra step of converting later on. Aaron Wells > Subject: RE: [R] ordering > Date: Wed, 11 Mar 2009 08:41:50 +1300 &g

[R] geometric mean of probability density functions

2009-03-18 Thread Aaron Spivak
there a way to set the x-values to be the same for multiple density plots? Also, there are no negative values in the dataset, so I'd like to bound the x-axis at 0 if at all possible? Is there a standard way to combine density functions? Thanks for the advice. -Aaron Spivak ps. I thoug

[R] Automatic command line options

2008-03-18 Thread Aaron MacNeil
Hello users, Does anyone know how to turn off the automatic double quoting and bracketing on the command line that appeared in R 2.6.x (OS X). It's driving me nuts! Many thanks, Aaron. ---- M. Aaron Ma

[R] read.table not clearing colClasses

2008-11-11 Thread Skewes,Aaron
is read in second. It appears that colClasses fails to adopt the new structure after the first iteration. Is there some way to clear colClasses of flush the write buffer between iterations? Thanks, Aaron for(i in 1:length(fullnames.in)) { cnames<- read.table(fullnames

[R] 2D density tophat

2008-11-29 Thread Aaron Robotham
le to find either don't support a simple tophat, or don't use the absolute distances between the sampling grid and the data. Should anyone know of more general 2D density routines that might support circular tophats, or know of a simple and efficient method for creating them, I would be v

[R] 2D density tophat

2008-12-01 Thread Aaron Robotham
od for creating them, I would be very grateful. Thanks for your time, Aaron PS: I tried sending this on Friday originally, but as far as I know that didn't work, so should another post appear from me asking the same thing I apologise in advance. [[alterna

[R] making sense of posterior statistics in the deal package

2008-12-04 Thread Aaron Tarone
eatly appreciated. Aaron Tarone __ Aaron Tarone Postdoctoral Research Associate Molecular and Computational Biology Program University of Southern California [EMAIL PROTECTED] (213) 740-3063 __ R-help@r-project.org mailing

Re: [R] 2D density tophat

2008-12-05 Thread Aaron Robotham
numeric(names(temp))]=as.numeric(temp) z <- matrix(pad, length(gx), length(gy)) if(density){z=z/(nx*pi*h^2)} return=list(x = gx, y = gy, z = z) } It works in more or less the same way as kde2d, but by default it returns counts not densities Aaron On 1 Dec 2008, at 11:46, Aaron Robotham wro

[R] very long integers

2008-12-11 Thread Aaron Robotham
internal storage up to LONG, or is the only solution to read it in as a character from the original data? In case anyone is curious, I didn't create the IDs, and in some form I must conserve all of the ID information for later use. Thanks,

[R] very long integers

2008-12-11 Thread Aaron Robotham
internal storage up to LONG, or is the only solution to read it in as a character from the original data? In case anyone is curious, I didn't create the IDs, and in some form I must conserve all of the ID information for later use. Thanks,

[R] Coxph frailty model counting process error X matrix deemed singular

2009-06-24 Thread Porter, Aaron T
Hello, I am currently trying to simulate data and analyze it using the frailty option in the coxph function. I am working with recurrent event data, using counting process notation. Occasionally, (about 1 in every 100 simulations) I get the following warning: Error in coxph(Surv(start, end,

[R] which rows are duplicates?

2009-03-29 Thread Aaron M. Swoboda
ate, but again it doesn't tell me which row it was a duplicate of (as far as I can tell). Am I missing something? How can I determine that row 3 is a duplicate OF ROW 1? Thanks, Aaron __ R-help@r-project.org mailing list https://stat.ethz.ch/m

[R] confusion over "names" of lm.influence()$hat

2009-04-14 Thread Aaron M. Swoboda
88 observations, but they don't appear to be the same observations. Thanks for anyone who can help me clear this up, Aaron __ 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] Save a graph file use jpeg(file=file)

2022-01-05 Thread Ebert,Timothy Aaron
The simple solution for Windows is to use (windows icon) + shift + s. You then select a portion of your screen and it gets copied to your clipboard. You can then paste that into your document. Of course this will not work if it is important that the reader is able to rotate the graphic. Tim --

Re: [R] NAs are removed

2022-01-14 Thread Ebert,Timothy Aaron
Hi Neha, You used a variable named "fraction" so we took a guess. However, as another pointed out 1/0 does not give NA in R. number/0 returns Inf except 0/0 which returns NaN. So 1/0 <= 1 returns FALSE and 0/0 <= 1 returns NA. A great deal of the behavior of your program hinges on what "fract

Re: [R] NAs are removed

2022-01-14 Thread Ebert,Timothy Aaron
y, January 14, 2022 5:11 PM To: Ebert,Timothy Aaron Cc: Jim Lemon ; r-help mailing list Subject: Re: [R] NAs are removed [External Email] I have a variable in dataset "CA", which has the following values: [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

Re: [R] [External] Weird behaviour of order() when having multiple ties

2022-01-31 Thread Ebert,Timothy Aaron
Dat1 <- c(0.6, 0.5, 0.3, 0.2, 0.1, 0.1, 0.2) print(order(Dat1)) print(sort(Dat1)) Compare output -Original Message- From: R-help On Behalf Of Martin Maechler Sent: Monday, January 31, 2022 9:04 AM To: Stefan Fleck Cc: r-help@r-project.org Subject: Re: [R] [External] Weird behaviour of

Re: [R] Convert a character string to variable names

2022-02-08 Thread Ebert,Timothy Aaron
will keep that in mind. I agree that it makes life much easier to enter the data in the way it will be analyzed. -Original Message- From: Jeff Newmiller Sent: Tuesday, February 8, 2022 10:10 PM To: r-help@r-project.org; Ebert,Timothy Aaron ; Richard O'Keefe ; Erin Hodgess Cc: r-

Re: [R] Convert a character string to variable names

2022-02-12 Thread Ebert,Timothy Aaron
How does “a value” differ from “an object?” From: Richard O'Keefe Sent: Friday, February 11, 2022 12:25 AM To: Ebert,Timothy Aaron Cc: Jeff Newmiller ; r-help@r-project.org; Erin Hodgess Subject: Re: [R] Convert a character string to variable names [External Email] You wrote "32

Re: [R] confusion matrix like detail with continuous data?

2022-02-16 Thread Ebert,Timothy Aaron
In your prediction you will have a target level of accuracy. Something like "I need to predict the slope of the regression to within 1%." You break your data into a training and testing data sets, then for the testing data set you ask is the prediction within 1% of the observed value. That is ab

Re: [R] Problem with data distribution

2022-02-17 Thread Ebert,Timothy Aaron
You pipe the filter but do not save the result. A reproducible example might help. Tim -Original Message- From: R-help On Behalf Of Neha gupta Sent: Thursday, February 17, 2022 1:55 PM To: r-help mailing list Subject: [R] Problem with data distribution [External Email] Hello everyone

Re: [R] Problem with data distribution

2022-02-17 Thread Ebert,Timothy Aaron
; > > > > > > This can be copied into an R session and the data set recreated > > > with > > > > > > data <- structure(etc) > > > > > > > > > Now the boxplots. > > > > > > (Why would you want to plot a vec

Re: [R] conditional filling of data.frame - improve code

2022-03-10 Thread Ebert,Timothy Aaron
You could try some of the "join" commands from dplyr. https://dplyr.tidyverse.org/reference/mutate-joins.html https://statisticsglobe.com/r-dplyr-join-inner-left-right-full-semi-anti Regards, Tim -Original Message- From: R-help On Behalf Of Jeff Newmiller Sent: Thursday, March 10, 2022 1

Re: [R] How important is set.seed

2022-03-21 Thread Ebert,Timothy Aaron
If you are using the program for data analysis then set.seed() is not necessary unless you are developing a reproducible example. In a standard analysis it is mostly counter-productive because one should then ask if your presented results are an artifact of a specific seed that you selected to g

Re: [R] How important is set.seed

2022-03-22 Thread Ebert,Timothy Aaron
:03 PM To: Ebert,Timothy Aaron Cc: Jeff Newmiller ; r-help@r-project.org Subject: Re: How important is set.seed [External Email] Thank you again Tim d=readARFF("my data") set.seed(123) tr <- d[index, ] ts <- d[-index, ] ctrl <- trainControl(method = "repeatedc

Re: [R] How important is set.seed

2022-03-22 Thread Ebert,Timothy Aaron
difference. Random.org might be another place for the OP to explore. Tim -Original Message- From: Bert Gunter Sent: Tuesday, March 22, 2022 12:12 PM To: Neha gupta Cc: Ebert,Timothy Aaron ; r-help@r-project.org Subject: Re: [R] How important is set.seed [External Email] OK, I'm somewhat

Re: [R] How important is set.seed

2022-03-22 Thread Ebert,Timothy Aaron
way to address this is to characterize the distribution of model outcomes from different runs with different random seeds. Examine that characterization and hope for understanding. Tim From: Bert Gunter Sent: Tuesday, March 22, 2022 2:03 PM To: Ebert,Timothy Aaron Cc: Neha gupta ; r-help@r

Re: [R] What is the intended behavior, when subsetting using brackets [ ], when the subset criterion has NA's?

2022-04-06 Thread Ebert,Timothy Aaron
I get an error with this: my_subset_criteria <- c( F, F, T, NA, NA) my_subset_criteria Tim -Original Message- From: R-help On Behalf Of Kelly Thompson Sent: Wednesday, April 6, 2022 4:13 PM To: r-help@r-project.org Subject: [R] What is the intended behavior, when subsetting using bracke

Re: [R] Error with text analysis data

2022-04-13 Thread Ebert,Timothy Aaron
Cc: Ebert,Timothy Aaron ; r-help mailing list Subject: Re: Error with text analysis data [External Email] Someone just told me that you need to pre process the data before model construction. For instance, make the text to lower case, remove punctuation, symbols etc and tokenize the text (give

Re: [R] Symbol/String comparison in R

2022-04-13 Thread Ebert,Timothy Aaron
https://en.wikipedia.org/wiki/ASCII There is a table towards the end of the document. Some of the other pieces may be of interest and/or relevant. Tim -Original Message- From: R-help On Behalf Of Kristjan Kure Sent: Wednesday, April 13, 2022 10:06 AM To: r-help@r-project.org Subject: [R

Re: [R] R Code Execution taking forever

2022-04-24 Thread Ebert,Timothy Aaron
1) Does it run perfectly with num_tirals_6 <- 100 ? 2) Rework the code to remove as much as possible from loops. Renaming column names each time through the loop seems pointless. Is the nested for loops converting the dice roll to person name necessary within the while loo

Re: [R] Confusing fori or ifelse result in matrix manipulation

2022-04-25 Thread Ebert,Timothy Aaron
A <- matrix(1:9,ncol=3) x <- c(0,1,0) M <- matrix(ncol=3,nrow=3) M<-A for(i in 1:3) { if(x[i]){ M[,i] <-0 } } } M The outcome you want is to set all of the middle column values to zero. So I used x as a logical in an if test and when true everything in that column is set to zero. Y

Re: [R] Convert one rectangle into two rectangles

2022-04-27 Thread Ebert,Timothy Aaron
One thought was that these were farm fields in satellite images. I just had no clue why one would subdivide them. Maybe a training and testing portion of these polygons? At a few dozen meters the curvature might be negligible on a sufficiently large spheroid. If 45-55% split is approximate then

Re: [R] Is there a canonical way to pronounce CRAN?

2022-05-04 Thread Ebert,Timothy Aaron
It would be nice in some ways if everyone would pronounce the same word in the same way, but then we could not argue over the correct pronunciation of words like tomato or aluminium/aluminum. I think of cran as "Kran". While I had German in high school I didnn't remember the German word for cran

Re: [R] R and .asc file extension

2022-05-20 Thread Ebert,Timothy Aaron
A google search returned a stack overflow page that might help. stackoverflow.com/questions/20177581/reading-an-asc-file-into-r (add the https part to get a functional link.) I would also try looking at the file using something like notebook, or any program that is a plain text editor. That way I

Re: [R] Suggestions as to how to proceed would be appreciated...............

2022-05-22 Thread Ebert,Timothy Aaron
Would lm, nls, or nlme work for what you need? Tim -Original Message- From: R-help On Behalf Of Bernard Comcast Sent: Sunday, May 22, 2022 3:01 PM To: Bert Gunter Cc: R-help@r-project.org Subject: Re: [R] Suggestions as to how to proceed would be appreciated... [External E

  1   2   3   >