similar to: Error using nls()

I'm solving 4 complex equations simultaneously. Code is below. The code returns only zero's for the solution though there should also be a non-zero result. I'm pretty confident that the equations are correct because they are straight from a published paper and I checked them pretty thoroughly. The parameter values I used are from the published paper as well. Any suggestions for how

Hello, I see in BugTraq that there's yet another problem with Wu-ftpd, but I see no mention of it in the freebsd-security mailing list archives...I have searched the indexes from all of June and July. Wu is pretty widely used, so I'm surprised that nobody seems to have mentioned this problem in this forum. The notice on BugTraq mentioned only Linux, not FreeBSD, but that's no

Hi, All I'm brand new to the list, and here's what's going on: I am running samba 1.9.18 on a RedHat Linux box (RH release 5.2). Samba was installed via RPM rather than built locally. While configuring smbd for a different purpose than I write about here, I discoverd smbtar. It struck me that this might prove a useful tool for backing up an NT mail-server (icky, I know) from a

Douglas Bates wrote: > > The GLMM function in the lme4 package allows you to specify crossed > random effects within the random argument without the need for the > pdBlocked and pdIdent constructions. Simply ensure that your grouping > factors are defined in such a way that each distinct group has a > different level in the grouping factor (this is usually not a problem

Dear helpeRs, a colleague of mine would like to give R a try. He uses econometric models which typically involve a large number of variables, esp. time series. Having no experience with handling very large data sets myself I turn to you. 1. Could you please describe your experiences to cope with these situations? 2. What kind of difficulties will he have to face? Are there special

Dear HelpeRs, I'm very fond of Sweave and I use it as often as possible. It'a a pity I can't use it for larger projects or can I? For instance suppose I have three files file1.rnw, file2.rnw and file3.rnw with Sweave code. Working on file2.rnw I whould like to exclude file1.rnw and file3.rnw temporarily and joining all of them later. This amounts to a mechanism similar to using

Dear HelpeRs, I wonder if anyone knows of ways to read EViews file types. I did not find a function in the package 'foreign' and a search query submitted to http://search.r-project.org was not successful. Any hint is very much welcome. Dietrich Trenkler -- Dietrich Trenkler c/o Universitaet Osnabrueck Rolandstr. 8; D-49069 Osnabrueck, Germany email: Dietrich.Trenkler at

I found the following strange behavior using qnorm() and pnorm(): > x<-8.21;x-qnorm(pnorm(x)) [1] 0.0004638484 > x<-8.22;x-qnorm(pnorm(x)) [1] 0.01046385 > x<-8.23;x-qnorm(pnorm(x)) [1] 0.02046385 > x<-8.24;x-qnorm(pnorm(x)) [1] 0.03046385 > x<-8.25;x-qnorm(pnorm(x)) [1] 0.04046385 > x<-8.26;x-qnorm(pnorm(x)) [1] 0.05046385 > x<-8.27;x-qnorm(pnorm(x))

I have some data which I would like to display with dotchart. The labels are very long, so the chart becomes too small. Setting cex=0.7 seems to be a good compromise, but the gap between the text and the chart still is too large. I did not find a "gap" parameter in the description of dotchart... Thanks for any help. D. Trenkler "a" <- structure(c(103.35, 36.73, 55.09,

Niels Waller wrote: > Does anyone know of an R (or S-PLUS) function for delineating all possible > combinations and permutations? The following function delivers all permutations of 1:n. all.perm <- function(n) { p <- matrix(1, ncol = 1) for (i in 2:n) { p <- pp <- cbind(p, i) v <- c(1:i, 1:(i - 1)) for (j in 2:i) { v <- v[-1]

I think I've discovered what went wrong. My workspace included a function wilcox.test formerly copied from ctest. Now ctest ist part of the distribution and because of that I always got the message: [Previously saved workspace restored] Error in autoload("wilcox.test", "ctest") : Object already exists I didn't take care of this message. Now I removed

Hallo everyone, excuse me if this is not a genuine R question but I do not know where to ask else. Referring to e.g. https://stat.ethz.ch/pipermail/r-help/2004-December/062114.html I wonder if these measurements of 3000 criminals (raw data) are available anywhere. At least I didn't find them in the R datasets package or by means of Google. What I did find was a table of frequencies of

I suspect that there is bug in sub when using "?": > string_"This is a bug!" > sub("!", ", or isn't it?", string) [1] "This is a bug, or isn't it?" > string_"This is a bug?" > sub("?", ", or isn't it?", string) [1] "This is a bug?" Regards, *** D.Trenkler ***

The following is a data frame > "jjd" <- structure(list(Observations = c(6.8, 6.6, 5.3, 6.1, 7.5, 7.4, 7.2, 6.5, 7.8, 9.1, 8.8, 9.1), LevelA = structure(c(1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3), .Label = c("A1", "A2", "A3"), class = "factor"), LevelB = structure(c(1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2), .Label =

Dear R-community, not to re-invent the wheel I wonder if someone of you has ever written a function to compute the GHK smooth recursive simulator to estimate multivariate normal probabilities. See for instance page 194 of @BOOK{Greene97, author = {William H. Greene}, year = 1997, title = {Econometric Analysis}, edition = {3rd}, publisher = {Prentice-Hall}, address = {New Jersey

Dear HelpeRs, I have some data: "ice" <- structure(c(0.386, 0.374, 0.393, 0.425, 0.406, 0.344, 0.327, 0.288, 0.269, 0.256, 0.286, 0.298, 0.329, 0.318, 0.381, 0.381, 0.47, 0.443, 0.386, 0.342, 0.319, 0.307, 0.284, 0.326, 0.309, 0.359, 0.376, 0.416, 0.437, 0.548, 41, 56, 63, 68, 69, 65, 61, 47, 32, 24, 28, 26, 32, 40, 55, 63, 72, 72, 67, 60, 44, 40, 32, 27, 28, 33,

Is this a bug? >"xx" <- c(50, 79, 120, 78, 90, 100, 78, 80, + 90, 80, 60, 39, 90, 85, 140, 100, 80, 80) > boxplot.stats(xx)$stats [1] 60 78 80 90 100 > boxplot.stats(0.1*xx)$stats [1] 7.8 7.8 8.0 9.0 10.0 I suppose the way the numbers in 0.1*xx are internally stored may cause this problem: > formatC(0.1*xx[2],format="f",digits=16) [1]

Hello all, I wonder if there is a package including a program to display a 3dim histogram. Thanks. --- D.Trenkler --- **************************************************************************** ***** Dr. Dietrich Trenkler (dtrenkler at nts6.oec.uni-osnabrueck.de) Statistik / Empirische Wirtschaftsforschung Universitaet Osnabrueck Rolandstrasse 8

I have two lists of equal lengths. Each one consists of vectors. Is there an elegant way to combine them elementwise? Elegant means: without loops. For instance, for > jja_list(1:3,4:7) > jjb_list(rep(0,4),(2,5)) I would like to get: > list(c(1:3,rep(0,4)),c(4:7,rep(2,5))). --- D.Trenkler --- ************************************************************************ *********

(Moved from r-help) On 3/13/2006 9:33 AM, Dietrich Trenkler wrote: > I have some data which I would like to display with dotchart. The > labels are very long, so the chart becomes too small. Setting cex=0.7 > seems to be a good compromise, but the gap between the text and the > chart still is too large. I did not find a "gap" parameter in the > description of dotchart...