search for: rsela

I am writing a program in which I would like to take in a formula, change the response (Y) variable into something else, and then pass the formula, with the new Y variable to another function. That is, I am starting with formula <- Y~X1+X2+X3 and I'd like to do something like Y <- formula$Y newY <- f(Y) lm(newY~X1+X2+X3) So far, it seems that my

When I use a model fit with LME, I get an error if I try to use "predict" with a dataset consisting of a single line. For example, using this data: > simpledata Y t D ID 1 -1.464740870 1 0 1 2 1.222911373 2 0 1 3 -0.605996798 3 0 1 4 0.155692707 4 0 1 5 3.849619772 1 0 2 6 4.289213902 2 0 2 7 2.369407737 3 0 2 8 2.249052533 4 0 2 9 0.920044316 1

I am trying to estimate models with subsets using the NLME package. However, I am getting an error in the case below (among others): > subset <- c(rep(TRUE, 107), FALSE) > fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1, subset=subset) Error in xj[i] : invalid subscript type 'closure' > fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1,

I am in the process of submitting a package to CRAN. R CMD check ran successfully on the package on my local computer, using R version 2.1.1. However, on the computers for CRAN (with version 2.10.0), the following errors occurred: Warning in parse_Rd("./man/predict.Rd", encoding = "unknown") : ./man/predict.Rd:28: unknown macro '\dquote' *** error on file

I am creating an R package. I ran R CMD check on the package, and everything passed until it tried to run the examples. Then, the result was: * checking examples ... ERROR Running examples in REEMtree-Ex.R failed. The error most likely occurred in: > ### * AutoCorrelationLRtest > > flush(stderr()); flush(stdout()) > > ### Name: AutoCorrelationLRtest > ### Title: Test for

...reshape/cast more than once to get to my answers so I haven't been against learning how to work with it.    Overall they are good inputs and I appreciate them. Thanks! Cheers, Mark ------------------------------ Message: 80 Date: Wed, 8 Jul 2009 14:31:49 -0400 (EDT) From: Rebecca Sela <rsela@stern.nyu.edu> Subject: [R] \dQuote in packages To: r-help <r-help@r-project.org> Message-ID:     <31158951.1855241247077909779.JavaMail.root@calliope.stern.nyu.edu> Content-Type: text/plain; charset=utf-8 I am in the process of submitting a package to CRAN.  R CMD check ran success...

...ies shall be understood as neither given nor endorsed by it. ______________________________________ The Shaw Group Inc. http://www.shawgrp.com      [[alternative HTML version deleted]] ------------------------------ Message: 76 Date: Wed, 27 May 2009 13:36:41 -0400 (EDT) From: Rebecca Sela <rsela@stern.nyu.edu> Subject: [R] "Error: package/namespace load failed" To: r-help <r-help@r-project.org> Message-ID:     <24291478.7471501243445801339.JavaMail.root@calliope.stern.nyu.edu> Content-Type: text/plain; charset="utf-8" I am writing my first R package, and...

I want to take the first row of each unique ID value from a data frame. For instance > ddTable <- data.frame(Id=c(1,1,2,2),name=c("Paul","Joe","Bob","Larry")) I want a dataset that is Id Name 1 Paul 2 Bob > unique(ddTable) Will give me all 4 rows, and > unique(ddTable$Id) Will give me c(1,2), but not accompanied by the name column.

I am trying to optimize a likelihood function using NLMINB. After running without a problem for quite a few iterations (enough that my intermediate output extends further than I can scroll back), it tries a vector of parameter values NaN. This has happened with multiple Monte Carlo datasets, and a few different (but very similar) likelihood functions. (They are complicated, but I can send them

According to the R documentation for NLMINB, the returned value of convergence is 0 for successful convergence. When I got another code (1), I looked up the PDF that linked from the documentation (http://netlib.bell-labs.com/cm/cs/cstr/153.pdf), which said that a return code under 3 was impossible. Is there other documentation that gives the correct meanings of the NLMINB convergence codes in

I just built a package in R (version 2.1.1 on a Linux machine) and then moved it to my computer (R version 2.8.1 on a Windows machine). When I loaded the package and tried to open help for one of the commands, I got the following warning: Warning message: In print.help_files_with_topic("F:/R/R-2.8.1/library/REEMtree/chm/RMSE") : No CHM help for 'RMSE' in package

I fit a random effects linear model to data, and then tried to use it to predict, but I got this error: > predict(lmeObject, newdata, level=0) Error in eval(mCall$fixed)[-2] : object is not subsettable This is a new error for me. It still occurs if I change the level to 1 or if I change the data for prediction back to the original dataset to which the lme model was fitted. Can anyone help?

I tried to use LME (on a fairly large dataset, so I am not including it), and I got this error message: Error in lme.formula(formula(paste(c(toString(TargetName), "as.factor(nodeInd)"), : nlminb problem, convergence error code = 1 message = false convergence (8) Is there any way to get more information or to get the potentially wrong estimates from LME? (Also, the page in the

I have two straightforward questions about linking in the man pages for R packages: First, is it possible to link from within parts of the documentation that are not the \seealso section? For example, I would like to have something like: \arguments{ \item{correlation}{an optional \code{corStruct} object describing the within-group correlation structure; the available classes are given in

I am writing my first R package, and I have been getting the following series of errors when I run R CMD check: * checking S3 generic/method consistency ... WARNING Error: package/namespace load failed for 'REEMtree' Call sequence: 2: stop(gettextf("package/namespace load failed for '%s'", libraryPkgName(package)), call. = FALSE, domain = NA) 1: library(package,