[R] empty cell when running GEE for binary data

some empty cells, so I can't simply use logit as a link function. I was wondering if you know any existing R solution for this? also, I tired to use "cauchit" as the link function, but it turned out that within the function geeglm "cauchit" is not available. Any idea on

Re: [R] empty cell when running GEE for binary data

p(A) would be 1/0, so it's not appropriate to use logit as the link function. I don't think, but not sure, gee would work, but i will try the other package you suggested! Thanks! Yue -- View this message in context: http://r.789695.n4.nabble.com/empty-cell-when-running-GEE-for-bina

[R] separation occuring in clustered data

andle such a problem using a penalised-likelihood method. However, it seems that logistf is only suitable when observations are completely independent, thus will be problematic if used on my clustered data. I was wondering if anyone had similar experience and has any suggestion for this? Thanks very

[R] exact logistic regression with correlated data?

hi all, does anyone know if R can do exact logistic regression with correlated binary data? Thanks!! Yue -- View this message in context: http://r.789695.n4.nabble.com/exact-logistic-regression-with-correlated-data-tp4635900.html Sent from the R help mailing list archive at Nabble.com

[R] firth's penalized likelihood bias reduction approach

tion. So, when involving the interaction term as a predictor, I also need to control for within-subject correlation. Has anyone experience a similar problem and how you solved it? or, any suggestion would be very much appreciated!!! Thanks very much!! Best, Yue -- View this message in context: http://r

[R] [R-pkgs] interplot: a new package for visualizing interactive effects

bout the package are available in http://cran.r-project.org/web/packages/interplot/vignettes/interplot-vignette.html. Please contact me with any questions, bug reports, and comments. Best, Hu, Yue Ph.D. Student, Political Science, University of Iowa, 313 Schaeffer Hall, 20E Washington St.,

[R] ffsave.image() error with large objects

Hi, I have been using ffsave.image() to save mixture of ff and normal objects in my workspace. e.g. ffsave.image(file = "C:\output\saveobjects", rootpath = "D:\fftempdir", safe = TRUE) It works fine but once my workspace has large (~4GB) objects, I get the error: Error in ffsave.image(file = "C:

Re: [R] ffsave.image() error with large objects

So it seems that I've hit the 4GB zipfile limit. However, even if I set compress = FALSE, it still tries to compress, and hence fails. YS On 2/11/10, Yue Sheng wrote: > Hi, I have been using ffsave.image() to save mixture of ff and normal > objects in my workspace. e.g. > >

[R] Save workspace with ff objects

individual ff objects separately? I've tried the naive way of just saving the workspace, only to find that ff objects are empty. Thanks for help. Yue Sheng [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] bigglm "update" with ff

Hi, since bigglm doesn't have update, I was wondering how to achieve something like (similar to the example in ff package manual using biglm): first <- TRUE ffrowapply ({ if (first) { first <- FALSE fit <- bigglm(eqn, as.data.frame(bigdata[i1:i2,,drop=FALSE]), chunksize = 1, family =

[R] Failed path analysis using sem package

Dear R users: I am trying to run a path analysis using sem package in R. But I have encountered one problem, below is my code: SEMEX<-read.csv("D:/Documents and Settings/z3409964/Desktop/Hospital 1.csv") library(sem) cov.matrixSEMEX<-cov(na.omit(SEMEX)) SEMEX<-specifyModel() CWB->CWB13,NA,1 CWB

Re: [R] using metafor for meta-analysis of before-after studies

if it is possible to calculate the ¡®sdi¡¯ (the standard deviation of the change scores) using the SD of percentage change. Thank you very much, I am looking forward to your reply. With best wishes. Qiang Yue M.D. Visiting Scholar of IMHR, University of Ottawa 1145 Carling Avenue, K1Z

Re: [R] using metafor for meta-analysis of before-after studies

regards. Qiang Yue From: Michael Dewey Date: 2013-04-27 07:28 To: qiangmoon; wvb CC: r-help Subject: Re: [R] using metafor for meta-analysis of before-after studies At 03:27 27/04/2013, Qiang Yue wrote: >Hello, Dr. Viechtbauer. > >I am trying to perform a meta-analyis on a group >of

Re: [R] using metafor for meta-analysis of before-after studies

Hello, Michael. Thanks very much, your suggestion is really helpful, and hopefully will address our issue. I will also read the suggested references to learn more about meta-analysis. With best wishes. Qiang Yue From: Michael Dewey Date: 2013-04-28 06:30 To: qiangmoon; Michael Dewey; wvb

Re: [R] using metafor for meta-analysis of before-after studies (escalc, SMCC)

), we can not use the percentage change to calculate SMCC, but have to get the raw change first? With best wishes. Qiang Yue From: Viechtbauer Wolfgang (STAT) Date: 2013-05-21 10:09 To: Moon Qiang; r-help Subject: RE: [R] using metafor for meta-analysis of before-after studies (escalc, SMCC

Re: [R] using metafor for meta-analysis of before-after studies (escalc, SMCC)

Dear Dr. Viechtbauer: Thanks so much! Now all these issues are clear. With best regards. Qiang Yue From: Viechtbauer Wolfgang (STAT) Date: 2013-05-23 05:06 To: qiangmoon; r-help Subject: RE: RE: [R] using metafor for meta-analysis of before-after studies (escalc, SMCC) The mean

[R] Fontconfig error in R 3.10 when running

default config file dev.off() The plot displays distorted text (attached), and I just couldn't figure out where the fix is after google for hours. Any helps will be very much appreciated. Many thanks, Yue __ R-help@r-project.org mailing list

[R] Fontconfig error in R 3.10 when running

default config file dev.off() The plot displays distorted text (attached), and I just couldn't figure out where the fix is after google for hours. Any helps will be very much appreciated. Many thanks, Yue __ R-help@r-project.org mailing list

[R] source file for rnorm

function? Thanks. Yue Zhang, CFA Cohen & Steers Capital Management, Inc. 280 Park Ave., 10th Floor New York, NY 10017 Tel: (212)796-9370 Fax: (212)319-8238 Email: yzh...@cohenandsteers.com The information transmitted is only for the intended recipient and may con

[R] How to let Rterm be vi style under windows?

Hi list, As title, under windows, the methord of ~/.inputrc won't work, dunno how to let Rterm can be vi style? -- Regards, Yue Wu Key Laboratory of Modern Chinese Medicines Department of Traditional Chinese Medicine China Pharmaceutical University No.24, Tongjia Xiang Street, Na

[R] Dealing with a lot of parameters in a function

e sum of 24 of these negative log-likelihood. Thanks. Yue Notice: This e-mail message, together with any attachme...{{dropped:11}} __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www

[R] Error msg in optim : (list) object cannot be coerced to type 'double'

e parameters. There also was a short dicussion about it a few years ago on this mailing list, but no final resolution was reached for it. Yue Notice: This e-mail message, together with any attachme...{{dropped:14}} __ R-help@r-project.org mailing lis

[R] How to let Rterm be vi style under windows?

Hi list, As title, under windows, the method of ~/.inputrc won't work, is it possible to let Rterm be vi style on windows? P.S. I've sent it serveral days ago, get no response, so I send it again to check if I've failed to send it. -- Regards, Yue Wu Key Laboratory of

[R] [w32] Bug? When trying to auto-complete, rest of the input from cursor is wiped out.

te the code, but the rest is gone: > write.table(file='foo',fileEncoding| No this issue when on *nix. -- Regards, Yue Wu Key Laboratory of Modern Chinese Medicines Department of Traditional Chinese Medicine China Pharmaceutical University No.2

[R] bigglm binomial negative fitted value

..$ : NULL Anyone can help on this case? Thank you in advance. Best, --Yue __ R-help@r-project.org mailing list 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.

Re: [R] bigglm binomial negative fitted value

Thank you very much. I do overlook something. On Thu, May 31, 2012 at 5:20 PM, Thomas Lumley wrote: > On Fri, Jun 1, 2012 at 1:17 AM, Yue Guan wrote: >> Hi, there >> >> Since glm cannot handle factors very well. I try to use bigglm like this: >> >> logit_mo

[R] Correlated Multivariate Distribution Generator

0.8 | | 0.8 0.8 1 0.9 | | 0.8 0.8 0.9 1 | Thank a lot. Best, Yue Yu __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

Re: [R] Correlated Multivariate Distribution Generator

be changed after transformation, won't keep the desired parameters r and p. Enrico's method is brilliant for rank correlation, is there any algorithm for linear correlation? Even some references for multivariate negative binomial generator will help. Your help are really appreciated Best,

[R] how to find the minimum

re to find b or p" Would you please give some suggestions? Thanks in advance. Best wishes, Bin Yue - Best regards, Bin Yue * student for a Master program in South Botanical Garden , CAS -- View this message in context: http://www.nabble.com/how-to-find-t

[R] Ｗhat is the "deviance " and the "ch ange in deviance" of a logistic regression?

ed due to missingness) AIC: 3618.5 Number of Fisher Scoring iterations: 5 > deviance(glm1) [1] 3612.516 is the "deviance(glm1)" the "deviance" ? If yes, what is the "change in deviance"? Is it "residual deviance-null.deviance"? Thanks all . Best wishes, Bin

[R] logistic regression using "glm",which "y" is set to be "1"

" . I don't know which level of the response was set to be 1. is the first element of the response set to be 1? Thank all in advance. Regards, - Best regards, Bin Yue * student for a Master program in South Botanical Garden , CAS -- View this message in context: ht

Re: [R] logistic regression using "glm",which "y" is set to be "1"

;. Thank you very much . Best regards, Bin Yue Marc Schwartz wrote: > > > On Wed, 2007-12-05 at 18:06 -0800, Bin Yue wrote: >> Dear friends : >> using the "glm" function and setting family=binomial, I got a list of >> coefficients. >>

Re: [R] logistic regression using "glm",which "y" is set to be "1"

oes anybody know how to obtain the transformation result ? It is hard to settle down before knowing the actual process R works . I have read some books and the "?glm" help file , but what they told me was not sufficient. Best wishes , Bin Yue Weiwei Shi wrote: > > Dear Bin:

[R] what actually a logistic regression fits

ix(small.glm)[,i]->res[,i] + } > apply(res,1,sum)->right > plot(left,right) #I got a straight line whose slope was about 1 Therefore ,I am almost sure that what a logistic regression fits is log(p/(1-p))~b0+b1*x. Is that right? Regards, - Best regards, Bin Yue *

[R] paradox about the degree of freedom in a logistic regression model

set affects the result of regression ,such as the degree of freedom. Dose anybody have any idea about this? Thank all who read this message. Regards, Bin Yue - Best regards, Bin Yue * student for a Master program in South Botanical Garden , CAS -- View this message in conte

[R] the observed "log odds" in logistic regression

nder how the residual is computed. Would you please help me ? Thank all very much again. Regards, Bin Yue - Best regards, Bin Yue * student for a Master program in South Botanical Garden , CAS -- View this message in context: http://www.nabble.com/the-observed-%2

Re: [R] the observed "log odds" in logistic regression

Dieter Menne: Thank you for your reply! I know that I don't have to do any logit , but I want to understand how R fit the glm models. I will read the examples your suggested . Best regards, Bin Yue Dieter Menne wrote: > > Bin Yue 163.com> writes: > >> After re

[R] re store pictures in an automatically named files using loops

env[[i]],main=names(spl.sp)[i]) legend("topleft",lty=1:4,col=1:4,legend=c("obs","theo","hi","lo")) text(2,500,cex=0.8, paste("n=",spl.sp[[i]]$n,sep="")) bmp(paste(names(spl.sp)[i],".bmp",sep="")) } } }

Re: [R] re store pictures in an automatically named files using loops

I put bmp() before the plot() syntax , then the code worked. Bin Yue wrote: > > Dear all: > I hope that the file in which the picture is stored has the same name > as the " main" title of the picture . But it turns out that , the name of > the file is :names