Code in constrOptim() (R-2.0.0, package stats): if (any(ui %*% theta - ci <= 0)) stop("initial value not feasible")
but the help page tells us in Section "Details": "The feasible region is defined by ui %*% theta - ci >= 0." Uwe Ligges --please do not edit the information below-- Version: platform = i386-pc-mingw32 arch = i386 os = mingw32 system = i386, mingw32 status = major = 2 minor = 0.0 year = 2004 month = 10 day = 04 language = R Windows NT 4.0 (build 1381) Service Pack 6 Search Path: .GlobalEnv, package:methods, package:stats, package:graphics, package:grDevices, package:datasets, package:utils, package:fortunes, Autoloads, package:base ______________________________________________ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-devel