[R] Unable to get the 64-bit version of R (3.2.1 Terminal or GUI) to start on a Windows 7 machine

I am unable to get the 64-bit version of R (3.2.1 Terminal or GUI) to start on a Windows 7 machine. I can get the 32-bit to start, just not the 64-bit. I am receiving a dialog box box that points me to four files but nothing seems usable when examining them. Installed to C:\Dev\R\R-3.2.1. I hav

[R] Unexpected/undocumented behavior of 'within': dropping variable names that start with '.'

a data frame you don't have to remove all extra variables at the end). I just wish it was documented. Cheers, Brian sessionInfo() R version 3.1.0 (2014-04-10) Platform: x86_64-pc-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_US.UTF-8

[R] Error: not 'a real'?

"numeric", "numeric", "numeric", "numeric") # Error: Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, (from #2) : scan() expected 'a real', got '894.6' How is '894.6

[R] supply methods to read.table

Hi List, I would like to optimize some data reading as well as clean up some code. The manual tells me to supply methods to colClasses but the manual and the methods documentation aren't helping... Can someone provide me an example please? Best, Brian R version 2.15.1 (2012-06-22) Pla

Re: [R] Error: not 'a real'?

Hi Jim, On 10/20/12 13:36, Jim Holtman wrote: how about supplying the context of the error. Show the lines in the file where the error occurred. > > Sent from my iPad > > On Oct 20, 2012, at 7:21, Brian wrote: > >> Hi List, >> >> when supplying a v

Re: [R] Error: not 'a real'?

Hi Milan, that's right, forgot about that one! Sorry to bother. Brian On 10/20/12 14:37, Milan Bouchet-Valat wrote: Le samedi 20 octobre 2012 à 14:25 +0200, Brian a écrit : Hi Jim, On 10/20/12 13:36, Jim Holtman wrote: how about supplying the context of the error. Show the lines in the

[R] apply and sapply oddity

quot;character" "character" "character" "character" "character" usp200s0 fve010z0 dkl010z0chron CO_sync CO2_dry_sync "character" "character" "character" "character" "character" &q

[R] Easy Uplift Tree Classify Error

Does anyone know if the error being generated when trying to predict test set data in the Easy Uplift Tree package is something fixable by the user or is this a bug in the program making the package essentially non-operable? This is from the package documentation and fails on the last step of a

[R] SVM Prediction and Plot

Hi I'm trying to predict using a model I fitted with SVM. I constructed the model (called Svm) using a training set, and now I want to use a test set (called BankTest) for prediction. The response variable is in the first column of BankTest. > SvmPred = predict(Svm, BankTest[,-1], probability=

Re: [R] SVM Prediction and Plot

Thanks for the reply. When I run this: > head(SvmPred) 3570 2361 5406 2041 3440 4123 NNYYYY Levels: N Y -- View this message in context: http://r.789695.n4.nabble.com/SVM-Prediction-and-Plot-tp3254421p3254570.html Sent from the R help mailing list archive at Nabble.co

Re: [R] SVM Prediction and Plot

There are a good few blanks in some variables in the data. There were in the training set too. Is that a problem? I don't have any na.action in my svm call either. So na.omit=na.fail causes this error: Error in na.fail.default(newdata) : missing values in object I tried the matrix. >SvmPred =

[R] warnings about factor levels dropped from predict.glm

version of R (not sure what number hers is) on her Mac, she gets these warnings constantly. I've checked some records manually by doing the algebra and the predict.glm() function is working correctly incorporating the factor levels on my machine. Any thoughts??? Brian Brian S. Cade, PhD

[R] glm$effects

stant quantity except for one of the values. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: ca...@usgs.gov tel: 970 226-9326 [[alternative HTML version de

Re: [R] question regarding the AICcmodavg package

among predictors and some don't. Stick to model averaging the predicted responses and you can do something that is sensible and that can be applied to any combination of predictor variables including those with interactions. Brian Cade Brian S. Cade, PhD U. S. Geological Survey Fort Co

[R] R 3.5 and RTools

I 'd like to confirm. Thanks in advance, Brian __ 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-guide.html and provide

[R] numeric comparison error

Hi, I am a little bit perplexed at why I am getting some values as FALSE: > cpgbins <- seq(0,1,0.05) > cpgbins [1] 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.00 > cpgbins[1] == 0.00 [1] TRUE > cpgbins[2] == 0.05 [1] TRUE > cpgbins[3] =

[R] grep

Hi, I was trying to find a pattern ("ABHD14A") in a character string ('xgen' in example below) using grepl. Note that the individual members may be separated by a semi-colon. The correct answer should return: "ABHD-ACY1 ; ABHD14A" "ABHD14A ; YYY" I have tried three approaches, but still seem a

[R] Extract cells and their adjacent cells that may appear anywhere in a dataframe.

appreciated. Brian [[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 posting guide http://www.R-project.org/posting-guide.html and

[R] Jagged ROC curves?

Hi, I was trying to draw some ROC curves (prediction of case/control status), but seem to be getting a somewhat jagged plot. Can I do something that would 'smooth' it somewhat? Most roc curves seem to have many incremental changes (in x and y directions), but my plot only has 4 or 5 steps even tho

Re: [R] Jagged ROC curves?

i', ylab='', xlab='') plot(smooth(roc_1),col="brown3", lwd=2, add=T, lty=1) I guess most ROCs that I've seen are somewhere in between, i.e. they have a little jaggedness, but not as much as in plot #1 above" thanks! On Mon, Jun 26, 2017 at 12:59 PM, Marc S

[R] ggplot2 geom_bar arrangement

Hi, I was trying to draw a geom_bar plot. However, by default, the bars are arranged according to the label, which I don't want. I want the bars to appear exactly as they appear in the data frame. For example in the code: Lab=c(letters[4:6],letters[1:3]) valuex = c(3.1,2.3,0.4,-0.4,-1.2,-4.4)

[R] ggplot2 geom_bar label justification

Hi, I was trying to make a horizontal bar plot. The barplot works when the text labels are of reasonable length, but not if some of them are slightly long. I think the long ones get 'squeezed' by default before the plot is flipped and keep the skew after the flip. Is there a way I can get around t

Re: [R] ggplot2 geom_bar arrangement

; > # set the factor levels to the same order as observed in the data frame > df$Lab <- factor(df$Lab, levels=unique(df$Lab)) > > px <- ggplot(df,aes(Lab,valuex,label=Lab)) + > geom_text(aes(y=0)) + > geom_bar(stat = "identity") > px > > Jean > >

[R] par(mfrow) for heatmap plots

Hi, I was trying to use par(mfrow) to put 4 heatmaps on a single page. However, I get one plot per page and not one page with 4 plots. What should I modify? Test code is given below: test = matrix(rnorm(60), 20, 3) pdf(file='test.pdf',width=10,height=8) par(mfrow=c(2,2)) heatmap(test) heatmap(te

Re: [R] Power Calculation:2-sided exact equivalence test for Binomial Proportions

MJ: I think the EnvStats package has various power functions for binomial applications (also confidence interval half-widths). Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: ca...@usgs.gov tel: 970

Re: [R] Linear model vs Mixed model

easily compared to the lmer() model estimates. They are close but will never be identical as the lmer() model estimates are based on assuming a normal distribution with specified variance. They rarely would be identical. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science

Re: [R] Linear model vs Mixed model

th an estimated variance to locate the subject specific intercepts. Why do you think those modeled intercepts should exactly coincide with your fixed effects intercepts that makes no such distributional assumption? Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center

Re: [R] glmmLasso with interactions errors

It has never been obvious to me that the lasso approach can handle interactions among predictor variables well at all. I'ld be curious to see what others think and what you learn. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C

[R] transposing x and y axes in xYplot

Is there a simple way to transpose the x and y axes with the xYplot() function in the Hmisc package, where y is a vector of point estimate and lower and upper confidence interval endpoints? What I'm looking for is something akin to coord_flip() used with ggplot(). Brian Brian S. Cade, PhD

Re: [R] transposing x and y axes in xYplot

Bert: Yes, with some fiddling of axes labels this looks like just what I needed. Thank you. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: ca...@usgs.gov tel: 970 226-9326 On Tue, Aug 23, 2016 at

Re: [R] R-square prob is not calculated by randomization in lmPerm::lmp

rcept term the permutation test for Ho: B1 = B2 = ... Bp = 0 (i.e., all coefficients other than the intercept = 0) is equivalent to a permutation test for Ho: R-squared = 0. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO

[R] box type in Hmisc xYplot

Does anyone know how to change the box type in Hmisc package function xYplot. I want only the left and bottom axes drawn, similar to what I would accomplish with bty="l" argument in plot() function. bty= argument did not do anything for me in xYplot(). Brian Brian S. Cade,

Re: [R] box type in Hmisc xYplot

No even in xyplot() the scales argument is not eliminating plotting of the top and right graph axes. It is not the scale of the axes I want eliminated but the actual lines. But thank you for trying. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave

Re: [R] Drawing the regression line and the 95% confidence intervals

The prediction intervals are likely to be much wider than the confidence intervals so you will need to be sure you scale the yaxis limits large enough to see them. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818

[R] plot: rug colors

Hi, I wanted the rug (in plot) to have different colors. For example: vals1 <- sample(1:100,5) vals2 <- sample(1:100,5) rugcols <- c("red","blue","brown","red","yellow") plot(vals1,vals2) rug(vals1,col=rugcols,lwd=2) However, with this code I only get 'red' for all the ticks. Is there a way I

Re: [R] plot: rug colors

Thanks Duncan! That works! On Sun, May 24, 2015 at 8:09 AM, Duncan Murdoch wrote: > On 24/05/2015 7:47 AM, Brian Smith wrote: > > Hi, > > > > I wanted the rug (in plot) to have different colors. For example: > > > > vals1 <- sample(1:100,5) > > vals2 &l

Re: [R] alternatives to KS test applicable to K-samples

samples but only for univariate responses. Both are included in the USGS Blossom package for R linked here: https://www.fort.usgs.gov/products/23735 (not yet distributed via CRAN). The MRPP may also be available in other R packages on CRAN (vegan ?). Brian Brian S. Cade, PhD U. S. Geological Survey

Re: [R] alternatives to KS test applicable to K-samples

Not sure what the issue is here. If I click on the link I provided, I go right to the USGS page where instructions for downloading the software are provided. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818

Re: [R] dplyr help

​dplyr solution: bevs %>% group_by(name, sex, drink) %>% summarise(​cost = sum(cost)) %>% select(name, drink, cost, sex) The last select statement puts the output in the column order you wanted in your result. I hope this helps. Brian On Wed, Jul 29, 2015 at 9:37 PM, Jon BR wrote:

[R] Unmarked

it or do i need a separate unmarkedFrameOccu for each species? Thanks -Brian Brian Ortman Quantitative Ecologist/Avian Specialist ATWELL, LLC 440.349.2000 Tel 334.524.7334 Mobile 440.349.2028 Fax 7100 E. Pleasant Valley Rd. | Suite 220 | Independence, OH 44131 www.atwell-group.com <http

Re: [R] Extracting elements out of list in list in list

On 1/16/15 9:34 AM, Bert Gunter wrote: Chee Hee's approach is both simpler and almost surely more efficient, but I wanted to show another that walks the tree (i.e. the list) directly using recursion at the R level to pull out the desired components. This is in keeping with R's "functional" progra

Re: [R] Percentage cover data with many zeros

of the cumulative distribution function about which you are able to make no inferential statements. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: ca...@usgs.gov tel: 970 226-9326 On Sat, Jan 24

[R] read.table - replaces 'T' with 'TRUE'

Hi, I had a table with 'T's in it. When I try to read the table, it replaces all the 'T's with TRUE. Is there a way that I can retain the Ts? thanks! Sample table (test.txt): xx1 xx2 0 0 2 2 A A G G A G A A C C G G G G G G A G A A A G A G A A G G A A A G G G A G A G G G A G A G G G G G A G G G

Re: [R] read.table - replaces 'T' with 'TRUE'

Thanks!! On Mon, Jan 26, 2015 at 5:23 PM, peter dalgaard wrote: > > > On 26 Jan 2015, at 23:10 , Duncan Murdoch > wrote: > > > > read.table(, colClasses = "character") > > > > (You might want "factor" instead of character.) > > Or maybe not. I'd expect trouble with getting the levels set to > c

[R] Fixed Width EBCDIC Files in R

I'm trying to read some mainframe data encoded as EBCDIC into R, and am at a loss. I'd like to avoid using an external program to convert the files, since I'm operating in a corporate environment. You can find the example files at at the link below, with both ASCII and EBCDIC versions. Note that t

Re: [R] Fixed Width EBCDIC Files in R

McKown wrote: > On Thu, Feb 5, 2015 at 2:08 PM, Brian Trautman > wrote: > >> I'm trying to read some mainframe data encoded as EBCDIC into R, and am at >> a loss. I'd like to avoid using an external program to convert the files, >> since I'm operating in

[R] MApply and SubStr

Hi! I'm trying to write a custom function that applies SubStr to a string, and then depending on the arguments, converts the output to a number. The substring part of my code works fine, but it's not converting the way I want to -- options('stringsAsFactors'=FALSE) require(data.table) substr_ty

Re: [R] MApply and SubStr

> > On Feb 10, 2015, at 3:58 PM, Brian Trautman wrote: > > > Hi! > > > > I'm trying to write a custom function that applies SubStr to a string, > and > > then depending on the arguments, converts the output to a number. > > > > The substring

[R] Noob: How to upgrade from 3.0.2 on Ubuntu 14.04 LTS?

One of my favorite blogs (AnalyzeCore) uses the dplyr library which fails in my environment (R 3.0.2.) I've tried upgrading my R package to get around this issue, to no avail. I know I'm making a simple mistake but haven't solved it. A simple "$ sudo apt-get install r-base --upgrade" didn't get the

[R] Merging Data.Tables on conditions other than equality

I have two tables that I would like to join together in a way equivalent to the following SQL. Note that I'm using a "greater than" statement in my join, rather than checking for equality. require(sqldf) require(data.table) dt <- data.table(num=c(1, 2, 3, 4, 5, 6), char=c('A', 'A', 'A', 'B', 'B

[R] Variable Column Names and Lists

This should be a simple question, but I am at my wits end. dt<-data.table(a=rep(1:10, 26), b=1:260, c=rep(1:2, 130)) sumvar <- 'mysum' bvar <- 'b' dt_min <- dt[, list(sumvar = sum(get(bvar))), by=list(a)] print(dt_min) I want the function to return two variables, "a" and "mysum". However, it i

Re: [R] Variable Column Names and Lists

This worked perfectly, thank you! On Fri, Feb 27, 2015 at 10:47 AM, Duncan Murdoch wrote: > On 27/02/2015 1:18 PM, Brian Trautman wrote: > > This should be a simple question, but I am at my wits end. > > > > dt<-data.table(a=rep(1:10, 26), b=1:260, c=rep(1:2, 130))

[R] Using dates in R

can help me with this. Let me know if you have any recommendations as to how this can be done relatively easily. Thanks! Appreciate it. Best, Brian [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more

[R] Tracing a variable name back

(p1, p2) { check_stuff(p1, p2) } bar <- function(x1, x2) { foo(x1, x2) } > check_stuff(df, "y") [1] "df" [1] "y" > foo(df, "y") [1] "p1" [1] "p2" > bar(df, "y") [1] "p1" [1] &quo

[R] ggplot2 boxplot points (size,color,shape)

Hi, I am trying to create a boxplot (with geom_jitter) such that the points from one set of values are shown as circles, and the second set of points also as circles, but with no fill. In other words, how can I control the shape and color for the points appearing in this boxplot? === library(

[R] hash - extract key values

Hi, I was trying to use hash, but can't seem to get the keys from the hash. According to the hash documentation ('hash' package pdf, the following should work: > hx <- hash( c('a','b','c'), 1:3 ) > class(hx) [1] "hash" attr(,"package") [1] "hash" > hx$a [1] 1 > keys(hx) Error in (function (classe

Re: [R] hash - extract key values

gt; keys(hx) > [1] "a" "b" "c" > > > Maybe restart your session? Clear your workspace? Upgrade? > > B. > > > > > > On Mar 27, 2015, at 7:39 PM, Brian Smith wrote: > > > Hi, > > > > I was trying to use hash, bu

Re: [R] hash - extract key values

ics grDevices utils datasets methods base other attached packages: [1] hash_3.0.1 loaded via a namespace (and not attached): [1] tools_3.1.2 > On Sat, Mar 28, 2015 at 5:02 AM, Uwe Ligges wrote: > Try to reinstall hash. Sounds like a broken installation. > > Best, > Uwe Ligge

Re: [R] hash - extract key values

t; Then there will be 26 items in the hash table and keys(yourhash) > will return the 26 lowercase letters. Is that what you want? > > Bill Dunlap > TIBCO Software > wdunlap tibco.com > > On Sat, Mar 28, 2015 at 2:56 PM, Brian Smith > wrote: > >> I deleted t

[R] ggplot: connect points with line (not in order)

Hi, I am trying to connect points, but not in a different order than the default value in ggplot. For example: xx <- sample(1:100,5) yy <- sample(1:100,5) mydat <- data.frame(xx,yy) print(mydat) ggplot(mydat,aes(xx,yy)) + geom_point() + geom_line() I want to connect the points as the

Re: [R] ggplot: connect points with line (not in order)

ts-in-ggplot2-in-r > > > > -- > | > > http://billyam.com || http://use-r.com || http://shinyserver.com (BETA) > > SAS Certified Base Programmer for SAS 9 > Oracle SQL Expert(11g) > > > > On Sun, Apr 5, 2015 at 9:03 PM, Brian Smith > wrote: > > Hi

Re: [R] Preventing loading of user packages

it were the only one or if it were the path to the base and stats packages. Fo me this reports correctly that my base package is not there: grepl("~/Library/R/3.0/library" , attr(packageDescription('base'), "file") ) [1] FALSE But this raises the question. Why is t

Re: [R] ddply question

eDate,Age),summarise,salvage=sum(TotalSalvage),Count=sum(TotalCount),pct=Count/sum(Count)) [[alternative HTML version deleted]] -- Brian S. Diggs, PhD Senior Research Associate, Department of Surgery Oregon Health & Science University __ R-help@r-project.org mai

[R] How best to model these datasets

Hello, I am a stats novice and I was wondering what kind of model would work best with the two datasets that are displayed. The graphs show an incremental area analysis for two feral cat home ranges where area increases with increasing number of GPS fixes. Thanks, Brian -- Brian Leo School

[R] package multicore: check progress?

Hi, I use multicore package quite a lot. However, I want to find a way to check on the progress of my job. For example: ftest <- function(x){ if(x %% 100 == 0) print(x) y <- 2x return(y) } res <- mclapply(1:1000,ftest) This would print the value of x in a for loop, but doesn't produce anythin

Re: [R] Comparing two populations based on the percentile values calculated from two independent samples

h quantile regression. Ecological Applications 21: 281-289) for some examples. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: ca...@usgs.gov tel: 970 226-9326 On Thu, Sep 24, 2015 at 12:31 PM, Mohamed A.Ab

Re: [R] algorithmic method quantile regression

han for the br (standard Barrodale and Roberts simplex linear program) algorithm. All conditions that could lead to different estimates. My recommendation for the sample sizes you are considering is to stick with the Barrodale and Roberts algorithm as it is the best understood, most reliable procedure.

Re: [R] monte carlo simulations in permanova in vegan package

Sean: There are only 20 possible combinations, 6!/(3! x 3!), so you just need to enumerate them completely (no Monte Carlo approximation required). I don't know if permanova() can do this but you can do it with the mrpp() functions and argument (,exact=TRUE) in Blossom package for R.

[R] ggplot2: Controlling width of line

Hi, I was trying to increase the size of certain lines in my plot (samples 'B' and 'D' in example below). However, when I try to modify the line size, I seem to screw up the linetypes. Also, is there a way to reflect the line size in the legend? Here is some sample code for illustration: library

[R] ggplot2: remove axis ticks

Hi, I was trying to remove the axis tick marks and their values using theme() but haven't had much success. Here is sample code: rx <- sample(1:100,10) ry <- sample(1:100,10) rz <- sample(letters[1:3],10,replace=T) rdf <- data.frame(rx,ry,rz) p <- ggplot(rdf,aes(x=rx,y=ry)) p1 <- p + geom_point(

[R] ChAMP: champ.runCombat error with methylation 450k data

Hi, I was trying to run COMBAT on methylation data, but keep on getting an error: Error in while (change > conv) { : missing value where TRUE/FALSE needed The error occurs irrespective of whether I give the entire or reduced (variation filter keeps only about 140k CpGs) datasets. Is there any o

Re: [R] Right censored data, abundant in zeros for regression analysis.

quantiles associated with the mass of zeros. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: ca...@usgs.gov tel: 970 226-9326 On Thu, Dec 24, 2015 at 6:41 AM, REES T. (706713) < t.rees.706...@swansea.ac

[R] Reducing of car package when loading

to not load all of the datasets? Thanks, Brian [[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 posting guide http://www.R

Re: [R] Reducing of car package when loading

I am not sure what functions I use in car actually. I am dependent on the drc package which imports "car". The main drc function I use is "drm", along with "LL.4" https://cran.r-project.org/web/packages/drc/index.html <https://cran.r-project.org/web/packa

Re: [R] Heteroscedasticity in a percent-cover dataset

, and incorporate predictor variables any way you would in any other linear (or generalized linear) model. There even are mixed-effects versions of quantile regression now (package lqmm) but I haven't used them enough to speak to their veracity and value. Brian Brian S. Cade, PhD U. S. Geolo

Re: [R] lm() with spearman corr option ?

I think you would just need to replace the lm() function call with cor(x,y,method="spearman". It would probably be more informative to actually plot by the magnitude of the correlation coefficient (all |r| >= 0.20 or something similar) rather than just by those with P <=0.05.

[R] x-axis tick marks on log scale plot

Hi, I have a plot with log scale on the axes. How do I add ticks and labels in addition to the ones provided by default? Can I specify where I want the ticks and labels? For example: set.seed(12345) x <- sample(1:1,10) y <- sample(1:1,10) plot(x,y,log="xy") For me, this plot has tick

Re: [R] x-axis tick marks on log scale plot

> -- > https://www.researchgate.net/profile/Ivan_Calandra > https://publons.com/author/705639/ > > > Le 19/05/2016 à 15:40, Brian Smith a écrit : > >> Hi, >> >> I have a plot with log scale on the axes. How do I add ticks and labels in >> addition to th

Re: [R] Confidence interval on quantile regression (quantreg)

, this typically happens for smaller sample sizes, more extreme taus, and more complex models. But unusual distributional characteristics of the data distributions can also contribute to this issue. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg

Re: [R] GAM with the negative binomial distribution: why do predictions no match with original values?

Well part of the issue is that the negative binomial estimates are for means and they can differ a fair bit from the raw counts, but I'm also guessing that part of the issue is that the offset may not be accounted for with the predict.gam() function. Brian Brian S. Cade, PhD U. S. Geolo

Re: [R] non-parametric manova with post-hoc test

uping structure. There are several R packages that I think implement MRPP but the Blossom package might be one of the better implementations in terms of alternatives provided (including permutation version of Hotelling's test). Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Scie

[R] A housekeeping question...

il. Just a pointer in the right direction would be much appreciated and I'll go quietly off and search on my own. Regards... -- Brian Lunergan Nepean, Ontario Canada "They know not the old ways. That shall be their undoing..." --- __ R-h

[R] [Fwd: ***HTML***R help]

g. I dont understand how thats possible? Sayeed: Wish I could help, but I have no idea why this would be happening. I've bounced this into the R Help group where there are more knowledgeable and experienced people who should be able to answer you without any problem -- Brian Lunergan Nepea

[R] gnls not optimizing values

85468073 -0.10985684 0.05285202 0.45446735 1.81225021 Residual standard error: 0.01055068 Degrees of freedom: 19 total; 13 residual Any thoughts as to why it's just spitting back the initial guesses, and how it might be avoided? Thanks. -Brian Doctrow [[alternative HTML ve

Re: [R] Problems with boxplot in ggplot2:qplot

df$t) + theme_bw() + scale_x_date(major = "months", minor = "weeks", format = "%b") > qplot(factor(df$t), df$x, geom = "boxplot") + theme_bw() + > scale_x_date(format = "%b") qplot(df$t, df$x, geom = "boxplot", group=df$t) + t

Re: [R] cut.POSIXt misconception/feature/bug?

os mingw32 > system i386, mingw32 > status Under development (unstable) > major 2 > minor 11.0 > year 2010 > month 03 > day09 > svn rev51229 > lang

Re: [R] cut.POSIXt misconception/feature/bug?

arification would be to say "Using both ‘right = TRUE’ and ‘include.lowest = TRUE’ will include both ends of the range of breaks". That makes it clear that the dates that are referred to there are the breaks, not x. > Thank you. > > Petr I've included a patch agains

Re: [R] cut.POSIXt misconception/feature/bug?

On 3/11/2010 11:52 AM, Brian Diggs wrote: I've included a patch against cut.POSIXt.Rd with these proposed changes. Apparently the patch didn't make it through, so I'm just pasting it below. -- Brian Diggs, Ph.D. Senior Research Associate, Department of Surgery, Oregon H

[R] Cross-validating two matrices

Hello, I am trying to help a colleague with an R problem (see below) to whit I can only generate a very inelegant solution. Any advice would be welcome. Thanks, Brian If you have two matrices, say a species by trait matrix (Matrix 1 below) and a plot by species matrix (Matrix 2 below

Re: [R] character vector -> numeric matrix ??

If your entire dataset looks like your example, try: do.call(rbind, sapply(zz, strsplit, split="")) Note that strsplit() requires a character vector as its first parameter. -Brian J. Koch Data Manager Decision Development Inc www.decisiondevelopment.com -Original Message- Fr

Re: [R] Re ferencing columns and pulling selected data

748 3201942 4413124 5 421789 6462960 7482247 8597957 9612876 10 700470 11 743862 12 940230 13 1011521 14 1269353 > JR > -- Brian Diggs, Ph.D. Senior Rese

[R] gnls not optimizing values

85468073 -0.10985684 0.05285202 0.45446735 1.81225021 Residual standard error: 0.01055068 Degrees of freedom: 19 total; 13 residual Any thoughts as to why it's just spitting back the initial guesses, and how it might be avoided? Thanks. -Brian Doctrow [[alternative HTML ve

Re: [R] [R-pkgs] formatR: farewell to ugly R code... Problem with too old ver. of Gtk

Hardy Heron) with all the called for updates. I've looked in synaptic but I'm not an 'under the hood' type of user so I have no clue what to look for in there to find out what I have or even if it is possible to update. Oh yes, I'm running v2.10.1 of R. Is there any likely s

[R] Follow up on installing formatR...

help *** installing help indices ** building package indices ... * DONE (formatR) The downloaded packages are in ‘/tmp/RtmpNL20Be/downloaded_packages’ Warning message: In install.packages("formatR") : installation of package 'RGtk2' had non-zero exit status > -- Brian Lunerg

Re: [R] Select single column, preserve name?

ase, I want "v4"). I guess I could kludge this by adding a dummy column that always passes and trimming it out of the output, later, but it seems like there ought to be an easier way. Am I misusing the bracket notation, or ignorant of some obvious way to subset just one column and r

Re: [R] coda sweave

-- Ista Zahn Graduate student University of Rochester Department of Clinical and Social Psychology http://yourpsyche.org [[alternative HTML version deleted]] -- Brian Diggs Senior Research Associate, Department of Surgery, Oregon Health & Science University _

Re: [R] packages gdata / gtools - installation in R 2.11.0

Oswaldo Cruz Foundation Rua Leopoldo Bulhoes 1480 - 8A 21041-210 Rio de Janeiro - RJ Brazil E-mail: pbar...@infolink.com.br E-mail: paulo.bar...@ensp.fiocruz.br -- Brian Diggs Senior Research Associate, Department of Surgery, Oregon Health & Science University

Re: [R] Fitting a quasipoisson distribution to univariate data

and computing novice. I start it with 'library(Rcmdr)' It has a quasipoisson model under the menu choice 'Statistics|fit models|Generalized Linear Models' in the family box, way down near the bottom. I am using Rcmdr version 1.4-10 and R version 2.8.1 Brian -- View this m

Re: [R] debug an error that incapacitates R?

ohn Fox on this. If it is an invalid sink, this would just be the first call that generates output. It would be a sink somewhere before that that is setting up a bad connection (which doesn't cause problems until something is to be written). > Thank you > Liviu > -- Brian Diggs, P

Re: [R] debug an error that incapacitates R?

> -Original Message- > From: Rolf Turner [mailto:r.tur...@auckland.ac.nz] > Sent: Monday, December 14, 2009 2:07 PM > To: Brian Diggs > Cc: r-help@r-project.org Help > Subject: Re: [R] debug an error that incapacitates R? > > > Wouldn't sink.number() gi

  1   2   3   4   5   6   7   8   9   10   >