similar to: Classic fail-safe N

I would suggest to post this to the (recently created) R-sig-meta-analysis mailing list. See: https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis Best, Wolfgang >-----Original Message----- >From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Naike Wang >Sent: Monday, June 26, 2017 08:32 >To: R-help at r-project.org >Subject: [R] Classic fail-safe N >

Hi all, I was conducting a meta-analysis of single proportions(i.e. without a control group) using the metafor package. When I performed a classic fail-safe N, I noticed that the result (the number of missing studies that would bring p-value to the alpha, to be exact)was different than that I got in Comprehensive Meta-Analysis Version 2.0. I wonder why R and CMA got different results. *Below is

Dear users of metafor, I am working on a meta-analysis using the metafor package. I have a excel csv database that I am working with. I am interested in pooling the effect measures for a particular subgroup (European women) in this csv database. I am conducting both sub-group and meta-regression. In subgroup-analyses, I have stratified the database to create a separate csv file just for European

Hi Dr. Viechtbauer, I'm doing meta-analysis using your package 'metafor'. I used the 'IRFT' to transform the incident rate. But when I tried to back-transform the summary estimates from function rma, I don't know what's the appropriate ti to feed in function transf.iirft. I searched and found your post about using harmonic mean for ni to back-transform the double

Hi, I am a student using R for my final year project, with the metafor package being particularly helpful. I have been following the steps laid out in the manual 'Conducting Meta-analysis in R with the Metafor package' (Viechtbauer, 2010) and applying it to my own data of infected Anopheles mosquitoes across Africa. In particular, I am attempting to apply a mixed-effect model to my data

Dear R community, I'm using the ,metafor'-package by Wolfgang Viechtbauer (Version: 0.5-5) to calculate random-effects meta-analyses using Correlations and Sample Sizes as the raw data. (By the way: Really a nice piece of work, Wolfgang! Thanks heaps.) I specified the "rma.uni' function so that it looks like this: MAergebnis<-rma.uni(ri=PosOutc, ni=N,

Dear metafor users, I am working on a meta-analysis of reliability and the correlation associations. I need some help about conducting categorical moderators variables. Questions 1: How to conduct the weighted ANOVAs assuming a mixed-effects model on the tranformed alpha coefficients/the tranformes correlation coefficients for the categorical moderator variables? Questions 2: How to

Hi all, I am not a very proficient R-user yet, so I hope I am not wasting people?s time. I want to run a linear mixed model with 3 random factors (A, B, C) where A and B are partially crossed and C is nested within B. I understand that this is not easily possible using lme but it might be using lmer. I encountered two problems when trying: Firstly, I can enter two random factors in lmer but

Hello, I am performing a meta-analysis using the metafor package. My data are proportions and I used the Freeman Tukey double arcine (FT) transformation to fit the random effects model. Now I want to create a forest plot with my estimates backtransformed to the original scale of proportions. Can this be done? Regards, Patricia

Bonjour, J'ai un petit souci avec mon petit couple rsync over ssh. Quand je fais un rsync d'une machine, et que celle-ci dispose d'une messagerie de type vpopmail, le backup m'indique ? la fin que certains fichiers n'ont pu ?tre transf?r?s. Ceci est du au fait que sachant que les services POP/SMTP sont actif, des fichiers ont ?t? supprim?s entre le temps ou la liste des

Dear All, I'm having trouble tweaking a forest plot made using the R meta-analysis package metafor. I did the analysis based upon the correlation coeff from studies and plotted the corresponding forest plot easily > q2<-rma(yi,vi,mods=cbind(grupo),data=qim) > q2 > forest (q2,transf=transf.ztor,digits=3, ... ,alim=c(0,1),refline=.5) > text(-1.55,42,"Esp?cie

Dear R experts, I am new to the list and R software. I need to convert arcview file to Winbugs having R has middle package. Got from friends that it is possible following the steps: 1. Converting arcview shapefile to "cgm clear text file" 2. Downloading "convert.r" to into R 2.1 use source("convert.r") 2.2 convert("filename_cgm") without the

Hello all, I am trying to do a Latin Hypercube Sampling (LHS) to a 5-parameter design matrix. I start as follows: library(lhs) p1<-randomLHS(1000, 5) If I check the distribution of each parameter (column), they are perfectly uniformly distributed (as expected).For example, hist(p1[,1]) Now the hard (maybe strange) question. I want the combination of the first three parameters to sum up to

I thought the curve function was a very flexible way to draw functions. So I could plot funtions like the following: # I created a function to produce functions, for instance: fp <- function(m,b) function(x) sin(x) + m*x + b # So I can produce a function like this ff <- fp(-0.08, 0.2) ff(1.5) # Is the same as executing sin(1.5) - 0.08*1.5 + 0.2 # Let's plot this

Hi All, I am trying to code an R script which gives me the time varying parameters of the NIG and GH distributions. Further, becasue I think these these time varying parameters should be more responsive to more recent observations, I would like to include a weighted likelihood estimation proceedure where the observations have an exponentially decaying weighting rather than the equal weighting

Dear Helpers, I just started working with R and I'm a bit overloaded with information. My data is from marsupials reindroduced in a area. I have weight(wt), hind foot lenghts(pes) as continues variables and origin and gender as categorial. condition is just the residuals i took from the model. > names(dat1) [1] "wt" "pes" "origin" "gender"

Hi all, I want to create two groups of random numbers to calculate proportions. The first group is to represent the number of cases in a study. The second group is to represent the sample size of the study. Apparently, the sample size is going to have to be bigger or equal to the number of cases, but the sample size of a study is not necessarily greater than the number of cases of another study.

Hello Helpers, I have some problems with fitting the model for my data... -->my Literatur says (crawley testbook)= Non-normality of errors-->I get a banana shape Q-Q plot with opening of banana downwards Structure of data: origin wt pes gender 1 wild 5.35 147.0 male 2 wild 5.90 148.0 male 3 wild 6.00 156.0 male 4 wild 7.50 157.0 male 5 wild 5.90

I have a vector Pes with 600 elements some of which are NA's. How can I form a vector of the indices of the NA's. for(i in 1:600) if(is.na(Pes[i])) print(i) prints the indices of the NA's but I can't figure out how to put the results in a vector. Thank you very much. Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de

Hi all, I just joigned this list, so I'm sorry if this problem has already been reported. I'm running Dovecot 2.0.13 on many servers, one for POP/IMAP access, others for LDA, others for authentification only, etc. All servers are accessing a shared file system, based on MooseFS (www.moosefs.org). The FS is mounted using FUSE. All my Dovecot servers have this configuration :