Re: [R] wanting to count instances of values in each cell of a series of simulated symmetric matrices of the same size

2021-06-01 Thread R. Mark Sharp via R-help
Bert, You are obviously correct about the diagonals. I was not thinking carefully. Typically they are expected to be at or near 0.5 in an outbred population but can theoretically go to 1.0 in completely inbred populations. This was subset from a madeup pedigree and I reused parents, hence the

Re: [R] wanting to count instances of values in each cell of a series of simulated symmetric matrices of the same size

2021-06-01 Thread Bert Gunter
Come again?! The diagonal values in your example are not all .5. If space is not an issue, a straightforward approach is to collect all the matrices into a 3d array and use indexing. Here is a simple reprex (as you did not provide one in a convenient form, e.g via dput()) x <- matrix(1:9, nr =

[R] wanting to count instances of values in each cell of a series of simulated symmetric matrices of the same size

2021-06-01 Thread R. Mark Sharp via R-help
I want to capture the entire distribution of values for each cell in a sequence of symmetric matrices of the same size. The diagonal values are all 0.5 so I need only the values above or below the diagonal. A small example with three of the structures I am wanting to count follows: F

Re: [R] replace NA into - for specific column

2021-06-01 Thread Jeff Newmiller
A unary negative sign with no number is not a number, so it has to be character. If you are done with computations you can format your numbers as character data and set the NA to "-", but such a conversion will prevent you from performing computations so it is only useful for creating report

[R] replace NA into - for specific column

2021-06-01 Thread Kai Yang via R-help
Hi List, I have a column MMR in a data frame proband_crc2. The column contents missing value and + (means positive).  I need to replace all missing value into - (means negative) I did try with difference ways.  Below are my testing, but not any one works. Can someone help me? Thanks, Kai

[R] [R-pkgs] r2rtf: Easily Create Production-Ready Rich Text Format (RTF) Table and Figure

2021-06-01 Thread Yilong Zhang
Dear all, I am happy to announce that {r2rtf} v0.3.0 is now on CRAN. r2rtf is an R package to create production-ready tables and figures in RTF format. The R package is designed to provide simple “verb” functions that correspond to each component of a table or figure. I hope you find it

Re: [R] DBDA2E-utilities.R file download

2021-06-01 Thread John Kane
I have added this to the R Help mailing list as it increases your chances of getting a helpful response. That error message can't "lazyload" indicates that you have some serious problem(s) upstream. Probably the best thing to do is copy and paste the entire installation output here for people

Re: [R] Binary data on JAGS

2021-06-01 Thread Massimiliano Tripoli
Thanks, for suggestions (To define p before). I thought it and I put this line of code inside the loop: p[i] <- y[i] / n[i] but I received a RUNTIME ERROR and: (Attempt to redefine node p[1]). Thanks Massimiliano Il giorno mar 1 giu 2021 alle ore 17:22 Bert Gunter ha scritto: > Where is p

Re: [R] Binary data on JAGS

2021-06-01 Thread Bert Gunter
Unless you have got reason not to, always reply to the list (included in this response). I cannot help, but someone else may be able to. Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom

[ESS] Broken parsing from compilation output in inferior-ess-r-mode

2021-06-01 Thread Stephen Berman via ESS-help
I have (define-key inferior-ess-r-mode-map "_" #'ess-insert-assign) in my init file. After installing R on a new system and starting ESS (built from the Github repo at commit e21ce9b260ff2ae0f162624ae3888c60bca3ae55) I called install.packages() and after that, typing `_' only inserted that

Re: [R] Binary data on JAGS

2021-06-01 Thread Bert Gunter
Where is p defined before it is used? (Is this part of what jags provides somehow?) Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Jun 1, 2021 at 7:57 AM

[R] Binary data on JAGS

2021-06-01 Thread Massimiliano Tripoli
Dear R users, I'm trying to reproduce the example 6.5.1 (Dobson (1983)) in BUGS book in linux using JAGS. Below the code as https://www.mrc-bsu.cam.ac.uk/software/bugs/the-bugs-project-the-bugs-book/bugs-book-examples/the-bugs-book-examples-chapter-6-6-5-1/ # By R code: library('rjags') jags

Re: [R] Plot GEE confindence interval band ggplot2

2021-06-01 Thread Rui Barradas
Hello, I don't know if the following is what you want but it gives confidence bars. The error in your code is to have Esp as id. dat$Id <- as.integer(factor(dat$Esp)) m5 <- geeglm( formula = tim ~ Pa*Pt, family = Gamma(link = log), data = dat, id = Id, corstr = "exchangeable" )

Re: [R] if statement and for loop question

2021-06-01 Thread Michael Dewey
Dear Kai When you ask again it is best to tell us what your input is and what output you were hoping for and what you actually got. If you can make a small data-set which shows all that then your post will be much more likely to get a helpful response. If you want to transfer the data-set to