Re: [R] Metafor -can't calculate heterogeneity with non-positive sampling variances

2014-08-27 Thread Owen, Branwen
much difference. I'm just beginning to wonder if there is another way. I'll think about it some more best wishes Branwen From: Viechtbauer Wolfgang (STAT) [wolfgang.viechtba...@maastrichtuniversity.nl] Sent: 27 August 2014 17:30 To: Owen, Bran

[R] Metafor -can't calculate heterogeneity with non-positive sampling variances

2014-08-27 Thread Owen, Branwen
e sampling variances. Is there any way around this? thanks Branwen From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] on behalf of r-help-ow...@r-project.org [r-help-ow...@r-project.org] Sent: 27 August 2014 13:36 To: Owen, Branwen Subject: M

Re: [R] how to make R plot under Linux

2014-05-14 Thread Owen
Vojtěch Zeisek opensuse.org> writes: > > On Linux You can use Rkward http://rkward.sourceforge.net/ - very nice and > good graphical user interface for R. > > Dne Pá 19. února 2010 19:39:53 xinwei stat.psu.edu napsal(a): > > Hi, I am using R in Linux environment. How can i make plot in Linux

[R] Moran's I

2013-07-18 Thread Owen Yancey
Dear list, I’m an undergraduate research assistant trying to create a single function that can calculate Moran’s I using the “ape” package. The example that I will use comes from online data about ozone levels. Here are the commands that I have been using: > library(ape) > ozone<-read.table("ht

[R] metafor matrix error

2013-05-21 Thread Branwen Owen
hi I've been using this code to set a reference level for uni and multivariate analysis rmix$interviewmethodcode<-relevel(mix$interviewmethodcode,ref="SAQ") summary(rma.1<-rma(yi,vi,mods=~interviewmethodcode,data=rmix,method="SJ",knha=F,weighted=F,intercept=T)) giving

[R] metafor - multivariate analysis

2013-03-15 Thread Owen, Branwen
Dear Metafor users, I'm conducting a metaanalysis of prevalence of a particular behaviour based on someone elses' code. I've been labouring under the impression that this: summary(rma.1<-rma(yi,vi,mods=cbind(approxmeanage,interviewmethodcode),data=mal,method="DL",knha=F,weighted=F,intercept=T

Re: [R] plot the y-axis upside down

2011-06-17 Thread Owen, Jason
That's fantastic! Thanks. > -Original Message- > From: Marc Schwartz [mailto:marc_schwa...@me.com] > Sent: Friday, June 17, 2011 3:47 PM > To: Owen, Jason > Cc: 'r-help@r-project.org' > Subject: Re: [R] plot the y-axis upside down > > On Jun

[R] plot the y-axis upside down

2011-06-17 Thread Owen, Jason
Hello, I need to create a scatterplot where the y-axis is upside down. If I have non-negative bivariate data in objects "x" and "y," then the operation plot(x, -y) gives me the figure I want -- a mirror image of plot(x, y) -- except that the y-values (coordinates) are negative, which I don't

[R] write.dta and export of variable labels

2011-01-15 Thread Andrew Owen
When I import a data file in stata format using read.dta, the full variable labels (e.g. variable= "inc2000"; variable label= "Total household income in year 2000"). When I then export the same data file using write.dta, the variable labels are not included in the new .dta file. Is there any way to

[R] Write.dta and export of variable labels

2011-01-14 Thread Andrew Owen
When I import a data file in stata format using read.dta, the full variable labels (e.g. variable= "inc2000"; variable label= "Total household income in year 2000"). When I then export the same data file using write.dta, the variable labels are not included in the new .dta file. Is there any way to

[R] Multiple Strata Sampling in R Survey

2010-07-27 Thread Owen Gallagher
But that really doesn't meet the needs of this project. I've also tried it with strata=(~strataOne + strataTwo), but that didn't work either. Does anyone know if its a) possible to do mutiple strata sampling in R survey, and b) how? I emailed the maintainer a week ago, but I still

Re: [R] problem with "dynformula" from "plm" package [RE-POST]

2009-11-27 Thread Owen Powell
, logs or diffs in my model. Any other ideas? ~Owen R> rm(list = ls()) R> options(prompt= "R> ") R> library("plm") R> data("EmplUK", package="plm") R> EmplUK <- plm.data(EmplUK, index = c("firm", "year")) R> log(

Re: [R] problem with "dynformula" from "plm" package [RE-POST]

2009-11-27 Thread Owen Powell
> plm(formula = formula(a), data=EmplUK) rather than just: R> plm(formula = a, data = EmplUK) ~Owen 2009/11/27 David Winsemius : > You might also note that the authors hint on the help page that one might > want to use the formula() operation on the result of dynformula. Following > that

[R] problem with "dynformula" from "plm" package [RE-POST]

2009-11-27 Thread Owen Powell
indeed, when I try to run a regression using that formula, it appears to not contain any lags or logs (output below). Any ideas? Thanks in advance, ~Owen -- Owen Powell http://center.uvt.nl/phd_stud/powell R> library("plm") R> data("EmplUK", package="plm") R&g

[R] problem with "dynformula" from "plm" package

2009-11-27 Thread Owen Powell
indeed, when I try to run a regression using that formula, it appears to not contain any lags or logs (output below). Any ideas? Thanks in advance, ~Owen -- Owen Powell http://center.uvt.nl/phd_stud/powell R> library("plm") R> data("EmplUK", package="plm") R&g

Re: [R] joint estimation of two poisson equations

2009-04-13 Thread Owen Powell
~ x1 + x2 out = lm.pb(eq1, eq2, data = d, zeroL3 = TRUE) print(out) I couldn't find out how to get standard errors and p-values from the package, so I bootstrapped them. ~Owen 2009/4/13 Tirthankar Chakravarty > You should probably try the -bivpois- package: > http://cran.r-projec

[R] joint estimation of two poisson equations

2009-04-13 Thread Owen Powell
that unless it's entirely necessary. I'll post my solution to the list when I've worked it out. Regards, ~Owen # CODE FOR "sur" OPTION rm(list = ls()) library(Zelig) y1 = c(1,2,3,4) y2 = c(0,2,0,2) x = c(2,3,4,8) d = data.frame(cbind(y1, y2, x)) eq1 = y1 ~ x eq2 = y2 ~

[R] Computing/Interpreting Odds Ratios for 3-way interactions from lmer

2009-01-01 Thread Owen, Amanda J
at there is improvement in the Cond when, second clause (as compared to condwhen-first clause-age-matched, right?). Again, my trouble is with what to use as the reference group and how to separate that out in terms of computing percentages. Thanks so much for any assistance you can

[R] Fill in NA values in vector with previous character/factor

2008-07-28 Thread Owen Jones
. should evaluate to: "A" "B" "B" "B" "C" "C" "C" "C" "C" "D" "D" "D" I tried to use rle() in a function to do this but have hit a brick wall. How would YOU do this? Many tha

Re: [R] strange timings in convolve(x,y,type="open")

2007-12-18 Thread owen
go quickly. -Art Quoting "Charles C. Berry" <[EMAIL PROTECTED]>: > On Tue, 18 Dec 2007, Art Owen wrote: > >> Dear R-ophiles, >> >> I've found something very odd when I apply convolve >> to ever larger vectors. Here is an example below >>

[R] strange timings in convolve(x,y,type="open")

2007-12-18 Thread Art Owen
like N or N log(N) , with the difference hard to discern in granularity and noise. The convolve function is not very big (see below). When type is not specified, it defaults to "circular". So my guess is that something mysterious might be happening inside the first else clause below,

Re: [R] Observations on SVD linpack errors, and a workaround

2007-10-17 Thread Art Owen
> Ph: (410) 502-2619 > > Fax: (410) 614-9625 > > Email: [EMAIL PROTECTED] > > Webpage: http://www.jhsph.edu/agingandhealth/People/Faculty/Varadhan.html > > > > ---- > > > > -Original Message- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On > Behalf Of Art

[R] Observations on SVD linpack errors, and a workaround

2007-10-16 Thread Art Owen
write this, it seems to be giving firefox a very bad time loading it. So proceed with caution (if at all) to the file badx in the web page stat.stanford.edu/~owen/ There is also a smaller one, called badx2 which illustrates the much rarer case where a skinny matrix makes svd choke, while its wide