Re: [R] stop a function

2015-01-30 Thread David Winsemius
_ > 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 commented, minimal, self-contained, reprod

Re: [R] Package corpcor: Putting symmetric matrix entries in vector

2015-01-30 Thread David Winsemius
t.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, reproducible code. David Winsemius Alameda, CA, USA __

Re: [R] error in submission

2015-01-31 Thread David Winsemius
"Heart.Failure"] <- > suppressWarnings(as.numeric(minHos[, "Heart.Failure"])) >minHosN <- na.omit(minHos) >minHospital <- minHosN[with(minHosN, order(Heart.Failure, > Hospital.Name)), ] >return(minHospital[[1]][[1]]) > }else if(outcome == "pne

Re: [R] Regression Overdispersion?

2015-02-01 Thread David Winsemius
/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 -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinf

Re: [R] How to enable https for R 3.1.2 on windows 8.1

2015-02-02 Thread David Winsemius
fo = >> "C:/users/john/documents/twitter/cacert.pem") >>>> Error in function (type, msg, asError = TRUE) : >>>> Protocol " https" not supported or disabled in libcurl >>>> >>>> [[alternative HTML version deleted]] &

Re: [R] How to enable https for R 3.1.2 on windows 8.1

2015-02-02 Thread David Winsemius
fo = >> "C:/users/john/documents/twitter/cacert.pem") >>>> Error in function (type, msg, asError = TRUE) : >>>> Protocol " https" not supported or disabled in libcurl >>>> >>>> [[alternative HTML version deleted]] &

Re: [R] rgl::writeWebGL( , prefix = , )

2015-02-03 Thread David Winsemius
at CRAN and it fails there, too. (Sometimes I am able to get packages to compile that report errors on CRAN.) The first error reported from efforts at installing both 1201 and 1208 versions is: checking for X... libraries , headers checking for glEnd in -lGL... no configure: error: missing

Re: [R] Removing 99% similar sequence help

2015-02-04 Thread David Winsemius
p > 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 __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and m

Re: [R] How to download and unzip data in a loop

2015-02-05 Thread David Winsemius
7-23285-%4i.gz', i)) >>>> url = sprintf(' >>>> ftp://ftp.ncdc.noaa.gov/pub/data/noaa/%4i/724927-23285-%4i.gz', i, i) >>>> #file = basename(url) >>>> download.file(url, file[i]) >>>> files = dir(tmpdir, '*.g

Re: [R] Superscript in legend without using expression function

2015-02-07 Thread David Winsemius
; > -- > Rolf Turner > Technical Editor ANZJS > Department of Statistics > University of Auckland > Phone: +64-9-373-7599 ext. 88276 > Home phone: +64-9-480-4619 > > ______ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and m

Re: [R] MApply and SubStr

2015-02-10 Thread David Winsemius
rgument: > str( mapply(substr_typeswitch, startdata, + startpos, endpos, typelist, SIMPLIFY=FALSE) ) List of 2 $ V1: chr [1:2] "aa" "bb" $ NA: num [1:2] 1 2 > > Can anyone help me understand what I'm doing wrong? > > Thank you! > > [[altern

Re: [R] Pass additional arguments to do.call(grid.arrange, plots)

2015-02-15 Thread David Winsemius
NSUBSCRIBE 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. David Winsemius Alameda, CA, USA

Re: [R] Image processing using R

2015-02-17 Thread David Winsemius
; > Javier > > [[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 provide comm

Re: [R] Substituting elements in vector

2015-02-17 Thread David Winsemius
st[c("B","E")] <- replist > unlist(my.vlist) A B1 B2 B3 B4 B5 C D E1 E2 E3 F G "A" "p" "q" "r" "s" "t" "C" "D" "x" "y" "z" "F" "G"

Re: [R] grepping out columns

2015-02-18 Thread David Winsemius
; "THUR" dat[ sapply(sads, grep, dat) ] #[1] "SA_TUES" "CH_TUES" "AR_TUES" "SA_MON" "CH_MON" "AR_MON" #[7] "SA_THUR" "CH_THUR" "AR_THUR -- David Winsemius Alameda, CA, USA

Re: [R] Numerical stability of eigenvalue and hessian matrix in R

2015-02-18 Thread David Winsemius
ian-matrix-in-R-tp4703443p4703456.html >>>> Sent from the R help mailing list archive at Nabble.com. >>>> >>>> __ >>>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >>>> https://stat.et

Re: [R] package 'corrplot'

2015-02-19 Thread David Winsemius
, "Greys"))) `corrplot.mixed`'s help page says any extra arguments get passed to corrplot, so look at the help page for corrplot where teh tl.col parameter is described: corrplot.mixed(cor(dat), col = rev(brewer.pal(10, "Greys")),

Re: [R] Subsetting a list of lists using lapply

2015-02-20 Thread David Winsemius
a == shas ) >>> [1] 194 >>> HTH, >>> Chuck >>> __ >>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >>> https://stat.ethz.ch/mailman/listinfo/r-help >>> PLEASE do read t

Re: [R] Chi-square test

2015-02-20 Thread David Winsemius
p11/sum(p11) )==1 [1] TRUE > But when I want to do the the Chi-square test, I get this error: > probabilities must sum to one. > Does anybody know the reason? Numerical accuracy. See R-FAQ 7.31 -- David Winsemius Alameda, CA, USA __ R-hel

Re: [R] Split a dataframe by rownames and/or colnames

2015-02-20 Thread David Winsemius
1.04 -2.16 $o1 c_1 d_1 e_1 a_p b_p c_p 1_o1 2_o1 3_o1 1_o1 -0.93 -0.02 0.69 -0.67 1.04 1.04 -1.50 -0.36 0.50 2_o1 0.02 -0.16 -0.09 -1.50 -0.02 -1.04 1.07 -0.45 1.56 3_o1 -1.42 0.88 -0.05 0.85 -1.35 0.21 1.35 0.92 -0.76 $p c_1 d_1 e_1 a_p b_p c_p 1_o1

Re: [R] glm.fit when family argument is not a "family" object

2015-02-22 Thread David Winsemius
lish_United States.1252 > > ## attached base packages: > ## [1] stats graphics grDevices utils datasets methods base > > ## loaded via a namespace (and not attached): > ## [1] tools_3.1.2 > > > -- > Mathew W. McLean > Research Assistant Professor >

Re: [R] Simultaneous Probit-Tobit

2015-02-22 Thread David Winsemius
... I reach for the sos-package: install.packages('sos') library(sos) findFn("tobit") # 98 links findFn("tobit probit") found 43 matches; retrieving 3 pages 2 3 Downloaded 25 links in 11 packages. >

Re: [R] Split a dataframe by rownames and/or colnames

2015-02-23 Thread David Winsemius
0.50 >> 2_o1 0.02 -0.16 -0.09 -1.50 -0.02 -1.04 1.07 -0.45 1.56 >> 3_o1 -1.42 0.88 -0.05 0.85 -1.35 0.21 1.35 0.92 -0.76 >> >> $p >> c_1 d_1 e_1 a_p b_p c_p 1_o1 2_o1 3_o1 >> a_p -1.35 0.91 -0.58 -0.63 0.94 -1.13 0.71 0.25 0.82 >&

Re: [R] Environment error in the car package Anova.lme function

2015-02-24 Thread David Winsemius
run outside of a function, it runs cleanly. > mod <- lme(Y~X,random=~1|Subject) > Anova(mod) > > mod1 <- lmer(Y ~ X +(1 | Subject) ) > Anova(mod1) > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list -- T

Re: [R] Saving Mean Relative Difference from all.equal()

2015-02-26 Thread David Winsemius
d that I got this as a result: > all.EQ( seq(1,2, by=0.1)*10, 10:20) [1] TRUE > all.EQ( seq(1,2, by=0.11)*10, 10:20) [1] "Numeric: lengths (10, 11) differ" > all.EQ(seq(1,2, by=0.11)*10, 10:19) [1] 0.03225806 — David Winsemius, MD Alameda, CA, USA __

Re: [R] integrate with vector arguments

2015-02-26 Thread David Winsemius
1)$value > fV<-Vectorize(f1) > outer( X=v1,Y=v2, FUN= fV) [,1] [,2] [,3] [,4] [,5] [1,] 0.833 1.17 1.5 1.83 2.17 [2,] 1.333 1.67 2.0 2.33 2.67 [3,] 1.833 2.17 2.5 2.83 3.17 [4,] 2.333 2.67 3.0 3

Re: [R] integrate with vector arguments

2015-02-27 Thread David Winsemius
nks a lot. > > Cheers, > > Marko > > > > On 02/27/2015 06:49 AM, David Winsemius wrote: >> >>> On Feb 26, 2015, at 1:49 PM, marKo wrote: >>> >>> v1<-c(1:5) >>> v2<-c(1:5) >>> >>> f1<-function (x) {

Re: [R] R 3.1.2 using a custom function in aggregate() function on Windows 7 OS 64bit

2015-03-05 Thread David Winsemius
cbind(). -- David. > > [[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

Re: [R] Dividing a map by electoral results with R?

2015-03-05 Thread David Winsemius
ividing-a-map-by-electoral-results-with-R-tp4704211.html > Sent from the R help mailing list archive at Nabble.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 http://ww

Re: [R] Confidence Intervals for survreg (survival)

2015-03-07 Thread David Winsemius
ant lines which I am now attempting to put before you (again) 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] question about the geometric distribution, please

2015-03-09 Thread David Winsemius
hank you! > > Sincerely, > Erin > > > -- > Erin Hodgess > Associate Professor > Department of Mathematical and Statistics > University of Houston - Downtown > mailto: erinm.hodg...@gmail.com > > [[alternative HTML version deleted]] -- David Wins

Re: [R] Extracting slots from an object (e.g.: object produced by unit root test function "urdfTest")

2015-03-12 Thread David Winsemius
" " > [27] " Critical values for test statistics: " > [28] "1pct 5pct 10pct" > [29] " tau3 -3.96 -3.41 -3.12" > [30] " phi2 6.09 4.68 4.03" > [31] " phi3 8.27 6.25 5.34" > That is being displayed as

Re: [R] write.csv to text string?

2015-03-12 Thread David Winsemius
a manner that breaks the string into useful or readable segments on a "page". If you want to have a text version of what the print.data.frame function returns, then wrap capture.output around print(dfrm). That will produce a character object. --\n David Winsemius Alameda, CA, USA _

Re: [R] Extracting slots from an object (e.g.: object produced by unit root test function "urdfTest")

2015-03-12 Thread David Winsemius
`summary` on the first argument and pull the p-value from the returned object. Many possibilities. -- David. > > An example for ur.df function: > > data(Raotbl3) > attach(Raotbl3) > lc.df <- ur.df(y=lc, lags=3, type='trend') > summary(lc.df) > > 2015-03-13

Re: [R] graphs

2015-03-13 Thread David Winsemius
damage >> caused by modifications of the e-mail or by delay with transfer of the >> email. >> >> In case that this e-mail forms part of business dealings: >> - the sender reserves the right to end negotiations about entering into a >> contract in any time, for an

Re: [R] need help on RATQ package

2015-03-13 Thread David Winsemius
power > of attorney is submitted to the recipient or the person represented by the > recipient, or the existence of such authorization is known to the recipient > of the person represented by the recipient. > __ > R-help@r-project.org

Re: [R] R Help Pops um in Webbrowser | how to dissable

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

Re: [R] Extract month and year as one column

2015-03-13 Thread David Winsemius
6.97181051627531, 5.74282686202818, 4.37018386569785, 12.5725962512824, >4.6583055309578, 8.76457542037641, 10.7070862034423, 12.84023567151, >5.78620621848167, 5.98643374478599, 13.0993210289842)), .Names = > c("date", > "temp"), row.names = c(NA,

Re: [R] Is RcppOctave package available for Mac to install?

2015-03-13 Thread David Winsemius
rew. > > Thanks, > > Mike > > [[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.

Re: [R] Persons with objects –> Objects and persons

2015-03-16 Thread David Winsemius
gt; them. As it is shown on the picture. >> >> <http://r.789695.n4.nabble.com/file/n4704655/Persons%2C_Objects.png> >> >> Can you please help me how to make this transformation in R? >> >> Thank you >> And don't forget that Nabble is not

Re: [R] How to filter data using sets generated by flattening with dcast, when I can't store those sets in a data frame

2015-03-16 Thread David Winsemius
="," ) ) ) >>>return( setequal( union( a, b ), a ) ) >>> } >>> >>> Here are some example calls: >>> >>> > includes2( '1,2,3', c(1) ) >>> [1] TRUE >>> > includes2( '1,2,3', c(1,2)

Re: [R] Returning to parent function

2015-03-16 Thread David Winsemius
like to return from telemStats > } > > telemStats <- function(a,b){ >b <- c(10,12) >continueIfTrue( {length(b) >=10 }, "progStats","00") >print("Since the above condition failed, I would not like this > message

Re: [R] Returning to parent function

2015-03-17 Thread David Winsemius
> > Ideally, calling telemStats() should return to the prompt and the > print in telemStats should not appear > > > On Mon, Mar 16, 2015 at 4:02 PM, David Winsemius > wrote: >> >> On Mar 16, 2015, at 3:08 PM, Saptarshi Guha wrote: >> >>> Hello, &g

Re: [R] remove outliers with scatterplotMatrix

2015-03-19 Thread David Winsemius
//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 __ R-help@r-projec

Re: [R] ordering a boxplot

2015-03-21 Thread David Winsemius
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 -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] UNSUBCRIBE

2015-03-21 Thread David Winsemius
st -- 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. David Winsemius, MD Alameda, CA, USA _

Re: [R] Joining two datasets - recursive procedure?

2015-03-22 Thread David Winsemius
;>>>>>>>> without >>>>>>>>>>> success. >>>>>>>>>>> >>>>>>>>>>> Following your suggestion, I am providing a sample of the >>>> dataset I >>>>>>>>

Re: [R] average wind direction calculation

2015-03-22 Thread David Winsemius
> On Mar 22, 2015, at 9:18 PM, Alemu Tadesse wrote: > > Dear All, > > I am wondering if you have an R script or know an R package for an average > wind direction > pkg:circular > Best, > > Alemu > > [[alternative HTML version deleted]] R is

Re: [R] Fast evaluation of functions in 3D domains

2015-03-23 Thread David Winsemius
__ > 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 commented, minimal, self-contained, reproducible code. David Winsemius Alame

Re: [R] Fast evaluation of functions in 3D domains

2015-03-25 Thread David Winsemius
> of that result divided by the volume? I cannot imagine any faster procedure . >> >> >>> But I need to do it over the whole volume. >>> >>> For 2 dimensions it can be done with >>> >>> A <- as.im(function,window,parameters) >>&g

Re: [R] Fast evaluation of functions in 3D domains

2015-03-25 Thread David Winsemius
On Mar 25, 2015, at 8:45 AM, David Winsemius wrote: > > On Mar 25, 2015, at 4:55 AM, > wrote: > >> Dear all, >> >> Finally I have tried three different options to integrate a function in a 3D >> volume. Here >> I show a test example. The volume

Re: [R] Converting charter string to vector of numbers

2015-03-25 Thread David Winsemius
r-string-to-vector-of-numbers-tp4705097.html > Sent from the R help mailing list archive at Nabble.com. > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE

Re: [R] matrix manipulation question

2015-03-27 Thread David Winsemius
> > > [[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:

Re: [R] From replicate to accesing sublists

2015-04-01 Thread David Winsemius
be wrapping around to later columns. -- David. > > Can you please advice me if there is a better approach to do that in R? > Regards > Alex > > > [[alternative HTML version deleted]] This is a plain text mailing list. Please reconfigure your email client to sent in

Re: [R] From replicate to accesing sublists

2015-04-02 Thread David Winsemius
hat loops in R are inefficient. It is usually the algorithm inside the loop that is inefficient. -- David. > > Regards > Alex > > > > On Wednesday, April 1, 2015 5:54 PM, David Winsemius > wrote: > > > > On Apr 1, 2015, at 3:34 AM, Alaios via R-help wrote:

Re: [R] Calculating Kendall's tau

2015-04-02 Thread David Winsemius
On Apr 2, 2015, at 7:06 AM, S Ellison wrote: >> I am analyzing trend using Mann-kendall test for 31 independent sample, >> each >> sample have 34 years dataset. I supposed to find Kendall “tau” for each >> sample. The data is arranged in column wise (I attached the data).To find >> Kendall

Re: [R] Cut breaks in descending order

2015-04-03 Thread David Winsemius
___ > 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 commented, minimal, self-conta

Re: [R] applying cumsum within groups

2015-04-03 Thread David Winsemius
2 1486 0 0 32 >11 1002 1 0 33 >12 1229 0 0 33 >13 1142 0 0 33 > 54 39 1797 0 0 3 3 > 54 40 1517 0 0 33 > 54 41 1527 0 0 33 > 54 42 1514 0 0 33 > > [[alte

Re: [R] Aggregating daily rainfall raster data to bimontly data

2015-04-05 Thread David Winsemius
[alternative HTML version deleted]] >>>>>>> >>>>>>> __ >>>>>>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >>>>>>> https://stat.ethz.ch/mailman/listinf

Re: [R] sort adjacency matrix

2015-04-06 Thread David Winsemius
. . . . 1 . . . . >> 5 . . 1 . . . . . . 1 >> 6 1 . . 1 . . . . 1 . >> 7 . 1 . . . . . 1 . . >> 8 . . 1 . . . 1 . . 1 >> 9 . . . . . 1 . . . 1 >> 10 . . 1 . 1 . . 1 1 . >> >> I want to sort it according to ones in each row ascending (where max number

Re: [R] Fast multiple match function

2015-04-06 Thread David Winsemius
you not post in HTML format. I have used gmail and I do know that it is fairly easy to post in plain text. -- David Winsemius Alameda, CA, USA __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listi

Re: [R] strip levels

2015-04-06 Thread David Winsemius
s and frustrated now. T_T > > Regards, > Christine > > __ > 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-pro

Re: [R] Change the position of label when using R package eulerr

2018-09-13 Thread David Winsemius
int 4 .. .. .. .. .. .. ..- attr(*, "unit")= chr "native" -- David. > > > Thank you, > > Aimin > > [[alternative HTML version deleted]] > > ______ > R-help@r-project.org mailing list -- To UNSUBSC

Re: [R] Question on Binom.Confint

2018-09-13 Thread David Winsemius
> > [[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://

Re: [R] Help with setting locale

2018-09-14 Thread David Winsemius
stions for correct command. > > With thanks > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do re

Re: [R] Help with setting locale

2018-09-14 Thread David Winsemius
> On Sep 14, 2018, at 8:44 AM, David Winsemius wrote: > > >> On Sep 14, 2018, at 1:02 AM, Kim Titcombe wrote: >> >> *Query or Set Aspects of the Locale* >> >> I have an issue with setting LOCALE in installation (new installation on >> new

Re: [R] New to R

2018-09-14 Thread David Winsemius
t; > >Sent from Mail<https://go.microsoft.com/fwlink/?LinkId=550986> for Windows > 10 > > > [[alternative HTML version deleted]] > >__ >R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see

Re: [R] Change the position of label when using R package eulerr

2018-09-14 Thread David Winsemius
lot() to `my.plot` you get: grid.ls(my.plot) euler.diagram GRID.gTree.11 diagram.grob.1 fills.grob.1 fills.grob.2 fills.grob.3 fills.grob.4 fills.grob.5 fills.grob.6 fills.grob.7 edges.grob labels.grob quantities.grob

Re: [R] Change the position of label when using R package eulerr

2018-09-15 Thread David Winsemius
= int 4 - attr(*, "unit")= chr "native" #--- str( fit1.plot$children$GRID.gTree.12$children$diagram.grob.1$children$quantities.grob$y ) # 'unit' num [1:7] -2.569native 15.889native 4.732native -0.788native -7.157native ... - attr(*, "valid.unit

Re: [R] as.Date and ylim in empty plot

2018-09-19 Thread David Winsemius
PDF's are accepted. PNGs are not. -- David. > > Thank you so much for your input. > > Best regards > Ogbos > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-

Re: [R] as.Date and ylim in empty plot: RESOLVED

2018-09-22 Thread David Winsemius
> Best > Ogbos > > On Wed, Sep 19, 2018 at 11:34 PM David Winsemius > wrote: > > > On Sep 19, 2018, at 7:55 AM, Ogbos Okike wrote: > > > > Dear Experts, > > I generated the plot attached. Every other thing is OK except the black > > horizontal

Re: [R] cut{base}: is it a bug?

2018-09-24 Thread David Winsemius
; ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ >>> Jose Claudio Faria >>> UESC/DCET/Brasil >>> joseclaudio.faria at gmail.com >>> Telefones: >>> 55(73)3680.5545 - UESC >>> 55(73)99966.9100 - VIVO >>> 55(73)98817.6159 - OI >>&g

Re: [R] Printing standard notation and scientific notation in the same column of a dataframe

2018-09-27 Thread David Winsemius
> [[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-guid

Re: [R] Strange result for strptime with %p

2018-10-02 Thread David Winsemius
eather Spray > > ______ > 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 commented

Re: [R] Strange result for strptime with %p

2018-10-02 Thread David Winsemius
> On Oct 2, 2018, at 10:11 AM, David Winsemius wrote: > > >> On Oct 2, 2018, at 6:30 AM, Marc Girondot via R-help >> wrote: >> >> Dear members... are these results normal ? >> >> For the first one, no problem. I expected this: >> >&

Re: [R] Polygon

2018-10-02 Thread David Winsemius
t;> # Add the shaded area. >>> polygon(cord.x,cord.y,col='skyblue') >>> >> > > -- > st...@ntu.edu.tw (S.T. Yen) > > > [[alternative HTML version deleted]] > > __ > R-help@r-project

Re: [R] Help to find ksl data set

2018-10-03 Thread David Winsemius
t.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 'Any technology distinguishable from magic is insufficiently

Re: [R] Defining Variables from a Matrix for 10-Fold Cross Validation

2018-10-09 Thread David Winsemius
ta = train) : > variable lengths differ (found for 'x') So the warning above is probably a great clue to the source of this error. Morale of the tale: Always read the warnings, even if your code proceeds. > > [[alternative HTML version deleted]] > > __ > R-help@

Re: [R] DNAbin

2018-10-16 Thread David Winsemius
___ > 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 commented, minimal, self-contained, reproducible code. D

Re: [R] "logical indexing, " [was] match() question or needle haystack problem for a data.frame

2018-10-28 Thread David Winsemius
> On Oct 25, 2018, at 11:28 PM, Knut Krueger wrote: > > Am 25.10.18 um 16:13 schrieb peter dalgaard: >> >> Yes: x[!(x$A %in% y$B),] > > Ok thats in my opinion a little workaround > why?: > > There is an > = and != > < and > > > > means the opposite is available between terms. > > why is

Re: [R] Corrupted package?

2018-10-31 Thread David Winsemius
> On Oct 31, 2018, at 11:30 AM, peter dalgaard wrote: > > Hm, a source install to r-devel gave me > > Peter-Dalgaards-MacBook-Air:BUILD pd$ ls -l library/pkgconfig/R/pkgconfig.rdb > -rw-r--r-- 1 pd staff 4515 Oct 31 19:15 library/pkgconfig/R/pkgconfig.rdb > Peter-Dalgaards-MacBook-Air:BUI

Re: [R] cox model

2018-11-03 Thread David Winsemius
It's also "well described" in the help materials for the obvious recommended package that ships with every copy of R. My copy sits at http://127.0.0.1:29434/library/survival/doc/timedep.pdf. Therneau's S package was first ported to R by Thomas Lumley and later Therneau took over maintenance.

Re: [R] Predict follow up time using parametric model in r

2018-11-03 Thread David Winsemius
On 11/3/18 12:59 PM, Israel Ortiz wrote: I am trying to predict follow-up time using several survival models, both parametric and semi-parametric. I achieve it for semi parametric models using predict.coxph function in R from survival package using type = "expected" as indicated in help. Howeve

Re: [R] Predict follow up time using parametric model in r

2018-11-03 Thread David Winsemius
On 11/3/18 12:59 PM, Israel Ortiz wrote: I am trying to predict follow-up time using several survival models, both parametric and semi-parametric. I achieve it for semi parametric models using predict.coxph function in R from survival package using type = "expected" as indicated in help. Howeve

Re: [R] cox model

2018-11-03 Thread David Winsemius
It's also "well described" in the help materials for the obvious recommended package that ships with every copy of R. My copy sits at http://127.0.0.1:29434/library/survival/doc/timedep.pdf. Therneau's S package was first ported to R by Thomas Lumley and later Therneau took over maintenance.

Re: [R] list contingency tables

2018-11-09 Thread David Winsemius
Seems like you are trying to recreate the calculations needed to perform an exact test. Why not look at the code for that or even easier, just use the function. -- David. On 11/8/18 8:05 PM, li li wrote: Hi all, I am trying to list all the 4 by 2 tables with some fixed margins. For exa

Re: [R] Line with linearly changing thickness – installation issues

2018-11-11 Thread David Winsemius
I would have imagined that drawing a polygon would be the way most people would have attempted. Regarding Murrell's package: I thought the package name was "vwline". My attempt to install was unsuccessful> > devtools::install_github("pmur002/vwline") Error in utils::download.file(url, path,

Re: [R] semiparametric manova

2018-11-12 Thread David Winsemius
On 11/12/18 12:37 PM, Yectli Huerta via R-help wrote: > Hello, > > I was wondering if there are other packages like MANOVA.RM that could be used > to analysis non normal distributions. I have to analyze data with more than 2 > predictor variables and a similar number of response variables. When

Re: [R] wBoot Package

2018-11-25 Thread David Winsemius
Look at the function's help page: No help there. The "parameter" argument is not defined in any substantive manner, and no examples other than `parameter=mean` appear in the help page. (Now) Look at the code. The parameter argument is expected to be a function. There is no function named `pr

Re: [R] Confirm Adding to mailing list

2018-12-08 Thread David Winsemius
I cannot tell what is being asked but it seems clear that you are confused about how mailing lists work.. I've looked in the pending posts in the moderation queue and none of them appear to be from you. I saw an earlier misdirected "confirm" message that should have gone to the mailserver-bot a

Re: [R] General information

2018-12-08 Thread David Winsemius
Most of those answers are in the R-FAQ. It ships with every distribution copy from CRTAN and is also available online. Estimating the number of R users would be difficult. There are probably 20,000+ subscriblers to Rhelp. Not certain what is meant by #5. David. On 12/8/18 2:46 AM, ghada

Re: [R] Recursive Feature Elimination with SVM

2019-01-01 Thread David Winsemius
On 1/1/19 4:40 AM, Priyanka Purkayastha wrote: I have a dataset (data) with 700 rows and 7000 columns. I am trying to do recursive feature selection with the SVM model. A quick google search helped me get a code for a recursive search with SVM. However, I am unable to understand the first part

Re: [R] Recursive Feature Elimination with SVM

2019-01-01 Thread David Winsemius
producible example. See the Posting Guide ... and learn to post in plain text. -- David On Tue, 1 Jan 2019 at 11:42 PM, David Winsemius mailto:dwinsem...@comcast.net>> wrote: On 1/1/19 4:40 AM, Priyanka Purkayastha wrote: > I have a dataset (data) with 700 rows and

Re: [R] Merge the data from multiple text files

2019-01-05 Thread David Winsemius
On 1/5/19 7:28 AM, Priya Arasu via R-help wrote: I have multiple text files, where each file has Boolean rules. Example of my text file 1 and 2 Text file 1: A = not(B or C) B = A and C C = D Text file 2: A = D and E B = not(D) I want to merge the contents in text file as follows A = not(B or C

Re: [R] mgcv : 3-way interaction and 3D-plots ?

2019-01-07 Thread David Winsemius
On 1/7/19 3:35 PM, varin sacha via R-help wrote: Dear R-experts, I have fitted a model with 2-way and 3-way interactions. I would like, for the 3-way interaction (year,age,by=education), to obtain 3D-plots. How could I do that ? Forget ggplot2. It has ignored this sort of visualization effo

Re: [R] Running list of drugs taken and dropped (via Reduce and accumulate = TRUE or by other means)

2019-01-10 Thread David Winsemius
On 1/10/19 7:00 AM, Paul Miller via R-help wrote: Hello All, Would like to keep a running total of what drugs cancer patients have taken and what drugs have been dropped. Searched the Internet and found a way to cumulatively paste a series of drug names. Am having trouble figuring out how t

Re: [R] SURVDIFF()

2019-01-12 Thread David Winsemius
On 1/12/19 1:18 AM, Medic wrote: How to note (in code) a few (!) adjusting covariates for cox regression. I had an example for one covariate, and tried (according to my own understanding) two variantes of code (pls, see below), and got ... a different p-value. What is the right code? Many than

Re: [R] [FORGED] Newbie Question on R versus Matlab/Octave versus C

2019-01-28 Thread David Winsemius
On 1/28/19 4:00 PM, Alan Feuerbacher wrote: On 1/28/2019 4:20 PM, Rolf Turner wrote: On 1/29/19 10:05 AM, Alan Feuerbacher wrote: Hi, I recently learned of the existence of R through a physicist friend who uses it in his research. I've used Octave for a decade, and C for 35 years, but wo

Re: [R] Extract the coordinates of a Polylines

2019-01-31 Thread David Winsemius
On 1/30/19 7:12 AM, javad bayat wrote: > Dear all; > Back to my previous question, I am trying to add X and Y coordinates > to every row of the data. > > topo = > readOGR("E:/New/Modelling_Water/MIKE/BathyMetry/GIS_Armator/Chitgar_Topo.shp")#Read > > shape file of the topo as polylines > plot

Re: [R] Help converting file to XTS

2019-02-05 Thread David Winsemius
On 2/5/19 6:06 AM, Johannes Møllerhagen wrote: Hello there! I am a master student working on my master thesis, and I am trying to convert some data to xts so I can apply a highfrequency package to it. At the moment I am trying to use a POSIXct function. I am quite new at this program and I a

<    3   4   5   6   7   8   9   10   11   12   >