similar to: methods: representation gives "unused arguments" on failure (PR#1799)

Full_Name: Markus Jantti Version: 1.5.1 OS: GNU/Linux Unstable (Sid) Submission from: (NULL) (193.166.59.44) I was trying to get the package yags 3.20 by Vince Carey <stvjc@channing.harvard.edu> to run, but - after successfull installation -- loading the library failed with the message > library(yags) Error in stop(paste("Duplicated slot names: ", paste(slots[duplicated(slo

i'm trying to learn R. i have no programing experience. i tried to add on a package and just picked "yags". i can't get it to work or get any help. --- > install.packages("yags") Warning in install.packages("yags") : argument 'lib' is missing: using 'C:\Documents and Settings\j\My Documents/R/win-library/2.6' trying URL

Hello R-users: I am using "yags" for fitting GEE which is giving me the same result as "Proc GENMOD". Now I have couple of questions related to yags output. (By the way, someone told me to run the geeglm for the same analysis and I did run but did not get the same result as of genmod and don't know how to correct the geeglm codes so that all three will be same!)

Hi A quick question. Standard errors reported by gee/yags differs from the ones in geeglm (geepack). require(gee) require(geepack) require(yags) mm <- gee(breaks ~ tension, id=wool, data=warpbreaks, corstr="exchangeable") mm2 <- geeglm(breaks ~ tension, id=wool, data=warpbreaks, corstr="exchangeable", std.err = "san.se") mm3 <- yags(breaks ~

I have read in a file (call it myData). The actual file is about 3000x30,000 columns and object.size() says myData takes: > 737910472/(1024^2) [1] 703.7263 Unfortunately, my program ends up using 40GB as indicated by maxvmem on Unix, which causes my department's cluster to stop working. Perhaps, I have some copying going on that I cannot find. I have created an example below that mimics

Here is an example given from ?yags library(methods) data(stackloss) Y1 <- yags(stack.loss~Air.Flow,id=1:21, data=stackloss) How can I access parts of the output. I tried: > str(Y1) Formal class 'yagsResult' [package "yags"] with 25 slots ..@ coefficients : num [1:2] -44.13 1.02 ..@ coefnames : chr(0) > Y1$coefnames Error in Y1$coefnames : $

Dear R users, I search the R archives and noted that the same problem has been posted but without solution. I know there IS instructions by the author of yags, but I just couldn't figured out. I know gee and geepack can also perform generalized estimating equation, but the reason why I need yags is I want to perform a Small Sample Adjustments by Michael P. Fay and Barry I. Graubard in

Dear R-ers, I am just a simple 'end-user' of R and am trying to analyse data with a binary response variable (dead or alive) in relation to weight and sex (of young birds). As some of the birds have the same biological mother, I am using mixed models with the identity of the mother as a random factor. (please, Mick Crawley, when are you going to write a chapter on mixed models with binary

Dear R-ers, I am just a simple 'end-user' of R and am trying to analyse data with a binary response variable (dead or alive) in relation to weight and sex (of young birds). As some of the birds have the same biological mother, I am using mixed models with the identity of the mother as a random factor. (please, Mick Crawley, when are you going to write a chapter on mixed models with binary

I have included data at the bottom of this email. It can be read in by highlighting the data and then using this command: dat <- read.table("clipboard", header = TRUE,sep="\t") I can obtain solutions with both of these: library(gee) fit.gee<-gee(score ~ chem + time, id=id, family=gaussian,corstr="exchangeable",data=dat) and library(yags) fit.yags <-

Dear list: since it took me a little while to figure out how to make use of naresid, I thought that the below R code might be useful as an example on the help page. Regards, Markus # generate some data x1 <- runif(20) y <- 10 + 5*x1 + rnorm(20) summary(lm.0 <- lm(y ~ x1)) # append some NA's to y y <- c(y, rep(NA, 5)) # generate some further x1s x1 <- c(x1, runif(5)) #

Hi Martin, John, Thanks for the responses! I've tidied up some of the notes from this mailing list thread and posted them on the bug tracker. John, in this case, I think namespaces are relevant because for non-exported S4 classes, the class information is made available through the '.__C__<package>' symbol in the package's namespace, but not the package environment that

Dear list: since it took me a little while to figure out how to make use of naresid, I thought that the below R code might be useful as an example on the help page. Regards, Markus # generate some data x1 <- runif(20) y <- 10 + 5*x1 + rnorm(20) summary(lm.0 <- lm(y ~ x1)) # append some NA's to y y <- c(y, rep(NA, 5)) # generate some further x1s x1 <- c(x1, runif(5)) #

I am having some trouble getting reftex, in particular the bibtex related features, to work properly in xemacs when editing text in Sweave files. I have added (defun Rnw-mode () (noweb-mode) (if (fboundp 'R-mode) (setq noweb-default-code-mode 'R-mode))) (add-to-list 'auto-mode-alist '("\\.Rnw\\'" . Rnw-mode)) (add-to-list 'auto-mode-alist

I am doing work on changes in establishment sizes and came across behavior that is quite understandable and easily worked around but a little surprising. On R 1.5.1 on Debian unstable (see below for R.version output): > cut.off <- c(-Inf, 0, Inf) > x <- c(-Inf, -10, 0, 10, Inf) > is.numeric(x) [1] TRUE > is.double(x) [1] TRUE > # but > cut(x, cut.off, include.lowest=T) [1]

I have a small idea as to what's going on now; at least, why exporting the class resolves this particular issue. Firstly, when an S4 class is not exported, the associated '.__C__<class>' object is not made part of the package environment. For example, I see: > getAnywhere(".__C__SubMatrix") A single object matching '.__C__SubMatrix' was found It

I need capabilities, for my data analysis, like the Pinheiro & Bates S-Plus/R package nlme() but with binomial family and logit link. I need multiple crossed, possibly interacting fixed effects (age cohort of twin when entered study, sex of twin, sampling method used to acquire twin pair, and twin zygosity), a couple of random effects other than the cluster variable, and the ability to

Hello -- I am plotting multiple densities in a single plot. If the max of a new density (added with 'lines()') is sufficiently larger than that of the one first plotted (using 'plot()'), the new density is cut or trucated off of the figure. Is there a simple way of adjusting the yscale after all lines have been added? Or is the procedure to be used to first estimate all

I am trying to use tapply to estimate the average wage in a set of Year*Education cells using sampling weights, as well as some other other statistics using the weights). I am having trouble trying to figure out how to pass the weights to the function, using: > tapply(wage,list(Educc,Year),weighted.mean, w=weight) whieh produces warnings of the type: 3: longer object length is not a

Dear All -- I am trying to use within a little table producing code an anova comparison of two gls fitted objects, contained in a list of such object, obtained using nlme function gls. The anova procedure fails to locate the second of the objects. The following code, borrowed from the help page of anova.gls, exemplifies: --------------- start example code --------------- library(nlme) ##