Re: [R] Setting up hypothesis tests with the infer library?

yourself but if you get stuck come back with your code and what went wrong. There are many other options in packages which you could install but what I suggest should work out of the box. Michael On 29/03/2025 16:09, Kevin Zembower via R-help wrote: Hello, all, We're now starting to

Re: [R] Remove all box around a plot except bottom line (base graphics)

Thank you Duncan, I will try that next. Michael On 28/11/2024 13:52, Duncan Murdoch wrote: On 2024-11-28 8:36 a.m., Michael Dewey wrote: To make the plot clearer I have removed the axes but I wish to remove all the boz except the bottom horizontal line. Using the bty parameter does not seem

[R] Remove all box around a plot except bottom line (base graphics)

d side and the base. I could also remove the entire box but that looks odd. Am I reading the documentation incorrectly or is it really impossible? -- Michael __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailma

Re: [R] JASP vs R

end the design beyond the standard suite for descriptives, linear models, etc --- Michael Friendly Email: friendly AT yorku DOT ca Professor, Psychology Dept. York University Voice: 416 736-2100 x66249 4700 Keele StreetWeb: http://www.datavis.ca<http://www.datavis.ca/> | @da

Re: [R] OFF TOPIC: Nature article on File Drawer Problem in Reserach

discussion on the mailing list dedicated to meta-analysis and interested readers may want to search its archives https://stat.ethz.ch/pipermail/r-sig-meta-analysis/ Michael On 25/07/2024 23:40, Robert Baer wrote: Chapter 9 might be of interest: https://bookdown.org/MathiasHarrer

Re: [R] Integration of functions with a vector argument

Dear Bert, Thank you very much for your suggestions! Yours sincerely, Michael Michael Levine Associate Professor, Statistics Department of Statistics Purdue University 250 North University Street West Lafayette, IN 47907 USA email: mlev...@purdue.edu Phone: +1-765-496-7571 Fax: +1-765-494

Re: [R] Integration of functions with a vector argument

Dear Ivan, Thank you very much again for your suggestions! I will try to use Cubature, then. yours sincerely, Michael Michael Levine Associate Professor, Statistics Department of Statistics Purdue University 250 North University Street West Lafayette, IN 47907 USA email: mlev...@purdue.edu

Re: [R] Regression performance when using summary() twice

Dear Christian Without knowing how big your datset is it is hard to be sure but confint() can take some time. Have you thought of calling summary once summ <- summary(model) and then replace all subsequent calls to summary with summ Michael On 21/06/2024 15:38, c.bu...@posteo.jp wr

Re: [R] Integration of functions with a vector argument

f any detailed discussion of this package beyond the two vignettes available on CRAN? Thank you very much again for your help! Yours sincerely, Michael Michael Levine Associate Professor, Statistics Department of Statistics Purdue University 250 North University Street West Lafayette, IN 47907 USA

Re: [R] Integration of functions with a vector argument

go about performing this task? Yours sincerely, Michael Michael Levine Associate Professor, Statistics Department of Statistics Purdue University 250 North University Street West Lafayette, IN 47907 USA email: mlev...@purdue.edu Phone: +1-765-496-7571 Fax: +1-765-494-0558 URL: www.stat.

[R] Integration of functions with a vector argument

with a vector argument x. At the very least, I haven't seen any good examples of this being done. Any suggestions? Yours sincerely, Michael Michael Levine Associate Professor, Statistics Department of Statistics Purdue University 250 North University Street West Lafayette, IN 479

Re: [R] R-help Digest, Vol 255, Issue 17

You might be interested in the `Rdatasets` package, https://vincentarelbundock.github.io/Rdatasets/ which lists over 2200 datasets from various packages. What is the context of the `lottery` dataset. I seem to recall smth to do with the NJ Lottery -Michael 1. Availability of Sdatasets

Re: [R] Strange variable names in factor regression

^5 17.406593 -12.036715 -4.968654 -1.852819 3.291477 4.263642 wday^6 2.591317 Yes, that is what ordered is supposed to do, fit polynomial contrasts. When you remove the intercept that breaks it so you get a different fit, in fact the same as you got when it was not ordered. Michael

Re: [R] Double buffering plots on Windows

: [R] Double buffering plots on Windows Message-ID: mailto:b74c68da-a0b2-47dd-b54f-6b318488c...@stat.auckland.ac.nz>> Content-Type: text/plain; charset="utf-8"; Format="flowed" Hi Take a look at dev.hold() and dev.flush() Paul --- Michael Friendly Email: friendly AT

[R] DescTools::Quantile

quantiles are equal to 100! Is this function working or am I not using it correctly? Michael Meyer __ 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

[R] Quantiles of sums of independent discrete random variables

/10, FUN.VALUE=0 ))) g <- function(s) dnorm(s,mean=0,sd=0.01) cat("\n\n") for(i in 1:5){       cat("Doing convolution ",i,"\n")   g <- convolve(g,f)$density } cat("\nConvolutions finished, plotting density.") s <- seq(0,100,length.out=1024) matplot(s,g(s)

[R] Use of geometric mean .. in good data analysis

ic mean all asymptotic results no longer apply. Michael Meyer [[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 h

[R] Quantiles of sums of independent discrete random variables

ideas how I could approach this problem with the tools of R? Thanks in advance for all replies, Michael Meyer [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/ma

Re: [R] Function with large nested list

he final parenthesis or given it an incomplete line. Michael As I said, i have checked the structure multiple times, but cant find an error. I have tried rebuilding it and testing each time i add a part. Each part functions isolated, but not together in the same function. I can't find any

Re: [R] [External] Re: Error in setwd(dir) when initializing R

free of problems. Have you asked around? Michael On 21/11/2023 12:41, Ana de las Heras Molina wrote: Hello, I am sorry for my ignorance, but what is Rgui.exe? El mar, 21 nov 2023 a las 12:06, Ivan Krylov () escribió: В Tue, 21 Nov 2023 10:51:59 +0100 Ana de las Heras Molina пишет: I

Re: [R] car::deltaMethod() fails when a particular combination of categorical variables is not present

Thank you very much, John. This has allowed us to move forward on several fronts and better understand our data. - Michael Cohn On Tue, Sep 26, 2023 at 8:39 AM John Fox wrote: > Dear Michael, > > My previous response was inaccurate: First, linearHypothesis() *is* able > to accommo

Re: [R] Question about R software and output

Dear Charity Since your organisation is a member of King's Health Partners you might like to ask colleagues in KCL for local support. Michael On 02/10/2023 08:48, Ferguson Charity (CEMINFERGUSON) wrote: To whom it may concern, My understanding is that the R software is downloaded f

[R] car::deltaMethod() fails when a particular combination of categorical variables is not present

to work as expected. I don't know of any theoretical reason that missing one interaction parameter estimate should prevent getting a confidence interval for a different combination of predictors. Is there a way to use deltaMethod() or some other function to do this without changing my data?

Re: [R] How to fix this problem

It looks here as though the E coli column has commas in it so will be treated as character. Michael On 25/09/2023 15:45, avi.e.gr...@gmail.com wrote: David, This may just be the same as your earlier problem. When the type of a column is guessed by looking at the early entries, any non

Re: [R] Odd result

Dear David To get the first 46 rows just do KurtzData[1:43,] However really you want to find out why it happened. It looks as though the .csv file you read has lots of blank lines at the end. I would open it in an editor to check that. Michael On 23/09/2023 23:55, Parkhurst, David wrote

Re: [R] prop.trend.test

Dear Thomas Are you looking for more than smokers / patients Michael On 07/09/2023 14:23, Thomas Subia via R-help wrote: Colleagues  Consider smokers  <- c( 83, 90, 129, 70 ) patients <- c( 86, 93, 136, 82 )  prop.trend.test(smokers, patients)  Output: Chi-squared Te

Re: [R] To UNSUBSCRIBE

You can find instructions on how to do this on the bottom of e-mailas from the list. On 25/06/2023 15:34, Siddharth Arun wrote: Please unsubscribe me from your mailing list. -- Michael http://www.dewey.myzen.co.uk/home.html __ R-help@r

[R] post tweets with media attachment - and with retirement of rtweet

All - I run an automated script that pulls in some market data, makes some charts, and posts them to twitter. A couple of them per day, and then once per month there's a little flurry around the CPI report. For a while rtweet() worked great. Then v2 happened and Oauth2.0. The package broke, wa

Re: [R] Error on using Confint function to calculate 95% CI

may be able to help you better. Michael On 26/04/2023 02:37, bharat rawlley via R-help wrote: Hello, I am trying to estimate 95% CI using the confint function for a generalized liner model. While I am able to estimate Odds ratio using the coef function but on using the confint function, I get

Re: [R] nlme varFixed

Dear Franz Only attachments of specific types are allowed on the R lists and yours did not come through. It might be necessary to show your code as well. Have you thought of trying the mailing list dedicated to mixed models? https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models Michael

Re: [R] dyn.load(now = FALSE) not actually lazy?

ymbol: work_mem Execution halted On Wed, Feb 1, 2023 at 9:05 PM Ivan Krylov wrote: > В Wed, 1 Feb 2023 14:16:54 +1100 > Michael Milton пишет: > > > Is this a bug in the `dyn.load` implementation for R? If not, why is > > it behaving like this? What should I do about it? > &

[R] dyn.load(now = FALSE) not actually lazy?

On Linux, if I have a .so file that has a dependency on another .so, and I `dyn.load(now=FALSE)` the first one, R seems to try to resolve the symbols immediately, causing the load to fail. For example, I have `libtorch` installed on my HPC. Note that it links to various libs such as `libcudart.so`

Re: [R] My forest plot is not fit to windows in R software

Dear Fahimeh Try preceding your call of forest.meta with pdf("myplot.pdf") forest.meta(..) dev.off() Why did you set xlim to c(0, 1) Michael On 10/01/2023 19:37, Fahimeh Alizadeh wrote: Dear Prof. Thank you for your kind response. I have only used: install.packages(

Re: [R] My forest plot is not fit to windows in R software

different graphics format like pdf? Then we may be able to see where the issue lies. Michael On 10/01/2023 12:30, Fahimeh Alizadeh via R-help wrote: Dear Prof. My forest plot is not fit to windows in R software, I have 49 studies but forest plot only show from 9 to 42. forest.meta(hidemeta

Re: [R] Functional Programming Problem Using purr and R's data.table shift function

Dénes, thank you for the guidance - which is well-taken. Your side note raises an interesting question: I find the piping %>% operator readable. Is there any downside to it? Or is the side note meant to tell me to drop the last: "%>% `[`"? Thank you, == Michael Lachansk

[R] Functional Programming Problem Using purr and R's data.table shift function

. Thank you for any help you can provide or general guidance. == Michael Lachanski PhD Student in Demography and Sociology MA Candidate in Statistics University of Pennsylvania mikel...@sas.upenn.edu [[alternative HTML version deleted]] __ R-he

[R] piecewiseSEM two errors: arguments imply different numbers of rows / object 'ret' not found

Dear Rlist members, We are trying to run a SEM with your package. Unfortunaley while running our code two error messages popped up. For the life of us we could not figure out how to solve these errors . We are hoping very much that somebody can help us: Why do we get these error messages? We built

Re: [R] Select dataframe row containing a digit

I suspect [[:digit:]] might have done what you want. Michael On 30/11/2022 13:04, Luigi Marongiu wrote: Thank you, I have been trying with [:digit:] but did not work. It worked with `df$val[grepl('[0-9]', df$val)] = "NUM"` On Wed, Nov 30, 2022 at 2:02 PM Ivan Krylov wr

[R] picewiseSEM undefined columns selected

Dear r-help list members, I would like to run a pathway analysis using the R-package piecewiseSEM (v.2.12). I used a very simple model structure that is similar to the example in the piecewiseSEM description. Unfortunatley, once I execute the "as.psem" function I get the error "Error in `[.data

Re: [R] Sensitivity and Specificity

So predict is a one-dimensional vector of predictions but you are treating it as a two-dimensional matrix (presumably you think those are the totals). Michael On 24/10/2022 16:50, greg holly wrote: Hi Michael, I appreciate your writing. Here are what I have after; > predict_test

Re: [R] Sensitivity and Specificity

Rather hard to know without seeing what output you expected and what error message you got if any but did you mean to summarise your variable predict before doing anything with it? Michael On 24/10/2022 16:17, greg holly wrote: Hi all R-Help , After partitioning my data to testing and

Re: [R] Forestplot, grid graphics Viewplot grid.arange

There are several other CRAN packages which provide forest plots (see CRAN Task View for details) and they do not all use grip graphics which I think the forestplot package does. It might be worth swapping to one of them. Michael On 14/10/2022 04:34, Jim Lemon wrote: Hi Mary, I didn'

Re: [R] High concurvity/ collinearity between time and temperature in GAM predicting deaths but low ACF. Does this matter?

ly, but it seems to me that, as usual, we have a mix of statistical and scientific questions here and you may need to rethink your sceientific model in the light of the issue with your data. Michael On 10/06/2022 10:30, jade.sho...@googlemail.com wrote: Hi Michael, I don't think my reply to

Re: [R] High concurvity/ collinearity between time and temperature in GAM predicting deaths but low ACF. Does this matter?

try. Michael On 08/06/2022 18:15, jade.shodan--- via R-help wrote: Hi Simon, Thanks so much for this!! I have two follow up questions, if you don't mind. 1. Does including an autoregressive term not adjust away part of the effect of the response in a distributed lag model (where the ou

Re: [R] Shiny app Deployment not Working

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. -- Michael http://www.dewey.myzen.co.uk/home.html _

Re: [R] Restoration of "rite" package of R as R-script editor or the like

Notepad++ is fine as a light weight editor with syntax highlighting but the plug-in for R does not seem to have been updated for years. Michael On 23/05/2022 17:37, Rui Barradas wrote: Hello, I was also thinking about Notepad++. And about an editor no one mentinoed, vim. I've learned

Re: [R] How to create density ellipses with R

ellipse draws mathematical ellipses, I believe For proper bivariate normal density ellipses (not simply contours), check out: car::dataEllipse heplots::covEllipses -Michael Message: 1 Date: Fri, 14 Jan 2022 10:12:50 -0500 From: Paul Bernal To: R Subject: [R] How to create density ellipses

Re: [R] how to run r biotools boxM terst on multiple groups?

iables (such as, death, fever, symptoms etc). > I would like to run a box matrix test and I am using biotools' boxM, > but it Michael Friendly Email: friendly AT yorku DOT ca Professor, Psychology Dept. & Former Chair, ASA Statistical Graphics Section York Universi

Re: [R] Issue about eRM package

I do not know that package so cannot help with that but the dataset did not come through. This mailing list is quite restrictive in what sorts of attachment it allows so I suggest trying something like a .csv file or a .txt file. Michael On 28/09/2021 10:37, 宋启发 wrote: I am using your

Re: [R] What are the pros and cons of the log.p parameter in (p|q)norm and similar?

Sent off-list Thanks Bill and Duncan. I only asked for advice but I got an education too. Michael On 03/08/2021 21:24, Bill Dunlap wrote: In maximum likelihood problems, even when the individual density values are fairly far from zero, their product may underflow to zero. Optimizers have

[R] What are the pros and cons of the log.p parameter in (p|q)norm and similar?

expanding the documentation. -- Michael http://www.dewey.myzen.co.uk/home.html __ 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

Re: [R] if statement and for loop question

us then using dput() will ensure that we have exactly the same as you have. Do not forget to tell us what libraries, if any, you have loaded too. Michael On 31/05/2021 17:26, Kai Yang via R-help wrote: Hi Jim, Sorry to post "same" question, because 1. I was asking to use plain t

Re: [R] calculating area of ellipse

Dear Stephen In that application the axes would be sensitivity and specificity (or their inverses) or some transformation of them like logits so the units would be the same. Whether the area has any scientific meaning I am not sure. Michael On 11/05/2021 15:20, Stephen Ellison wrote: In

[R] mgcv: bam warning messages and non-convergence

eason why it does not onverge? And maybe a zero inflated negative binomial would better but I’m not sure how to undertake that. I know there’s a lot here but any help would be appreciated. Many thanks, Mike Michael Williamson London NERC DTP Candidate Email: michael.william...@kcl.ac.uk<

Re: [R] readline in function call with space in prompt.

ut the UI may be important for a function like readline() which does UI .. Martin __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] make new collumns with conditions

Dear Krissie I think you misunderstood Rui's response. He was generating some fake data to test the code not suggesting you rebuild your data frame. Michael On 25/01/2021 16:01, krissievdh wrote: Hi, Thanks for your response. I do get what you're doing. However, the table I sent

Re: [R] make new collumns with conditions

Dear Krissie I think you may be looking for the aggregate command. Note that this is a plain text list so if you post in HTML we do not see what you see. In this case we did not see any red numbers. Michael On 25/01/2021 13:57, krissievdh wrote: Hi, I have a dataset (d_vigi)with this kind

Re: [R] [SPAM] Re: Different results on running Wilcoxon Rank Sum test in R and SPSS

See comments inline On 19/01/2021 10:46, bharat rawlley wrote: Thank you for the reply and suggestion, Michael! I used dput() and this is the output I can share with you. Simply explained, I have 3 columns namely, drug_code, freq4w_n and PFD_n. Each column has 132 values (including NA). The

Re: [R] Different results on running Wilcoxon Rank Sum test in R and SPSS

posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Michael http://www.dewey.myzen.co.uk/home.html __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https:/

Re: [R] Help with connection issue for R (just joined, leading R for our agency)

Just to add to Petr's comment there are other basic editors with syntax highlighting like Notepad++ which are also OK if you want a fairly minimalist approach. Michael On 14/12/2020 08:16, PIKAL Petr wrote: Hallo Alejandra Although RStudio and ESS could help with some automation (each

Re: [R] Find the ideal cluster

Dear Jovani If you cross-post on CrossValidated as well as here it is polite to give a link so people do not answer here when someone has already answered there, or vice versa. Michael On 12/12/2020 15:27, Jovani T. de Souza wrote: So, I and some other colleagues developed a hierarchical

Re: [R] scott-knot ESD effect size test

There seems to be a package on CRAN dedicated exclusively to this test. It is called ScottKnotESD rather unoriginally. Michael. On 02/12/2020 10:32, Neha gupta wrote: I have the following data from resample svm= svm$resample$RMSE nn= nn$resample$RMSE we perform the statistical tests like

Re: [R] WG: Packegs for mathematics

be relevant. Note that the task views tell you what exists but they do not offer best buy information. For that you may need to come back here and ask about a specific subset of packages. Michael On 16/11/2020 10:02, ELLENBERGER Mihaela via R-help wrote

Re: [R] Error message when using "revtools" package

. Hopefully that will be a single entry which you can then include in your post. It also increases the chances of someone looking at it as a large .bib file is going to deter people. Michael On 05/11/2020 02:38, John wrote: Hi, I tried to read a bib file by "read_bibliography" f

Re: [R] build a literature database

Dear John If you are doing a systematic review have you thought of investigating some of the packages in the CRAN Task View on MetaAnalysis like metagear or revtools? I have not used any of them but they claim to support that part of the process. Michael On 04/11/2020 09:22, John wrote

Re: [R] Error: variable not found

Dear Hannah I think the problem is that attach() is not doing what you think it is. It does seem to make it easy to make mistakes. I would suggest switching to using with() instead or using the data = parameter to functions which support it. Michael On 30/10/2020 08:15, Hannah Van Impe

Re: [R] Need help in programming R on functional data

You have already asked this and people gave you a variety of answers. Just asking again without clarifying why those answers did not help you is not going to solve your problem. Tell us what you tried and why it failed might help. Michael On 29/10/2020 07:50, Ablaye Ngalaba wrote: Hello, I

Re: [R] Dear R experts, I have a question about meta-analysis

the archive and for instructions on registering for the list. Michael On 27/10/2020 13:32, 21803...@zju.edu.cn wrote: Dear R experts, Greetings from China! I'm Zhang in the College of Education, Zhejiang University, and I am recently running a meta-analysis. Since research using the rando

[R] vanilla session in R Gui or RStudio

1.2-18 reprex_0.3.0 assertthat_0.2.1 [55] minqa_1.2.4 rmarkdown_2.4rstudioapi_0.11 R6_2.4.1 boot_1.3-25 nnet_7.3-14 [61] nlme_3.1-149 compiler_3.6.3 > Michael Friendly Email: friendly AT yorku DOT ca Professor, Psychology Dept. & Former

Re: [R] R for mac

installing R on mac and one student is struggling to install R on Chromebook. Is there someone who would be willing to provide some pointers? Thank you for considering this request, Michael On Thu, Sep 24, 2020 at 9:10 AM Dirk Eddelbuettel wrote: > > On 23 September 2020 at 15:32, Michael Jo

Re: [R] Stats help for dissertation project

first port of call. If they fail to meet their obligations then there is a site Cross Validated where you may have better luck. Michael On 18/09/2020 18:19, Raija Hallam wrote: Hello, I am a conservation Masters student who is new to R and in need of some confirmation of my methods/ stats help

Re: [R] question including crossover trials in meta-analysis

are also likely to get more immediate help if you frame your question in terms of either meta or metafor the two packages which dmetar uses. The authors of meta and metafor all frequent the list Michael On 14/09/2020 20:41, Marc Schwartz via R-help wrote: Hi, Bert has pointed you to some R

Re: [R] truth[(truth[, 1]=="G3" & truth[, 2]=="G2") | (truth[, 1]=="G2" & truth[, 2]=="G3"), 3]<-1

I am afraid this is completely unreadable because you posted in HTML ad this is a plain text list. Best to resend it having set your mailer to send plain text as HTML gets mangled here. Michael On 06/09/2020 10:58, Hesham A. AL-bukhaiti via R-help wrote: helloout<-read.csv("outbr.cs

[R] R 4.0.2 not running on macOS 10.15.6

Greetings, I am a brand-new user to R and want to install and run this on my iMac running macOS 10.15.6. I successfully downloaded the correct package, installed it with no problems (I was not given any options for Tcl/Tk or Texinfo), and started it up as any other normal app. Some introducto

Re: [R] hist from a list

readable because you are posting in HTML and that makes the post unreadable as this is a plain text list. Michael Or barplot(bwchist) Error in barplot.default(bwchist) : 'height' must be a vector or a matrix library(PerformanceAnalytics) library(dplyr) library(tibble) librar

Re: [R] Fortune nomination .... Re: Looping through a dataframe

Dear Chris Just send it to the maintainer Achim Zeileis to whom I have cc'ed this in case he is not reading the list at the moment. Michael On 25/07/2020 06:21, Chris Evans wrote: I really don't want to put too many Emails to the list but I had the same reaction to that lovely l

Re: [R] Question about citing R packages in an academic paper

although I am not a lawyer, but it is a kindness to the authors. Michael On 22/07/2020 06:26, Dr. Mehdi Dadkhah wrote: Hi, I hope you are doing well! I have a question. I and my colleagues wrote a paper by using R language and its packages. We also used some tutorials. We have words count

Re: [R] Reading help functions

name. Of course ultimately you may prefer to read the documents in English because the program help and vignettes are usually in English and familiarity with the English terminology may help. Michael On 22/07/2020 18:20, Pedro páramo wrote: Hi all, I am trying all time to use ?? help

Re: [R] output from R to simple html

I can't tell what kind of structure you want to output, but one simple thing that I've done along these lines is to use knitr::kable to output a dataframe in either HTML or LaTeX format (for direct inclusion on a web page or to convert to PDF for other uses). -- Mike On Wed, Jul 8, 2020 at 12:02

Re: [R] how to calculate odd ratios with R?

Dear Luigi You could try the epitools package which gives a large number of ways of doing this. I would have thought that using Wald intervals for the log odds ration was not optimal with small frequencies. Michael On 06/07/2020 14:01, Luigi Marongiu wrote: Hello, Is it possible to

Re: [R] how to summarize results from studies?

Dear Frederik There is also a mailing list dedicated to meta-analysis in R https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis// Michael On 01/07/2020 16:40, Marc Schwartz via R-help wrote: Hi, It sounds like you will want to engage in a meta-analysis. There is a CRAN task view here

Re: [R] Using OpenBLAS on Ubuntu

Or maybe R-devel. On Sun, Jun 28, 2020 at 8:51 PM Jeff Newmiller wrote: > > Wouldn't this question make more sense in r-sig-debian? > > On June 28, 2020 7:54:49 PM PDT, Erin Hodgess wrote: > >Hello! > > > >Hope everyone is doing well. > > > >I'm not sure if this is the correct place to post, but

Re: [R] Error message in meta-analysis package Metafor-weights =""

in this case as the message told you. Michael On 23/06/2020 16:29, K Amoatwi wrote: Dear Wolfgang, Yes! The "metafor::forest(result.md)" works Thank you very much. Any reason why "forest(result.md)" was not working? Kobby On Tue, Jun 23, 2020 at 11:18 AM Viech

Re: [R] comparing variances/distributions

Dear Ana This really depends on your scientific question. The two techniques you have shown do different things and there must be many more which could be applied. Michael On 17/06/2020 20:57, Ana Marija wrote: Hello, I have p values from two distributions, Pold and Pnew head(m) CHR

Re: [R] Convergence in Monte Carlo Simulation

Dear Edward Every time you call your function powercrosssw() it resets the seed so you must be calling it multiple times in some way. Michael On 14/06/2020 13:57, Phat Chau wrote: Thank you Michael. I will clarify some more. The function in the first part of the code that I posted

Re: [R] Convergence in Monte Carlo Simulation

failure to converge must be meaningful about the scientific question whatever that is. At the moment your count appears to be in the correct place to count successes. Michael On 14/06/2020 02:50, Phat Chau wrote: Hello, I put together the following code and am curious about its correctness. My

Re: [R] Error in ordinal mixed effects

where you define na unless it is really NA and you did not copy it correctly. Michael On 09/06/2020 12:33, SUSANA ALBERICH MESA wrote: Hi, I'm trying to run an ordinal mixed effects model with Mixor command. I have 65 cases and repeated visits in 0, 6, 9, 12 and 18 months. My code i

Re: [R] a question of etiquette

://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Michael http://www.dewey.myzen.co.uk/home.html __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman

Re: [R] struggling with apply

pply/seq_along. sapply(seq_along(UB), function(i) pmin(UB[i], somematrix[,i])) # [,1] [,2] [,3] [,4] #[1,] 1.0 5.5 8.5 7.0 #[2,] 2.5 3.0 8.0 10.5 #[3,] 2.5 5.5 5.0 10.5 Hope this helps, Rui Barradas Às 18:46 de 27/05/20, Michael Ashton escreveu: Hi - I have a matrix of n rows and

Re: [R] struggling with apply

Always amazes me how many ways there are to do these things, none of which I was able to find myself. Thanks! I think the key here was ‘pmin,’ which I didn’t know before. Michael Ashton, CFA Managing Principal Enduring Investments LLC W: 973.457.4602 C: 551.655.8006 Schedule a Call: https

Re: [R] struggling with apply

That's it! Thanks. Learn something new every day! Michael Ashton, CFA Managing Principal Enduring Investments LLC W: 973.457.4602 C: 551.655.8006 Schedule a Call: https://calendly.com/m-ashton -Original Message- From: Rui Barradas [mailto:ruipbarra...@sapo.pt] Sent: Wednesday, M

[R] struggling with apply

[,1] [,2] [,3] [,4] [1,] 1.0 2.5 2.5 2.5 [2,] 2.5 2.5 2.5 2.5 [3,] 2.5 2.5 2.5 2.5 I'm sure there's a simple and elegant solution but I don't know what it is! Thanks in advance, Mike Michael Ashton, CFA Managing Principal Enduring Investments LLC W: 973.457.46

Re: [R] text annotation on Manhattn plot in qqman

Dear Ana That looks like something is hard coded in manhattan(). The simplest thing might be to contact the maintainer of the package and ask. You can make a copy of manhattan or textxy and modify them but I think the maintainer is the simplest course of action. Michael On 20/05/2020 16:10

Re: [R] text annotation on Manhattn plot in qqman

eters. You may also need to play with the sep =, and collapse = parameters to paste() to get the precise layout you want. Michael On 19/05/2020 17:21, Ana Marija wrote: Hi Michael, can you please send me code how that would be done? Thanks Ana On Tue, May 19, 2020 at 11:18 AM Michael Dewey

Re: [R] text annotation on Manhattn plot in qqman

Dear Ana Perhaps paste together SNP and GENE using paste() and then supply that as the snp parameter. Michael On 19/05/2020 17:12, Ana Marija wrote: Hello, I am making manhattan plot with: library(qqman) manhattan(a, chr="CHR", bp="BP", snp="SNP", p="P&q

[R] the volcano orientation

x27;t(volcano[,ncol(volcano):1])'. If you are interested in a roughly georeferenced version I have code here: https://gist.github.com/mdsumner/20fe3ffa04421bf8e0517c19085e5fd8 (Also see fortunes::fortune("conventions") ) Best, Mike -- Michael Sumner Software and Database

Re: [R] NA command in a 'for' loop

Just a thought Helen but is x being treated as a real and what you think are zero and are printed as zero are in fact some very small number? If so you need to alter your test appropriately. Michael On 20/04/2020 17:25, Helen Sawaya wrote: Thank you for your reply. I tried d[] <- lappl

Re: [R] how to get all strings in a data frame that start with a particular string

Dear Ana Would it not be possible to use grep instead of grepl and get the values using the value = TRUE parameter? Michael On 10/04/2020 17:15, Ana Marija wrote: Hello, Hello, I have a data frame (tot) with about 2000 columns. How can I extract from it all strings that start with E14? I

Re: [R] Same results but different functions ?

The documentation suggests that the rlm method for a formula does not have psi as a parameter. Perhaps try using the method for a matrix x and a vector y. Michael On 23/03/2020 12:39, varin sacha via R-help wrote: Dear R-experts, The rlm command in the MASS package command implements

Re: [R] how do I add text lables on QQ plot

Dear Ana You can specify the first three parameters to text() as vectors so it is all done in one call. That may or may not answer your question. Michael On 12/03/2020 14:08, Ana Marija wrote: HI David, thank you for getting back to me. Is there is a way for qq() to pick up text label

  1   2   3   4   5   6   7   8   9   10   >