Re: [R] Real Basic Question

2019-09-27 Thread David Winsemius
Sent from my iPhone > On Sep 27, 2019, at 12:55 AM, Phillip Heinrich wrote: > > Just when I think I’m starting to get the hang of R I run into something that > sends me back to Go without collecting $200. > > The working directory seems to be correct when I load an .rda file but it is >

Re: [R] Error: package 'stats' could not be loaded

2019-09-27 Thread David Winsemius
I’m having difficulty understanding why you installed R v 3.1.1. That’s five major versions ago. — David Sent from my iPhone > On Sep 26, 2019, at 2:51 AM, Carlos Weiss wrote: > > I'm getting an error when using R for the first time. It's a web service > written in C#, .NET 4.5, IIS 10.0

Re: [R] Plot two values on same axes

2019-09-26 Thread David Winsemius
Instead of trying to mix lattice and base functions, you might try using the formula: maxtemp+mintemp ~ sampdate And then: col= c(“red”, “blue”) Sent from my iPhone, so make sure those quotes are ordinary double quotes. — David > On Sep 27, 2019, at 6:27 AM, Rich Shepard wrote: > > I

Re: [R] & statement within an ifelse Loop

2019-09-21 Thread David Winsemius
On 9/21/19 12:57 PM, Phillip Heinrich wrote: Still putzing around trying to increment a count vector when the date changes. Date count 1 2018-03-29 1 2 2018-03-29 1 3 2018-03-29 1 81 2018-03-30 1 82 2018-03-30 1 83 2018-03-30 1 165 2018-03-31 1

Re: [R] Problem loading packages in R 3.6.1

2019-09-19 Thread David Winsemius
On 9/19/19 1:11 AM, Jeff Newmiller wrote: Since no answer has been forthcoming, my advice is to read the Posting Guide. There is a more appropriate list for your question, and you should take the warning about posting using plain text format on all r-project mailing lists seriously. On

Re: [R] bi-directional bar chart with a central axis

2019-09-17 Thread David Winsemius
On 9/17/19 3:14 PM, Sabatier, Jennifer F. (CDC/DDPHSIS/CGH/DGHP) via R-help wrote: Hi R-help, I have this data: my.dta <-data.frame(matrix(c( 26.3, 21.4, 20.1, 13.4, 7.9,3.9, 16.5, 14.6, 5.3,3.6, 38.6, 25.6, 34.4, 21.6, 77.4, 79.5, 58.2, 56.1, 80.5, 84, 37.7, 31.9,

Re: [R] Not the same length

2019-09-17 Thread David Winsemius
On 9/17/19 2:08 PM, David Winsemius wrote: On 9/17/19 1:35 PM, varin sacha wrote: Many thanks David, it perfectly works. Now, one last think. If I want my R code here below to run let's say B=500 times and at the end I want to get the average for the MSE_GAM and for the MSE_MARS. How can I

Re: [R] Not the same length

2019-09-17 Thread David Winsemius
z)+s(a)) mars_model<-earth(y_obs~x+z+a) MSE_GAM<-mean((gam_model$fitted.values - y_model)^2) MSE_MARS<-mean((mars_model$fitted.values - y_model)^2) MSE_GAM MSE_MARS Le mardi 17 septembre 2019 à 22:27:54 UTC+2, David Winsemius a écrit : On 9/17/19 12:48 PM, varin sacha via

Re: [R] Not the same length

2019-09-17 Thread David Winsemius
On 9/17/19 12:48 PM, varin sacha via R-help wrote: Dear R-helpers, Doing dput(x) and dput(y_obs), the 2 vectors are not the same length (1800 for y_obs and 2000 for x) How can I solve the problem ? Here is the reproducible R code   #  #  #  #  #  #  #  #  #  # library(mgcv)

Re: [R] about paired or unpaired t.test or wilcox.test in R

2019-09-11 Thread David Winsemius
On 9/11/19 1:14 PM, Bogdan Tanasa wrote: Dear all, if would be great if you could please advise on the use of PAIRED or UNPAIRED T.TEST or WILCOX.TEST in R : let's say shall we have 2 samples : -- CONTROL : where we measure the expression of 100 genes G1 ... G100 in one million cells C1

Re: [R] [R-devel] Source Code for function

2019-09-10 Thread David Winsemius
On 9/10/19 11:13 AM, Golden, Shelby wrote: Wow, this is great! Thank you, Richard, for going so in-depth with this explanation. My ultimate goal is to create a vector that stores logical evaluations to be used later in a function that subsets data. I’ve coded multiple projects, now, that

Re: [R] Moving 2nd column into 1st column using R

2019-09-09 Thread David Winsemius
On 9/9/19 7:19 PM, smart hendsome via R-help wrote: Hi R-user, I have a problem regarding R.  How can I move my 2nd column into 1st column.  For example I have data as below: mydf <- data.frame(matrix(1:6, ncol = 2))  mydf   X1 X2   1   4   2   5   3   6 I want move my 2nd column become

Re: [R] [R-devel] Source Code for function

2019-09-06 Thread David Winsemius
On 9/6/19 1:07 PM, Golden, Shelby wrote: Thank you all for your reply. I should clarify, that I am looking to understand why the keyword function can take a logical argument (eg: x<4) and use that later inside the function's definition for logical evaluations. Consider this example, which

Re: [R] Change Y-axis labels

2019-09-04 Thread David Winsemius
On 9/4/19 7:35 AM, Medic wrote: The Y scale is divided (by default) as: 0.0 ... 0.2 ... 0.4 ... 0.6 ... 0.8 ...1.0 But I would like so: 0 ... 20 ... 40 ... 60 ... 80... 100 (with rotating axis labels) When I use par function (marked as comment here) it turns out correctly for ONLY ONE

Re: [R] New Work Based on R

2019-08-28 Thread David Winsemius
On 8/28/19 6:13 AM, Ogbos Okike wrote: Dear Contributors, Some of us that use the end product of R could pay for your services if asked to do so. While all your help is free, I am a little disturbed that there is no archive or repository where publications employing R tools are deposited for

Re: [R] A goodness of fit test for two discrete distributions with unequal variance?

2019-08-23 Thread David Winsemius
On 8/23/19 2:52 PM, Serena De Stefani wrote: I have a computer simulation in which a virtual agent end up in different areas of a layout based on several factors. There are 18 conditions in total. If I collapse the datapoint into bins, where each bin is one of the areas, the data would look

Re: [R] Unable to load GRanges

2019-08-17 Thread David Winsemius
Looks like you do not have permissions to write to that directory, but it's not clear that this is a major problem. The updating of boot and foreign is not something that regular non BioC users would have usually needed.  Did you recently update your R installation? Library locations sometimes

Re: [R] Trying to understand how to sort a DF on two columns

2019-08-12 Thread David Winsemius
On 8/12/19 7:20 PM, Sorkin, John wrote: I want to sort a DF, temp, on two columns, patid and time. I have searched the internet and found code that I was able to modify to get my data sorted. Unfortunately I don't understand how the code works. I would appreciate it if someone could explain

Re: [R] Help

2019-08-11 Thread David Winsemius
BSCRIBE 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, US

Re: [R] ggplot2 fill problem in colour shading

2019-08-11 Thread David Winsemius
ernative 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 > an

Re: [R] Using read.table for importing gz file

2019-08-10 Thread David Winsemius
read system.time( z <- fread('https://TCGA.xenahubs.net/download/TCGA.GBMLGG.sampleMap/HumanMethylation450.gz', sep="\t")  ) #--- [100%] Downloaded 597770433 bytes...    user  system elapsed  20.682   3.322  29.292 > dim(z) [1] 485577    686 -- David. On 8/10/19 5:

Re: [R] Using read.table for importing gz file

2019-08-10 Thread David Winsemius
o R so that I can proceed with a survival analysis. Also, what particular rule of the mailing list did I break? I apologize in advance, as I thought that code specific queries like the one I asked were acceptable. Many thanks, Spencer On Sat, Aug 10, 2019 at 5:51 PM David Winsemius mailto

Re: [R] Using read.table for importing gz file

2019-08-10 Thread David Winsemius
Have you tried using readLines in the manner illustrated on the ?gzfile help page? David. On 8/10/19 12:29 PM, Spencer Brackett wrote: Hello, I am trying to read the following Xena dataset into R for data analysis:

Re: [R] broken installation on ubuntu 16.04 for 3.4.4-1xenial0

2019-08-08 Thread David Winsemius
The proper place for Ubuntu/Debian questions is r-sig-deb...@r-project.org -- David On 8/8/19 12:52 PM, k...@breadfinance.com wrote: Sorry for unintentionally sending a rich text email. I have confirmed that installation follows ubuntu guide lines in your link and I reproduce the output of

Re: [R] Extract row as NA with no matching name

2019-08-08 Thread David Winsemius
I don't know a clean way of delivering that result but if you use logical indexing you can get an empty matrix with three columns: str( mdat["nope" %in% rownames(mdat), ] )  num[0 , 1:3]  - attr(*, "dimnames")=List of 2   ..$ : NULL   ..$ : chr [1:3] "C.1" "C.2" "C.3" # it prints thus to the

Re: [R] CoxPH multivariate frailty model with coxph (survival)

2019-07-28 Thread David Winsemius
On 7/19/19 10:19 AM, Denise b wrote: Dear R users, I am interested in estimating the effects of a treatment on two time-to-event traits (on simulated data), accounting for the dependency between the two time-to-event outcomes. I precise that the events are NOT recurrent, NOT competitive,

Re: [R] How to rename file names with condition

2019-07-23 Thread David Winsemius
On 7/23/19 6:30 AM, danielmessay--- via R-help wrote: Could somebody please help me on this? I have many files on my computer and would like to change its names so that it would be more meaningful for me. The filenames are in the following format: cvcvcv198307.xlsx I want to change it in

Re: [R] List of tables needed

2019-07-20 Thread David Winsemius
On 7/20/19 5:17 PM, Jackson Rodrigues wrote: Dear all, My name is Jackson and I need a help in applying functions to a list of tables. Could anyone help me to use loop/array on a list of tables ? I need to apply a group of functions (listed below) to a list of tables I am not that good

Re: [R] Call `rlang::last_error()` to see a backtrace

2019-07-10 Thread David Winsemius
Did you read the bug report from three weeks ago and the suggested fix documented on the webpage you cited? — David Sent from my iPhone > On Jul 10, 2019, at 1:09 PM, Patrick (Malone Quantitative) > wrote: > > First response: The ID column in your data is labeled "eid" but your > function

Re: [R] Output for pasting multiple vectors

2019-06-25 Thread David Winsemius
On 6/25/19 6:29 AM, Spencer Brackett wrote: Mr. Barradas, My apologies for the delayed response. No, (meth) is not a dataset within CRAN. I’m not sure why my supervisor wrote that in as the object for the bit of script I shared previously. Assuming that the correct object for this

Re: [R] How to extract data based on first column ID in R

2019-06-23 Thread David Winsemius
t how I can get the row values along with the header? > > > Thanks > > Yogesh > > [[alternative HTML version deleted]] > > ______ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/

Re: [R] Problem witth nnet:multinom

2019-06-23 Thread David Winsemius
he > same issue (just a different number of iterations). > > Thanks in advance. > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posti

Re: [R] Output for pasting multiple vectors

2019-06-21 Thread David Winsemius
ut together a small example and the show what your goals are. — David. >> On Fri, Jun 21, 2019 at 9:47 PM David Winsemius >> wrote: >> All the +-signs are telling you that the expression is not complete. Please >> read the posting guide. I’m pretty sure you’ve alr

Re: [R] Output for pasting multiple vectors

2019-06-21 Thread David Winsemius
All the +-signs are telling you that the expression is not complete. Please read the posting guide. I’m pretty sure you’ve already been warned NOT to use html. — David. Sent from my iPhone > On Jun 21, 2019, at 6:11 PM, Spencer Brackett > wrote: > > Hello, > > I am attempting to paste

Re: [R] Recovering former projects on R

2019-06-21 Thread David Winsemius
On 6/21/19 2:34 PM, Spencer Brackett wrote: Thank you, I called them but was able to load what appears to be a reprex of my Rhistory. If I remember correctly, there is a way to to reimplement this history into your console so to regenerate previous work. How is this done? Assuming you still

Re: [R] Plotting more than one regression line in ggplot

2019-06-05 Thread David Winsemius
On 6/5/19 9:57 AM, David Winsemius wrote: On 6/5/19 9:37 AM, rain1...@aim.com wrote: Hi David (and everyone), Thank you for your response. I changed the column names to x and y, but the error/warning persists: Warning message: Computation failed in `stat_smooth()`: 'what' must be a function

Re: [R] Plotting more than one regression line in ggplot

2019-06-05 Thread David Winsemius
.1 I'm also running: RStudio Version 1.1.463 – © 2009-2018 RStudio, Inc. You should now restart a clean session, try again with just the required packages and report back with full code and data. Best; David > > > -----Original Message- > From: David Winsemius > To: r-help

Re: [R] Plotting more than one regression line in ggplot

2019-06-05 Thread David Winsemius
On 6/5/19 8:04 AM, rain1290--- via R-help wrote: Hi Jeff (and everyone), Thank you for your response and feedback. Yes, I know what you mean - it was a blind and quick choice to use "lm" as my object name. Unfortunately, changing the object name to something else does not eliminate that

Re: [R] R demographic facts

2019-06-03 Thread David Winsemius
On 6/3/19 7:15 AM, Paul Bernal wrote: Dear friends, Hope you are all doing great. I would like to know if there is any place where I can find R demographic facts, such as, number of R users worldwide, by region, by country, by sector, etc. Over the years there have been several efforts to

Re: [R] how to separate string from numbers in a large txt file

2019-05-16 Thread David Winsemius
x, it happened in the switch from 84 to 85, I guess. It went on like [990] "//1,//2,//3,//4" [991] "//1,//2,//3,//4" [992] "//1,//2,//3,//4" [993] "//1,//2,//3,//4" [994] "//1,//2,//3,//4" [995] "//1,//2,//3,//4" [996] "

Re: [R] how to separate string from numbers in a large txt file

2019-05-16 Thread David Winsemius
aybe, the "seriously" will be in a column all to itself, as will be the "I've'"and the "never" etc. I will use a regular expression if I have to, but it would be nice to keep the dates and times on there. Originally, I thought they were meaningless, but I've since changed my

Re: [R] how to separate string from numbers in a large txt file

2019-05-15 Thread David Winsemius
On 5/15/19 4:07 PM, Michael Boulineau wrote: I have a wild and crazy text file, the head of which looks like this: 2016-07-01 02:50:35 hey 2016-07-01 02:51:26 waiting for plane to Edinburgh 2016-07-01 02:51:45 thinking about my boo 2016-07-01 02:52:07 nothing crappy has happened, not

Re: [R] Trying to understand the magic of lm (Still trying)

2019-05-10 Thread David Winsemius
On 5/10/19 12:53 PM, Sorkin, John wrote: A number of people have helped me in my mission to understand how lm (and other fucntions) are able to pass a dataframe and then refer to a specific column in the dataframe. I thank everyone who has responded. I now know a bit about

Re: [R] Question concerning Loading large Data Sets

2019-05-05 Thread David Winsemius
I wouldn’t have expected an message, but it does raise the question: why are you making two different copies of the same text file if you are concerned about size issues? — David Sent from my iPhone > On May 5, 2019, at 5:14 PM, Spencer Brackett > wrote: > > Also, > > In case there is

Re: [R] Passing formula as parameter to `lm` within `sapply` causes error [BUG?]

2019-04-30 Thread David Winsemius
Try using do.call — David Sent from my iPhone > On Apr 30, 2019, at 9:24 AM, Jens Heumann > wrote: > > Hi, > > `lm` won't take formula as a parameter when it is within a `sapply`; see > example below. Please, could anyone either point me to a syntax error or > confirm that this might be

Re: [R] Problems w/ creating object

2019-04-22 Thread David Winsemius
On 4/22/19 11:49 AM, Spencer Brackett wrote: Hello R users, I am trying to create an object out of some data a colleague sent my way, so to duplicate the following code... library(data.table) anno = as.data.frame(fread(file = "/rsrch1/bcb/kchen_group/v_mohanty/data/TCGA/450K/mapper.txt", sep

Re: [R] Help installing netReg

2019-04-14 Thread David Winsemius
On 4/14/19 3:20 PM, Spencer Brackett wrote: Good evening, I am having problems with downloading the package used to generate regression models on R. The following is the error message I received. I tried installing BiocManager instead as suggested, but this too did not work. Any ideas? The

Re: [R] [FORGED] Re: Quantile Density Contours

2019-03-28 Thread David Winsemius
has a "sliced" approach to disply of 2d densities that I found informative and attractive. -- David. Bernard McGarvey Director, Fort Myers Beach Lions Foundation, Inc. Retired (Lilly Engineering Fellow). On March 28, 2019 at 1:40 PM David Winsemius wrote: On 3/27/

Re: [R] [FORGED] Re: Quantile Density Contours

2019-03-28 Thread David Winsemius
On 3/27/19 3:43 PM, Bernard Comcast wrote: To follow on Jeff, is there a function to do 2-D (double) numerical integration in R? Packages pracma and cubature offer a variety of solutions to that task. -- David. Bernard Sent from my iPhone so please excuse the spelling!" On Mar 27,

Re: [R] Regression Modeling Strategies and the R rms Package Short Course 2019

2019-03-24 Thread David Winsemius
Hi Graeme; I took the course about ten years ago. I did so after getting a Masters in Epidemiology from the University of Washington and doing very well in all my stats courses and submitting my thesis work on solving regression problems with stratified sampling using bootstrap methods. So

Re: [R] counting unique values (summary stats)

2019-03-21 Thread David Winsemius
On 3/21/19 3:31 PM, reichm...@sbcglobal.net wrote: r-help I have the following little scrip to create a df of summary stats. I'm having problems obtaining the # of unique values unique=sapply(myData, function (x) length(unique(x), replace = TRUE)) I just looked

Re: [R] Rank ANCOVA

2019-03-21 Thread David Winsemius
If you are looking for robust multivariate or multivariable methods, then review the Robust Methods Task View: https://cran.r-project.org/view=Robust -- David. On 3/13/19 4:06 AM, Jackson, Daniel wrote: Hi Frank and Dennis I am in a similar situation but I would prefer to use a

Re: [R] A general question about using Bayes' Theorem for calculating the probability of The End of Human Technological Civilisation

2019-03-19 Thread David Winsemius
On 3/19/19 12:49 PM, Jeff Newmiller wrote: Highly off topic. Try StackOverflow. As it stands it's off-topic for SO. (You would just be making more work for those of us who know the rules but need 4 close votes for migration.)  Better would be immediately posting at CrossValidated.com

Re: [R] A general question about using Bayes' Theorem for calculating the probability of The End of Human Technological Civilisation

2019-03-19 Thread David Winsemius
), and if these 'statistics discussion' questions are inappropriate to the mission (as described), then...why isn't the 'moderator' (him/her/they) blocking on submission? On 3/19/2019 1:59 PM, David Winsemius wrote: Rhelp is not a forum for discussions of statistics. Instead it is for persons

Re: [R] A general question about using Bayes' Theorem for calculating the probability of The End of Human Technological Civilisation

2019-03-19 Thread David Winsemius
Rhelp is not a forum for discussions of statistics. Instead it is for persons who have specific questions about the use of R. Please read the list info page where you started the subscription process. And do read the Posting Guide. Both these are linked at the bottom of this response. There

Re: [R] About wilcox test

2019-03-12 Thread David Winsemius
On 3/12/19 9:23 AM, Thomas Subia via R-help wrote: Javid wrote: "I have two set of data in excel: A column( 16.38, -31, -16.77, 127, -57, 23.44 and so on) B column ( -12, -59.23, -44, 34.23, 55.5, -12.12 and so on) If those are in an excel spreadsheet than you need to indicate exactly how

Re: [R] How to change the number of bins?

2019-03-10 Thread David Winsemius
one to specify a value such as 999 to be missing but R needs to have it changed to NA is.na(Table$pdays) <- Table$pdays == 999 -- David > > David Winsemius <mailto:dwinsem...@comcast.net>> 於 2019年3月10日 週日 下午10:48寫道: > > Seems rather likely that 999 is not r

Re: [R] How to change the number of bins?

2019-03-10 Thread David Winsemius
Seems rather likely that 999 is not really a measured value but rather is a missing value indicator. -- David. On 3/10/19 1:54 PM, wong bowie wrote: I wish to calculate the weight of evidence of a variable x, which is positively skewed, with over 6000 of the observations are 999 but only

Re: [R] add points to lattice cloud plot (3D scatter)

2019-02-28 Thread David Winsemius
On 2/28/19 7:56 AM, Jeff Newmiller wrote: You are missing the point... lattice assembles the entire data set at once so it can adjust and synchronize all of the scales and then it generates an object that can be printed to a device. This approach is entirely incompatible with the base

Re: [R] Help formulating gamm with repeated measures and spatial autocorrelation using mgcv package

2019-02-26 Thread David Winsemius
Crossposting is explicitly advised against in the rhelp posting guide. Furthermore, this is the wrong list for this sort of question within the r-project hierarchy of mailing lists. I do not believe an answer is possible in the absence of the data (and note the request some would say

Re: [R] Saving and reloading function in a package

2019-02-15 Thread David Winsemius
Isn't this much more on topic with the package development list? -- David > On Feb 15, 2019, at 6:37 AM, Dimitrios Stasinopoulos > wrote: > > I would like to put a graphic background to a model diagnostic plot. > The background is created with plot()/lines() but it takes time. > My

Re: [R] Using ggplot2 geom_path() in a grouped variable

2019-02-12 Thread David Winsemius
To Norberto; Your code _probably_ would have succeeded if you had used color as a grouping argument, or you could have used group or linetype or probably others, but I find locating the listing of ggplot2 "behavioral" parameters rather frustrating. These details are not to be found in any of

Re: [R] Using ggplot2 geom_path() in a grouped variable

2019-02-12 Thread David Winsemius
> On Feb 12, 2019, at 2:09 PM, Norberto Hernandez > wrote: > > Hi! I am trying to make a scatter/path graph in one variable that is > divided in two groups (variable Control), but when I use the > geom_path() option, the line continues from the group one to the group > two, and I wasn't able

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

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 >

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

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

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

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

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

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] 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] 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,

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

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

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.

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] 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

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.

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.

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] 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 >

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] 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] Defining Variables from a Matrix for 10-Fold Cross Validation

2018-10-09 Thread David Winsemius
ble 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@r-project.org mailing list -- To UN

Re: [R] Help to find ksl data set

2018-10-04 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 advan

Re: [R] Polygon

2018-10-02 Thread David Winsemius
shaded area. >>> polygon(cord.x,cord.y,col='skyblue') >>> >> > > -- > st...@ntu.edu.tw (S.T. Yen) > > > [[alternative HTML version deleted]] > > ______ > R-help@r-project.org mailing list -- To UNSUBS

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] Strange result for strptime with %p

2018-10-02 Thread David Winsemius
r 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, min

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-g

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 >>> ///\\\///\\\///\\\///\\\///\\

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

2018-09-22 Thread David Winsemius
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 lines which sho

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

2018-09-19 Thread David Winsemius
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-help >

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

2018-09-15 Thread David Winsemius
it1.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")= int 4 - attr(*, "unit")= chr "native" -- David > t

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

2018-09-14 Thread David Winsemius
ot;), ) But I don't think that will succeed since you never assigned the value of the plot operation to a name. Instead you pulled out part of the grid object that was sitting "free" and unassigned to a name. If you assign that value of plot() to `my.plot` you get: grid.ls(my

Re: [R] New to R

2018-09-14 Thread David Winsemius
t; for Windows > 10 > > > [[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 postin

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] Help with setting locale

2018-09-14 Thread David Winsemius
ns 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 read t

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://www.R-project.o

<    1   2   3   4   5   6   7   8   9   10   >