Re: [R] Quantile

2014-09-15 Thread David Winsemius
This is what the code shows: > quantile function (x, ...) UseMethod("quantile") That may be informative in your case. Please read the Posting Guide. It asks that you not crosspost. If you post a followup to rhelp, then the reading of the Posting guide will tell you that much

Re: [R] plots on log="y" scale with smooths

2014-09-17 Thread David Winsemius
ry(ggplot2) ggplot(CPpos, aes(x=length, y=intensity)) + geom_jitter(position=position_jitter(height=.1), alpha=0.25) + scale_y_log10(breaks=c(1,2,5,10,20,50,100, 200)) + stat_smooth(method="loess", color="red", size=1.5) + stat_smooth(method="lm&q

Re: [R] Failure with .Rprofile on Mac OS X

2014-09-18 Thread David Winsemius
Info() and startup settings), over on the R-MAC-SIG mailing list. On Sep 18, 2014, at 11:18 AM, Gang Chen wrote: When R starts in GUI (e.g., /Applications/R.app/Contents/MacOS/R) on my Mac OS X 10.7.5, the startup configuration in .Rprofile works fine. However, when R starts on the terminal (e

Re: [R] R/Ubuntu, “package ‘stats’ in options(”defaultPackages“) was not found”

2014-09-18 Thread David Winsemius
] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] graphics grDevices utils datasets methods base loaded via a namespace (and not attached): [1] tcltk_3.1.1 tools_3.1.1 Any idea? Thanks! -- Davide 2014-09-16 17:15 GMT-04:00 David Winsemius : On Sep 16, 2014,

Re: [R] Failure with .Rprofile on Mac OS X

2014-09-19 Thread David Winsemius
vironment variables present in one launch and absent from > another, but what I was really looking for is whether something in his shell > is changing a path. Because mac environment variables are funky that way.) >> David Winsemius September 19, 2014 at >> 12:57 AM >&g

Re: [R] Converting xy plot arguments into a point argument and fixed x-axis scaling

2014-09-20 Thread David Winsemius
ylab="Value [m]") > axis(1, at = seq(1000, 2050, 5), cex.axis=1, labels=FALSE, tcl=-0.3) > axis(2, at = seq(-10, 10, 500), cex.axis=1, labels=FALSE, > tcl=-0.3) > x <- x[,-1] > segments(x[,2],x[,3],x[,4],x[,5],lwd=2) > dev.off() > } > > __

Re: [R] A new environment within the main function

2014-09-21 Thread David Winsemius
" as well. Thanks and Regards Billy [[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 comm

Re: [R] How to combine character month and year columns into one column

2014-09-23 Thread David Winsemius
t;, "NA-1999", "NA-1999", "NA-1999", >> "NA-1999", "NA-1999")), .Names = c("month", "Year", "views", >> "MonthDay"), row.names = 109:120, class = "data.frame") >>> >> > >

Re: [R] error in rownames

2014-09-25 Thread David Winsemius
233.18, 666.481, 457.721, 1553.81, >>>> 679.505, 1115.53, 515.162, 692.974, 498.604, 204.075, 388.138, >>>> 885.474, 343.097), elevation = c(1901.69, 1992.82, 1911.9, 1985.14, >>>> 1979.67, 1870.83, 1909.5, 2111.45, 1913.09, 1922.76, 1996.68, >>>> 209

Re: [R] modify function in a package

2014-09-25 Thread David Winsemius
ce the NanoStringNorm function but still able to call other > function in the package? > I still have package NanoStringNorm in my search path but somehow it can't > find the other function. Take a look at assignInNamespace. It's also p

Re: [R] Help with conversion of variable labels (Hmisc and Epicalc)

2014-09-26 Thread David Winsemius
dt to use the attributes function to display hte names and structures of the "labels" in the two packages. Firts look at help(pack=Hmisc) and help(pack=Epicalc) to see how the package-specific functions are documented. -- David Winsemius Alameda, CA, USA _

Re: [R] Help with conversion of variable labels (Hmisc and Epicalc)

2014-09-26 Thread David Winsemius
coln in seq_along(.data) ) { label(.data[[coln]]) <- attr(.data, "var.labels")[coln]} > str(.data) 'data.frame': 8 obs. of 2 variables: $ sbp:Classes 'labelled', 'numeric' atomic [1:8] 120 100 110 120 140 120 NA NA .. ..- attr(*, "label"

Re: [R] c() and dates

2014-10-03 Thread David Winsemius
Xlt" objects converts them to the current time zone, and on "POSIXct" objects drops any "tzone" attributes (even if they are all marked with the same time zone)." > Could someone illuminate? > Agree it's "unexpected behavior" and worthy of a featur

Re: [R] c() and dates

2014-10-03 Thread David Winsemius
DT" "1969-12-31 19:04:41 PST" [4] "1969-12-31 19:04:42 PST" c.POSIXct <- function (..., recursive = FALSE) .POSIXct(c(unlist(lapply(list(...), as.POSIXct > c(temp1, temp2) [1] "2000-05-04" "2000-05-05" "2000-05-04" "2000-05-0

Re: [R] Removing description from lm()

2014-10-06 Thread David Winsemius
t at the console or output to a file: xcf <- cat(coef(lm(c(1:3, 7, 6) ~ x))[1]) -0.7 > xcf NULL So `unname` would be preferable if you were trying to store that value. -- David > Thanks! > Billy > > > > On Mon, Oct 6, 2014 at 8:42 AM, David Winsemius > wrote: > ?co

Re: [R] aggregating a zoo object (monthly data to quarterly data)

2014-10-06 Thread David Winsemius
> " 8.750", " 7.060", " 8.590", " 5.490", " 5.660", "10.680", "13.030", > "13.670", "11.620", " 9.160", "11.170", "11.680", "10.040", " 6.670", >

Re: [R] lattice add a fit

2014-10-08 Thread David Winsemius
put a fit from lm (but not a > simple straight line) and the fit has to be done for each panel separately, > not one fit for the full data set, so sth like an lm equivalent of > panel.locfit (there is no panel.lmfit) > Thank you. > > Stephen B > > > [[alter

Re: [R] Obtain a list of data tables from a main data table

2014-10-08 Thread David Winsemius
gt; > Thank > you very much to all the members! > > > [[alternative HTML version deleted]] This is a plain text mailing list. > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman

Re: [R] jitter function when length(x) = 1

2014-10-09 Thread David Winsemius
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. David Winsemius Alameda, CA, USA

Re: [R] Windows RT

2014-10-09 Thread David Winsemius
from contamination by the outside world. http://en.wikipedia.org/wiki/Windows_RT#Differences_from_Windows_8 http://en.wikipedia.org/wiki/Windows_RT#Software_compatibility Welcome to the Hotel Microsoft. -- David. > > Duncan Murdoch > >> >> Ivania Andrade >> > David Wi

Re: [R] understanding the no-label concept

2014-10-11 Thread David Winsemius
>> need it to. >> >> [code] >>> table(x) >> >> 1 2 3 4 >> 17 6 6 2 >> [/code] >> >> I understand that R itself is powerful and well developed. So I try to >> understand why it doesn't support labe

Re: [R] rbind in array of lists

2014-10-12 Thread David Winsemius
gt; __ >> 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, rep

Re: [R] how to overwrite a Unary operator ?

2014-10-17 Thread David Winsemius
Rolf Turner -- Rolf Turner Technical Editor ANZJS __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-co

Re: [R] how to overwrite a Unary operator ?

2014-10-17 Thread David Winsemius
thod dispatch. I was able to emulate that example to create a C-like, unary `+` operator for a new class, but I'm not willing to put it in print for fear that my karmic account might be depleted. -- David. > -- > > PO SU > mail: desolato...@163.com > Majored in Statisti

Re: [R] [survMisc]: error message in examples of comp()

2014-10-17 Thread David Winsemius
ortunes_1.5-2 > > loaded via a namespace (and not attached): > [1] colorspace_1.2-4 digest_0.6.4 gam_1.09.1 grid_3.1.1 > [5] gtable_0.1.2 km.ci_0.5-2 KMsurv_0.1-5 MASS_7.3-33 > [9] munsell_0.4.2plyr_1.8.1 proto_0.3-10 Rcpp_0.11.2 >

Re: [R] Comparing values of different columns: Error: level sets of factors are different

2014-10-17 Thread David Winsemius
ASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- David Winsemius Alameda, CA, USA __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo

Re: [R] Smoothing Data-dplyr

2014-10-18 Thread David Winsemius
e(svel = prod/Z) %>% > ungroup() > > The code works but takes 1 hour. I think the delay is caused by > do(data.frame(kv = .$imDv:.$ipDv)). Is there any faster way to do this? > > [[alternative HTML version deleted]] > >

Re: [R] Build a package on one Mac but unable to load this on another Mac

2014-10-18 Thread David Winsemius
ist. > __ > 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

Re: [R] pasteFromExcel

2014-10-20 Thread David Winsemius
ocale (might > "timezone" suffice?), I'd appreciate that help too. > ?Sys.getlocale # perhaps "LC_MONETARY" ?options # look for OutDec > [[alternative HTML version deleted]] David Winsemius Alameda, CA, USA __ R-help@r-pro

Re: [R] matching genes to a list of gene groups an built binary data frame

2014-10-20 Thread David Winsemius
p20 1 0 0 1 0 > group30 01 1 0 0 > > sapply(Genes, function(x) as.numeric(sapply(ListGroup, '%in%', x=x) ) ) ACACA BAX BCL2 BID BAX MAPK9 [1,] 1 00

Re: [R] making a plot

2014-10-20 Thread David Winsemius
=c(min(ano), max(ano+ranges))) >> segments(ano,1:6,ano+ranges,1:6,pch=25,lty=1,lend=4) >> >> I need to put an asterix (failure) by the end of the three first lines and >> a small circle (censoring) >> by the end of the last three. >> >> Someone can help me? >> >&g

Re: [R] pasteFromExcel

2014-10-20 Thread David Winsemius
" ) > gsub(rmchar, "", c("$1,000", "1,200", "800")) [1] "1000" "1200" "800" > Thank you, David. > -Dan > > On Mon, Oct 20, 2014 at 10:04 AM, David Winsemius > wrote: >> >> On Oct 19, 2014,

Re: [R] format negative numbers

2014-10-20 Thread David Winsemius
for me) murky regular expression waters I >> try to ask community about changing format of negative numbers. >>> >>> For some reason I get a file with negative numbers formatted with >> negative sign at end of number. >>> >>> something like >>> >>>

Re: [R] pasteFromExcel

2014-10-20 Thread David Winsemius
quot;, "$1,000", "(124)", "$(123)", "($123)", > " 1,000 ", "NA") > > # grep will tell you whether elements of x, trimmed of > beginning/ending whitespace, match the currencypattern > grep(currencypattern, trim(x)) >

Re: [R] pasteFromExcel

2014-10-21 Thread David Winsemius
*(\\.\\d{0,2})?|[1-9]{1}\\d{0,}(\\.\\d{0,2})?|0(\\.\\d{0,2})?|(\\.\\d{1,2}))$|^\\-?\\$?([1-9]{1}\\d{0,2}(\\,\\d{3})*(\\.\\d{0,2})?|[1-9]{1}\\d{0,}(\\.\\d{0,2})?|0(\\.\\d{0,2})?|(\\.\\d{1,2}))$|^\\$?\\(([1-9]{1}\\d{0,2}(\\,\\d{3})*(\\.\\d{0,2})?|[1-9]{1}\\d{0,}(\\.\\d{0,2})?|0(\\.\\d{0,2})?|(\\.\\d{1,2}))\\)

Re: [R] assignment 2 part 1,2,3

2014-10-22 Thread David Winsemius
.R-project.org/posting-guide.html ^^^ > and provide commented, minimal, self-contained, reproducible code. ^^^^^^^^ David Winsemius Alameda, CA, USA __ R-h

Re: [R] Fwd: as.POSIXlt() Function

2014-10-22 Thread David Winsemius
ime())); >> [1] "2014-10-21 20:59:23 PDT" >>> >>> as.POSIXlt(Sys.time()); >> [1] "2014-10-21 20:59:23 PDT" >>> c(as.POSIXlt(Sys.time())); >> [1] "2014-10-21 20:59:23 PDT" >>>

Re: [R] big datasets for R

2014-10-27 Thread David Winsemius
elp > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD Alameda, CA, USA __ R-help@r-project.org mailing list https://stat.ethz.ch/m

Re: [R] Order of boxs in box plots

2014-10-27 Thread David Winsemius
quot;, > c("Up","Down"), fill=terrain.colors(2), horiz=TRUE, cex = 0.7,) > title(main="Lower Mudstone Ammonia Box Plots") > [[alternative HTML version deleted]] > You should learn to post in plain text. > ______ >

Re: [R] Variable selection from given data

2014-10-30 Thread David Winsemius
u learn to post in plain text. > and provide commented, minimal, self-contained, reproducible code. … and provide an example to get tested code. David Winsemius, MD Alameda, CA, USA __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listi

Re: [R] problem with posting question: This post has NOT been accepted by the mailing list yet

2014-11-01 Thread David Winsemius
g guide. > and provide commented, minimal, self-contained, reproducible code. David Winsemius Alameda, CA, USA __ 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

Re: [R] treating missing values in timeSeries package

2014-11-02 Thread David Winsemius
, IIT KGP > KGP > __ > 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, re

Re: [R] split a field in dependence of the semicolon

2014-11-03 Thread David Winsemius
>ghi >> 2 2 ; abc; abc >> 3 3 ;def; def >> >> Maybe anyone can help me. Thank you. >> >> Best regards. >> >> Mat >> >> > > _

Re: [R] Variance of multiple non-contiguous time periods?

2014-11-04 Thread David Winsemius
; __ >> snipped >> > > Importing into R wasn't an issue; some of the fields contain spaces & > symbols, but all the fields are tab separated so I can simply use; > > foo <- read.csv("bar",header=T,

Re: [R] Variance of multiple non-contiguous time periods?

2014-11-04 Thread David Winsemius
On Nov 4, 2014, at 9:16 AM, CJ Davies wrote: > On 04/11/14 17:02, David Winsemius wrote: >> >> On Nov 4, 2014, at 8:35 AM, CJ Davies wrote: >> >>> On 04/11/14 16:13, PIKAL Petr wrote: >>>> Hi >>>> >>>>> -Original Message

Re: [R] Variance of multiple non-contiguous time periods?

2014-11-04 Thread David Winsemius
On Nov 4, 2014, at 3:41 PM, CJ Davies wrote: > On 04/11/14 17:42, David Winsemius wrote: >> On Nov 4, 2014, at 9:16 AM, CJ Davies wrote: >> >>> On 04/11/14 17:02, David Winsemius wrote: >>>> On Nov 4, 2014, at 8:35 AM, CJ Davies wrote: >>>> >

Re: [R] upgrade from 3.1.1 > 3.1.2

2014-11-10 Thread David Winsemius
your code, then why did you not provide the code > why can't 'second' run...? thx in advance... R told you. You gave it too many arguments. I seriously doubt this has anything to do with the upgrade. This is not a behavior that changed in R 3.1.2 > -r -- David Winsemiu

Re: [R] Two dimensional likelihood surface plot

2014-11-11 Thread David Winsemius
f Guelph > Guelph, ON > > [[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 pr

Re: [R] Two dimensional likelihood surface plot

2014-11-11 Thread David Winsemius
df, contour( x=unique(V1), y=unique(V2), z= matrix( V3, length(unique(V1)), length(unique(V2) ) ) )) > Gyanendra Pokharel > University of Guelph > Guelph, ON > > On Tue, Nov 11, 2014 at 3:53 PM, David Winsemius > wrote: > > On No

Re: [R] how to use the rpart model to predict new data frame?

2014-11-11 Thread David Winsemius
de http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD Alameda, CA, USA __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Help with ddply/summarize

2014-11-13 Thread David Winsemius
"M", "F"), size = 29, replace = TRUE), age = runif(n = 29, min = 18, max = 54) ) ddply(dfx, .(group, sex), summarize, mean = round(mean( get('age') ), 2)) group sex mean 1 A F 34.81 2 A M 33.38 3 B F 35.89 4 B M 33.67 5

Re: [R] Data Import to R

2014-11-13 Thread David Winsemius
the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. David Winsemius Alameda, CA, USA __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-

Re: [R] "R CMD Rd2txt" generate "_^H" for in all section titles

2014-11-13 Thread David Winsemius
ted. > > Thanks, > Da > > __ > 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

Re: [R] Bootstrap CIs for weighted means of paired differences

2014-11-14 Thread David Winsemius
95% (-0.3276, 0.1594 ) (-0.4781, -0.3477 ) * > > weighted.mean(x=dataset[,1]-dataset[,2], w=dataset[,3]) > > *[1] -0.07321734* > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list &g

Re: [R] boot strapping poisson getting warnings and negative values

2014-11-14 Thread David Winsemius
tinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mail

Re: [R] how to compure R-squared in glm

2014-11-14 Thread David Winsemius
text mailing list > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help Please: > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code

Re: [R] Bootstrap CIs for weighted means of paired differences

2014-11-15 Thread David Winsemius
On Nov 14, 2014, at 3:18 PM, David Winsemius wrote: > > On Nov 14, 2014, at 12:15 PM, ivan wrote: > >> Hi, >> >> I am trying to compute bootstrap confidence intervals for weighted means of >> paired differences with the boot package. Unfortunately, the weig

Re: [R] how to compure R-squared in glm

2014-11-15 Thread David Winsemius
rch/jbes96preprint.pdf — David. > > Thanks > Aravindhan > > -----Original Message- > From: David Winsemius [mailto:dwinsem...@comcast.net] > Sent: Saturday, November 15, 2014 1:16 PM > To: Aravindhan, K > Cc: R-help@r-project.org > Subject: Re: [R] how to compure R

Re: [R] Newbie question: ROC function in TTR package

2014-11-16 Thread David Winsemius
the discrete and the continuous types?? Thanks. It's rather simple to look at the code. Just type: TTR::ROC > wizardchef Reading the code is considered an important first step in answering such questions. -- David Winsemius Alameda, CA, USA

Re: [R] Rose Diagrams for Geology

2014-11-18 Thread David Winsemius
oject.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. David Winsemius Alameda, CA, USA ___

Re: [R] R 3.1.2 (x64) Programming Assignment 1: Air Pollution(corr) do not understand the task

2014-11-18 Thread David Winsemius
ase, tell me where Im wrong. > > Thank you, > Eva > > Best Regards, > Evgeniia-Liza Sheliukhina > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http

Re: [R] Bootstrap CIs for weighted means of paired differences

2014-11-19 Thread David Winsemius
t all. -- David. > > I guess this is rather a statistical question and hence I will have to look > further into it. > > In any case, thanks a lot for your help. > > Best > > > > On 15 Nov 2014, at 17:27, David Winsemius wrote: > >> >> On Nov

Re: [R] Symbolic equations to R code?

2014-11-19 Thread David Winsemius
t; [39] "RcmdrPlugin.SurvivalT" "RcmdrPlugin.NMBU" [41] "RcmdrPlugin.pointG" "RcmdrPlugin.MAc" [43] "RcmdrPlugin.MAd""RcmdrPlugin.mosaic" [45] "RcmdrPlugin.Export"

Re: [R] Bootstrap CIs for weighted means of paired differences

2014-11-20 Thread David Winsemius
set.seed(.) before the runs. > > Level PercentileBCa > 95% (-0.6714, 0.4661 ) (-0.6747, 0.4559 ) > Calculations and Intervals on Original Scale > > On 19 Nov 2014, at 17:49, David Winsemius wrote: > >>>> vw_m_diff <- function(da

Re: [R] [R-pkgs] rms 4.1-0

2014-11-20 Thread David Winsemius
> pas de méthode pour 'calibrate' applicable pour un objet de classe "c('cph', > 'rms', 'coxph')" > > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org

Re: [R] SVYPLOT

2014-11-20 Thread David Winsemius
gt; >>> >>> >> >>[[alternative HTML version deleted]] >> >> __ >> R-help@r-project.org mailing list >> https://stat.ethz.ch/mailman/listinfo/r-help >> PLEASE do read the posting gui

Re: [R] SVYPLOT

2014-11-20 Thread David Winsemius
On Nov 20, 2014, at 12:56 PM, David Winsemius wrote: > > On Nov 20, 2014, at 8:10 AM, Anthony Damico wrote: > >> survey:::svyplot.default with style="grayhex" calls >> hexbin:::gplot.hexbin >> >> an internet search turns up lots of people asking

Re: [R] Bootstrap CIs for weighted means of paired differences

2014-11-21 Thread David Winsemius
S. I think you should consult the code, first. And you should also look at the `stype` parameter where "w" is one option. -- David. > > On Thu, Nov 20, 2014 at 8:19 PM, David Winsemius > wrote: > > On Nov 20, 2014, at 2:23 AM, i.petzev wrote: > > > Hi David

Re: [R] Help with Time

2014-11-21 Thread David Winsemius
outputs 22:30:00. > > if Sys.time is 12:13:22 then I would like to get 12:00:00 etc. There already are round and trunc functions for POSIXt objects so why not take a look a that code and hack it into a shape that you find useful. -- David Winsemius Alameda, CA, USA

Re: [R] Help with Time

2014-11-21 Thread David Winsemius
ke a function that outputs 22:30:00. >> >> if Sys.time is 12:13:22 then I would like to get 12:00:00 etc. >> >> Any help would be appreciated. >> >> Many thanks and regards, >> Raghu > > __ > R-help@r-project.org maili

Re: [R] Help with Time

2014-11-21 Thread David Winsemius
On Nov 21, 2014, at 3:19 PM, David Winsemius wrote: > > On Nov 21, 2014, at 2:55 PM, Raghuraman Ramachandran wrote: > >> Sorry I forgot to mention it clearly. I like to round it to the >> nearest 30th minute that is past. So 12:28:59 will be again 12:00:00 >> and &

Re: [R] Generate random numbers under constrain

2014-11-22 Thread David Winsemius
t. Please respond to the mailing list if appropriate. For those >> needing to send personal or professional e-mail, please use appropriate >> addresses. >> >> >> FREE ONLINE PHOTOSHARING - Share your photos

Re: [R] R 3.1.2 GUI 1.65 Mavericks build (6833) freezes in Yosemite

2014-11-23 Thread David Winsemius
ailman/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. David Winsemius, MD Alameda, CA, USA __ R-help@r-project.org mailing list ht

Re: [R] How to print labels and how to add missing label values?

2014-11-24 Thread David Winsemius
me. > not.png>__ > 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, repr

Re: [R] More elegant way of stacking the data

2014-11-24 Thread David Winsemius
2 c2 7 2 b 3 c2 8 3 c 4 c2 9 4 d 5 c2 10 5 e 6 c2 11 1 a 3 c3 12 2 b 4 c3 13 3 c 5 c3 14 4 d 6 c3 15 5 e 7 c3 16 1 a 4 c4 17 2 b 5 c4 18 3 c 6 c4 19 4 d 7 c4 20 5 e 8 c4 -- David Winsemius Alameda, CA, USA ___

Re: [R] plot.hclust point to older version

2014-11-26 Thread David Winsemius
__ >>> R-help@r-project.org<mailto: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] Survival Analysis and Predict time-to-death

2015-08-18 Thread David Winsemius
It depends on several factors. You need answers to all these questions: How many events occurred, ... and was the period of observation long enough to cover a significant fraction of the life expectancy, …. and is there external evidence or theory that will help establish that this process shou

Re: [R] Bayesian data analysis recommendations

2015-08-19 Thread David Winsemius
If you need the link to the latest email address of a package maintainer just use the maintainer function: ?maintainer > On Aug 15, 2015, at 10:59 AM, Jeff Slagle wrote: > > I have a question about AtelieR out on stackoverflow.com. Perhaps you could > forward the link to: > > Yvonnick Noel

Re: [R] Newbie question: error message with install.packages

2015-08-19 Thread David Winsemius
> On Aug 19, 2015, at 9:13 PM, Peter Wicher wrote: > > Hi, > > I’m starting to work my way through “Machine Learning With R” by Brett Lantz. > > > Running on Mac OS X 10.10.4 > > I’ve downloaded and installed R and the R Console comes up fine. > > Whenever I use the install.packages com

Re: [R] glm help

2015-08-20 Thread David Winsemius
> On Aug 19, 2015, at 8:54 AM, Joaquín Aldabe wrote: > > Dear All, I´m running a glm with poisson errors and have a doubt when > ploting the predicted values. One of my variables has a positive slope in > the summary output, but when I plot the predicted values on the original > plot it draws a

Re: [R] glm help

2015-08-21 Thread David Winsemius
tinfo/r-help >>> PLEASE do read the posting guide >> http://www.R-project.org/posting-guide.html >>> and provide commented, minimal, self-contained, reproducible code. >> > > > > -- > *Joaquín Aldabe* > > *Grupo Biodiversidad, Ambiente y Sociedad*

Re: [R] Generalised poisson regression

2015-08-21 Thread David Winsemius
n plain text rather than HTML. -- David Winsemius Alameda, CA, USA __ 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.ht

Re: [R] Generalised poisson regression

2015-08-22 Thread David Winsemius
p a page for "Undocumented and Internally Used Functions and Classes". Thus endeth today's meanderings among the help pages for pkg:VGAM. Please read the Posting Guide and also go back and read the general information webpages. You seem to have missed some important details. -- D

Re: [R] Error while running swirl package in R

2015-08-23 Thread David Winsemius
nstalled. You did not include the parameter dependencies=TRUE when you installed swirl. > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help &

Re: [R] Calculate the area under a curve

2015-08-24 Thread David Winsemius
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 commented, minimal, self-contained, reproducible code. > >

Re: [R] Choosing columns by number

2015-08-25 Thread David Winsemius
articulated. > > I feel like there must be a better way to do this so I wanted to ask > the collective wisdom here what people do to accomplish this. > Obviously this is a trivial example, but the issue really becomes > problematic when you have a large dataframe. > > Thanks i

Re: [R] how to generate 2 dimensional mesh and hyper-cubes ?

2015-08-25 Thread David Winsemius
quot;spargraphs". I did find some genetics oriented tutorials for students that mentioned producing "spar graphs" with the limma-package. You might address this question to the BioConductor mailing list if your area of interest is genetics/

Re: [R] Choosing columns by number

2015-08-25 Thread David Winsemius
>> ** >> >> Let's not spend our time and resources thinking about things that are so >> little or so large that all they really do for us is puff us up and make us >> feel like gods. We are mammals, and ha

Re: [R] xyplot colour points and layout

2015-08-26 Thread David Winsemius
l=culr,strip=strip.custom(bg='white')) Try putting in a groups argument (after adding the group identifier to the dataframe: xyplot(Abun~Date1|Station, groups=culr, col=c("grey","black"), data=Raw, type="p", xlab=list("Month",cex=1.5), yla

Re: [R] Change the maximum likelihood of multinomial logic model in R

2015-08-28 Thread David Winsemius
Dear all, >> >> Can anyone help me to change the maximum likelihood of multinomial logic >> model in R? >> >> Thanks > > ______ > David Winsemius Alameda, CA, USA __ R-hel

Re: [R] Multiple Integrals

2015-08-29 Thread David Winsemius
this case dlnorm, that is positive definite versus an expectation simply of the sum of such values? -- David Winsemius Alameda, CA, USA __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Fw: Multiple Integrals

2015-08-30 Thread David Winsemius
1+X2+X3)/3-X4) and then take the absolute, then, we will get a lower bound > of the expectation I want to find. > I understood the error in my thinking when Jeff Newmiller pointed out the minus sign that I had missed. Thanks; David. > On Saturday, August 29, 2015 7:24 PM

Re: [R] Multiple Integrals

2015-08-30 Thread David Winsemius
1+X2+X3)/3-X4) and then take the absolute, then, we will get a lower bound > of the expectation I want to find. > > On Saturday, August 29, 2015 7:24 PM, David Winsemius > wrote: Using the adaptIntegrate function in package:cubature I seem to be getting convergence near 5.35

Re: [R] modify strip labels with given text using lattice package

2015-08-31 Thread David Winsemius
"light grey"), >superpose.polygon=list(col= c("yellow", "blue")) >), >scales = list( >alternating = FALSE >), > key = list( >space="top", >columns=2, >text=list(c("Negative", &q

Re: [R] Conditional replacement and removal of data frame values

2015-08-31 Thread David Winsemius
rus 11 2.6 Adenovirus 12 2.7 Rotavirus 02 2.8 Norovirus 12 2.10 Sapovirus 02 3.11 Adenovirus 03 3.12 Rotavirus 03 3.13 Norovirus 03 3.15 Sapovirus 03 > Thank you. > Best regards, > Luigi > >

Re: [R] Median on second group of CSV file produce Na

2015-09-01 Thread David Winsemius
instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml> >> >>>> >>> >>> >>> >>> >>> -- >>> View this message in context: >> http://r.789695.n4.nabble.com/Median-on-second-group-of-CSV-file-p

Re: [R] Urgent query

2015-09-02 Thread David Winsemius
ut it throws an error because the is not a syntactically correct R expression, so I switched to using text versions of the symbols. -- David Winsemius Alameda, CA, USA __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https:/

Re: [R] Use of contrasts in ANOVA

2015-09-02 Thread David Winsemius
84) > Fax: 1-866-PLessThan (1-866-753-7784) > www.PLessThan.com > > __ > 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 htt

Re: [R] simulation in vector autoregressive model (VAR)

2015-09-04 Thread David Winsemius
lp@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 commented, minimal, self-contained, reproduc

Re: [R] Reformatting text inside a data frame

2015-09-07 Thread David Winsemius
> On Sep 7, 2015, at 1:20 PM, Jon BR wrote: > > Hi John, > Thanks for the reply; I'm pasting here the output from dput, with a > 'df <-' added in front: > > df <- structure(list(rowNum = c(1, 2, 3), first = structure(c(NA, 1L, > 2L), .Label = c("AD=2;BA=8", "AD=9;BA=1"), class = "factor"),

<    5   6   7   8   9   10   11   12   13   14   >