Message-ID: <AS8PR08MB91930282442674DAD517F2C38B332@AS8PR08MB9193.eurprd08.prod.outlook.com>
Date: 2024-03-20T22:34:16Z
From: Wolfgang Viechtbauer
Subject: [R-meta] Calculation of p values in selmodel
In-Reply-To: <23ca01da7b12$19b808d0$4d281a70$@gmail.com>
See: https://www.metafor-project.org/doku.php/installation#development_version
> -----Original Message-----
> From: R-sig-meta-analysis <r-sig-meta-analysis-bounces at r-project.org> On Behalf
> Of Will Hopkins via R-sig-meta-analysis
> Sent: Wednesday, March 20, 2024 23:01
> To: 'R Special Interest Group for Meta-Analysis' <r-sig-meta-analysis at r-
> project.org>
> Cc: Will Hopkins <willthekiwi at gmail.com>
> Subject: Re: [R-meta] Calculation of p values in selmodel
>
> Wolfgang, I figured that I might be able to implement your selmodel(...,
> pval=dat$pValue) by reverting to for-loop processing. I solved the previous
> problems I had with this approach by stripping out the NULLs from the list
> objects, so the simulations are running nicely.
>
> However, it appears that the latest version of metafor, which I have just
> installed, does not have the pval functionality yet, because I get this warning
> with every invocation of selmodel:
> Extra argument ('pval') disregarded.
>
> I am still a newbie with R, so I may have done something wrong with the update
> process. I got this warning, and nine others like it, but I presume this doesn't
> apply to metafor:
> Warning: package 'boot' in library 'C:/Program Files/R/R-4.3.2/library' will not
> be updated.
>
> I then got this, again presumably irrelevant:
> Warning: cannot remove prior installation of package ?lattice?
> Warning: restored ?lattice?
>
> The only warning for loading metafor was this:
> package ?Matrix? was built under R version 4.3.3
>
> The other possibility is that my for-loop processing has some kind of mistake,
> but it seems to be OK. Each loop produces a meta-analysis object resultmeta from
> invocation of ram.uni with a dataset meta_data (the i-th subset of my full
> dataset). I then apply selmodel to resultmeta within a function, where the x is
> resultmeta:
> result <- metafor::selmodel(x, type="step", steps=(0.025),
> pval=meta_data$pValue)
>
> I give up. Help, please!
>
> Will