similar to: Proper / Improper scoring Rules

Hi All, Could someone help me decode what this error means ? > BIC(nb.80) Error in log(attr(object, "nobs")) : Non-numeric argument to mathematical function > BTW, nb.80 is a negative binomial glm model created using the MASS library with the call at the bottom of the message In the hopes of trying to figure this out I tried the following workaround but it did not work

Dear R listers, I am learning the MLE utility optim() in R to program ordered logit models just as an exercise. See below I have three independent variables, x1, x2, and x3. Y is coded as ordinal from 1 to 4. Y is not yet a factor variable here. The ordered logit model satisfies the parallel regression assumption. The following codes can run through, but results were totally different from what I

Dear All, I have a pair of .cpp and .def file can be compiled using VC++ and works perfectly well in S-PLUS. I wanted to do the same for R; so I followed the guidline given in "Building from a source-code library under Windows" as much as possible and manage to compile them using VC++ and call it from R. But it gives different answer from the one called from S-Plus. I know that I did

Dear Friends, I have postfix mail server When I try to install amavisd-new, clamAV, SpamAssassin follow this link "http://wiki.centos.org/HowTos/Amavisd". I have done all the configuration after that when I restart clamav So It is giving my this error message Stopping Clam AntiVirus Daemon: [FAILED] Starting Clam AntiVirus Daemon: LibClamAV Error:

Earlier this month I reported a problem in that ClamTK couldn't virus-scan a file. This was cured by following the dev's instructions <quote> If you open up a terminal window and type "ls /var/clamav" (without quotes of course), I'm betting you'll see a variety of files and/or directories in there. If you're up for it, as root type rm /var/clamav/* -rf which

Dear list members, I'm running some tests looking at differences between means for various levels of a factor, using Tukey's HSD method. I would like to plot the data as boxplots or dotplots, with horizontal significance lines indicating which groups are statistically significantly different, according to Tukey HSD. Here's a nice image showing an example of such a graphical

Hello, I'm estimating an ordered logit model on a probability weighted survey sample. polr permits case weights with the "weights" option, but I cannot figure out from existing documentation what it actually does with these weights. I'm concerned about this because I get somewhat different results using Stata's ologit command with the pweights option and very

Hi, ?? Are there any packages which prints out p-values for OLR's (like `ologit' from Stata)? I want to run a bunch of OLRs and print the p-value for the first coefficient from each of them. ? I checked polr() under MASS and it doesn't. ?There's a lrm() function under Design which does print out p-values but I couldn't extract p-values from the output. ? Thanks, ? Debs

>>>>> Winston Chang <winstonchang1 at gmail.com> >>>>> on Fri, 23 Jan 2015 10:15:53 -0600 writes: > I think you can simplify a little by replacing this: > pkg %in% loadedNamespaces() > with this: > .getNamespace(pkg) almost: It would be !is.null(.getNamespace(pkg)) > Whereas getNamespace(pkg) will load the

All- Thank you in advance for any help you might be able to lend. Here is my issue. I am trying to open a fairly large .dat file. The file originally was downloaded as a GZ file but I unzipped it (with 7-zip) into it's current 1.86 gig .dat format. I know that the data is "just a plain ASCII file with 720 columns and 360 rows per time step (month). It should be readable by

I fixed the country code issue but that did not resolve the problem. I also removed all files in /var/lib/clamav and reran freshclam (without rebooting), that also did not fix the problem. Jay > Am 15.07.2018 um 00:13 schrieb Jay Hart: >> ClamAV update process started at Sat Jul 14 15:10:48 2018 >> Using IPv6 aware code >> Querying current.cvd.clamav.net >> TTL: 1232

I want to test whether the proportional odds assumption for an ordered regression is met. The UCLA website points out that there is no mathematical way to test the proportional odds assumption (http://www.ats.ucla.edu/stat//R/dae/ologit.htm), and use graphical inspection ("We were unable to locate a facility in R to perform any of the tests commonly used to test the parallel slopes

(Apologies if this is better suited for R-help.) On my system (macOS Sierra, late 2014 MacBook Pro; R 3.4.1, Homebrew build), I found that it is faster to construct a function using eval(call("function", ...)) than using as.function(list(...)). Example: make_fn_1 <- function(a, b) eval(call("function", a, b), env = parent.frame()) make_fn_2 <- function(a, b)

>>>>> Michael Lawrence <lawrence.michael at gene.com> >>>>> on Mon, 26 Jan 2015 05:12:55 -0800 writes: > A isNamespaceLoaded() function would be a useful thing to > have in general if we are interested in readable code. An > efficient implementation would be just a bonus. Good point (readability), and thank you for the support! Note

I decided to do a direct comparison of transpose and sweep. library(microbenchmark) NN <- matrix(c(1, 2, 3, 4, 5, 6), nrow = 2, byrow = TRUE) # Example matrix lambda <- c(2, 3, 4) # Example vector colNN <- t(NN) microbenchmark( sweep = sweep(NN, 2, lambda, "/"), transpose = t(t(NN)/lambda), colNN = colNN/lambda ) Unit: nanoseconds expr min lq

Hello, I get this error when compiling this code with llvm-gcc: error: unsupported inline asm: input constraint with a matching output constraint of incompatible type! #include <stdio.h> #include <stdlib.h> #include <sys/time.h> #include <sys/types.h> #include <unistd.h> int main(void) {

Hi, I am running an analysis with an ordinal outcome and I need to run a test of the parallel regression assumption to determine if ordinal logistic regression is appropriate. I cannot find a function to conduct such a test. >From searching various message boards I have seen a few useRs ask this same question without a definitive answer - and I came across a thread that indicated there is no

These macros will be used by multiple VMWare modules for handling host communication. v2: * Keeping only the minimal common platform defines * added vmware_platform() check function v3: * Added new field to handle different hypervisor magic values Signed-off-by: Sinclair Yeh <syeh at vmware.com> Reviewed-by: Thomas Hellstrom <thellstrom at vmware.com> Reviewed-by: Alok N Kataria

These macros will be used by multiple VMWare modules for handling host communication. v2: * Keeping only the minimal common platform defines * added vmware_platform() check function v3: * Added new field to handle different hypervisor magic values Signed-off-by: Sinclair Yeh <syeh at vmware.com> Reviewed-by: Thomas Hellstrom <thellstrom at vmware.com> Reviewed-by: Alok N Kataria

IVR Number :17275691533 When I try it from xlite configuring my provider directly, it works perfectly. When I try to dial out from dialer , it doesnt work. [sip8] type=peer username=user fromuser=user authuser=user secret=password host=8.14.146.111 nat=no canreinvite=yes insecure=very disallow=all allow=g729 allow=ulaw context=default dtmfmode=rfc2833 What cld be the reason ? --------------