similar to: parameter-restrictions in OPTIM

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

Hello What is the appropriate method to test if a given distribution is a weibull thank you meriema email meriema.belaidouni at int-evry.fr

I wonder if something like this has been reported before: > pexp(85:86,0.438) [1] 1 -Inf --please do not edit the information below-- Version: platform = i386-pc-mingw32 arch = i386 os = mingw32 system = i386, mingw32 status = major = 1 minor = 5.0 year = 2002 month = 04 day = 29 language = R Windows NT 4.0 (build 1381) Service Pack 6 Search Path: .GlobalEnv,

Hi all, Suppose I have a vector like this: [1] "STAT1" "STAT1" "STAT1" "STAT1" "GAPDH" "GAPDH" "GAPDH" "ACTB" "ACTB" [10] "ACTB" "DDR1" "RFC2" "HSPA6" "PAX8" "GUCA1A" "UBE1L" "THRA" "PTPN21" [19]

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 ***

Dear "Debuggers", please have look a the following: > plot(0:1,0:1,type="n") > points(c(0,0.5,1),c(0,0.5,1)) > dev.print(device = postscript, width=5 , height=5, file="test1.ps") windows 2 ># This worked > dev.print(device = pdf, width=5 , height=5, file="test1.pdf") windows 2 ># That also >

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]

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've got a Solaris 2.6 SPARC system running Samba 2.0.6 and an NT-4.0sp5 PDC system. If Samba is running on the Sun system when the NT PDC is rebooted it (the NT system) complains that there is already another PDC in the domain. This should not be as far as I can see as I've tried specifically to configure the Samba systems not to do this and in fact prior to upgrading to Samba version

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

Hi Timo, when logging out like a001 logout the imap child dies from signal 11. The back trace looks like this: Program received signal SIGSEGV, Segmentation fault. 0xb7ed4991 in strcasecmp () from /lib/tls/i686/cmov/libc.so.6 (gdb) bt #0 0xb7ed4991 in strcasecmp () from /lib/tls/i686/cmov/libc.so.6 #1 0x0806ab6c in command_unregister (name=0x815b9ab "LOGOUT") at commands.c:83 #2

Hi there colleagues, atapci3: <nVidia nForce MCP55 SATA300 controller> port 0xbc00-0xbc07,0xb880-0xb883,0xb800-0xb807,0xb480-0xb483,0xb400-0xb40f mem 0xefcb3000-0xefcb3fff irq 23 at device 5.2 on pci0 atacontrol detach ata7 - insert ATA disk (ad14) atacontrol attach ata7 pinics with Fatal trap 12: page fault while in kernel mode (kgdb) bt #0 doadump () at pcpu.h:196 #1 0xc0533227

Dear experts, when trying to estimate an kernel density function with density(x) I get the following error message with imported data from either EXCEL or text files: Error in density(spr) : argument must be numeric. Other procedues such as truehist work. If I generate data within R density works fine. Does anybody have an idea? Yours -- Christoph Hanck Wissenschaftliche Hilfskraft

Kevin writes: > Of course, I can fix it by myself and check every terms length, but > that will add more overhead to big data computing. How is the overhead different whether your code checks it or Xapian does? Best regards, Adam -- "Oh, we all like motorcycles, to some degree." Adam Sj?gren asjo

Hello, I am looking for a documentation or an example to use the synonym function. I tried this db.add_synonym("omega","xapain"); and this works by adding the flag FLAG_AUTO_SYNONYMS. If i try to use the db.add_synonym("omega","xapain is search engine "); it fails why? Can xapian use synonym for Abbreviations like MBA => Master of business

Hi there, I have a xapian index whose results are being sorted by a value, with (PHP bindings): $enquire->set_sort_by_value($sort_data_value); This is because I want the results returned in chronological order of publication date. However, I now have a need to have certain results be 'sticky' at the top of the resultset, regardless of their publication date. Obviously there are

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))

Hi, I've been working on a patch for dovecot 1.2 from the Kolab branch (http://hg.intevation.org/kolab/dovecot-1.2_kolab-branch/) that implements listing of shared namespaces. I've got something that works in some basic way but is still missing some pieces. See the attached patch, which also contains some installation and configuration notes. Implementation notes: One of the main

Dear R users/experts, I've heard it's possible to estimate the above kinds of models in R. However, after (an admittedly brief) survey of the packages, I haven't found an obvious candidate. Can you offer any help? (Yes I'm a newby.) Yours sincerely Christoph Hanck Studentische Hilkskraft Lehrstuhl f?r Empirische Wirtschaftsforschung, Prof. Dr. Wilfling