[R] lmekin() function in kinship package

Hi, I had a problem with the lmekin() in kinship package: lmekin() can not be wrapped into another function library(kinship) #creat an example dataset xx-rnorm(100) yy-rnorm(100) id-1:100 test.dat-as.data.frame(cbind(xx,yy,id)) rm(xx,yy,id)

[R] unbiased sandwich variance estimator for GEE estimates

Hi, Anyone knows any existing package/program that has implemented unbiased (or bias-reduced) sandwich variance estimator (Wu (1986) and Carroll (2001) for GEE estimates? Thanks Qiong __ R-help@stat.math.ethz.ch mailing list

[R] plot groupedData in nlme

Hi, Does anyone know how to make the color of the lines all black when plotting groupedData with an outer factor: For example, library(nlme) plot(Dialyzer, outer=~QB, key=F) This generated colored curves in R.2.4.1. How to make all the curves black ? (or how to alter the color (type) of lines

[R] How to exclude a level from a factor

: factor(x,exclude=factor(c,levels=c(a,b,c))) [1] a b c Levels: a b c What's wrong with my codes? Thanks for any help Qiong Yang __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] How to exclude a level from a factor

Problem solved. Thanks a lot for your replies! x [1] a b c Levels: a b c factor(as.character(x),exclude=c) [1] abNA Levels: a b exclude= option may not work on factors. One has to convert the factor to character first. Qiong On Tue, 20 Sep 2005, Sundar Dorai-Raj wrote: Qiong