search for: lrt

configure uses clock_gettime to check whether -lrt is needed - and don''t check other functions. With glibc 2.17 clock_gettime is part of libc, so use timer_gettime instead, which is in -lrt in old and new versions of glibc. Signed-off-by: Olaf Hering <olaf@aepfle.de> --- configure | 2 +- 1 file changed, 1 insertion(+), 1 deletion(...

HI, Suppose I have the following data structure. LRT tp 1 1.50654010 522 2 0.51793929 522 3 0.90340299 522 4 1.20293325 522 5 1.05578774 523 6 0.01617942 523 7 0.68183543 523 8 0.43820244 523 9 1.14123995 524 10 0.05809550 524 11 0.93061597 524 12 1.39739700 524 13 1.05220953 525 14 0.03471461 525 15 0.63168798 525 16 1.4059...

Dear List, having searched the help and archives, I have the impression that there is no automatic model selection procedure implemented in R that includes/excludes predictors in logistic regression models based on LRT P-values. Is that true, or is someone aware of an appropriate function somewhere in a custom package? Even if automatic model selection and LRT might not be the most appropriate methods, I actually would like to use these in order to simulate someone else's modeling approach... Many thanks fo...

Hi there, This is not a strictly related R question; I apologize for this. I'm interest in simulate the sampling distribution of the LRT, testing for trend among the levels of some categorical variable X in a regression model. (in practice this is achieved by assigning scores to the levels of X and fitting such numeric variable). To simulate the sampling distribution the steps are: 1. Simulate the model *without* trend 2. for every...

I would like to try a likelihood ratio test in place of waldtest. Ideally I'd like to provide two glm models, the second a submodel of the first, in the style of lrt (http://www.pik-potsdam.de/~hrust/tools/farismahelp/lrt.html). [lrt takes farimsa objects] Does anyone know of such a likelihood ratio test? Chris Elsaesser, PhD Principal Scientist, Machine Learning SPADAC Inc. 7921 Jones Branch Dr. Suite 600 McLean, VA 22102 703.371.7301 (m) 703.6...

...> Makes sense. > > Acked-by: Matt Wilson <msw@amazon.com> > > > On Wed, Jan 16, Olaf Hering wrote: > > > > > Ping > > > > > > On Wed, Dec 05, Olaf Hering wrote: > > > > > > > configure uses clock_gettime to check whether -lrt is needed - and don''t > > > > check other functions. With glibc 2.17 clock_gettime is part of libc, so > > > > use timer_gettime instead, which is in -lrt in old and new versions of > > > > glibc. > > > > > > > > Signed-off-by: Ol...

On my setup with glibc-2.8, benchmark_residual linkage fails with undefines references to clock_gettime(). Adding -lrt fixes that. The following is a small patch for it. Regards. -- O.S. diff --git a/configure.ac b/configure.ac index 993ac33..392485e 100644 --- a/configure.ac +++ b/configure.ac @@ -437,6 +437,11 @@ if test x$enable_stack_smash_protection = "xyes" ; then XIPH_GXX_STACK_PROTECTOR...

My CentOS box with: $ uname -r 2.6.18-164.6.1.el5 gets the file order wrong with: ls -lrt * . Two files that were last modified in the same second are shown in order that is the reverse of the order in which they were actually last modified. Note also that the later file was first created after the earlier file was last modified. Is the precision of this only to the second? Mike.

Dear list, I am fitting a logistic multi-level regression model and need to test the difference between the ordinary logistic regression from a glm() fit and the mixed effects fit from glmer(), basically I want to do a likelihood ratio test between the two fits. The data are like this: My outcome is a (1,0) for health status, I have several (1,0) dummy variables RURAL, SMOKE, DRINK, EMPLOYED,

...ff -r f35234b6636a -r 7e326d27d899 tools/blktap2/drivers/Makefile --- a/tools/blktap2/drivers/Makefile Tue Mar 15 16:33:59 2011 +0000 +++ b/tools/blktap2/drivers/Makefile Wed Mar 16 10:38:52 2011 +0000 @@ -23,11 +23,7 @@ CFLAGS += -fPIC CFLAGS += -fPIC endif -LIBS += -lrt -lz - -LBLIBS_img := $(LDLIBS_libxenctrl) $(CRYPT_LIB) -lpthread -lz -lm - -LIBS += -L$(LIBVHDDIR) -lvhd +VHDLIBS := -L$(LIBVHDDIR) -lvhd REMUS-OBJS := block-remus.o REMUS-OBJS += hashtable.o @@ -88,26 +84,26 @@ all: $(IBIN) lock-util qcow-util all: $(IBIN) lock-util qcow-util -tapdis...

...ing to this value of Chi. I looked up in the current litterature but cannot really find a proper answer to that dilmena. So, dear evolutionary biologists R users, how would you present it ? Thank you very much, Julien. -- View this message in context: http://old.nabble.com/How-shall-one-present-LRT-test-statistic-in-a-scientific-journal---tp26532480p26532480.html Sent from the R help mailing list archive at Nabble.com.

...-Wmissing-declarations -Wpointer-arith -Wchar-subscripts -Wformat=2 -Wbad-function-cast -fno-builtin-strftime -Wstrict-aliasing=2 -Wl,--as-needed -o .libs/anvil main.o anvil-connection.o anvil-settings.o connect-limit.o penalty.o -Wl,--export-dynamic ../../src/lib-dovecot/.libs/libdovecot.so -ldl -lrt -Wl,-rpath -Wl,/usr/local/lib/dovecot ../../src/lib-dovecot/.libs/libdovecot.so: undefined reference to `cdb_read' ../../src/lib-dovecot/.libs/libdovecot.so: undefined reference to `cdb_find' ../../src/lib-dovecot/.libs/libdovecot.so: undefined reference to `cdb_free' ../../src/lib-dove...

...ove the DrugPair term and model4 in which I only keep the fixed effect of MatingPair. 3. I want to perform the log likelihood ratio test to compare these models and that's why I have REML=F. However the code anova(model1, model2, model3, model4) gives me a chisq estimate and a p-value, not the LRT values. How do I get LRT (L.Ratio) while using lmer? 4. I am under the impression after reading a few posts that LRT is not usually obtained with lmer but it is given if I use lme (the old model). 5. I could not find how to input the random interaction term while using lme? Is it the following wa...

Hi all -- I have to calculate a likelihood ratio test for a logit model. I found logLik, but I need to calculate the log likelihood for the model without any predictors. How can I specify this in glm? If the full model is glm(y ~ x1), is the one without predictors (y ~ 0)? Or (y ~ 1)? Is there a more direct way of getting this? -- Chris

....2l.32/lib ./.libs/libhttp.a > ../lib-dns/.libs/libdns.a ../lib-ssl-iostream/.libs/libssl_iostream.a > ../lib-master/.libs/libmaster.a ../lib-settings/.libs/libsettings.a > ../lib-test/.libs/libtest.a ../lib/.libs/liblib.a > ../lib-ssl-iostream/.libs/libssl_iostream_openssl.so -lcrypto -lrt -lnsl > -lsocket -lsendfile -lcrypt -R/Applications/dovecot-2.2.31.32/lib/dovecot I simulated your compile arguments as best as I could, and compared the libtool command in my compile logs with your's. (mine) Mine: -lssl -lcrypto -lrt -lnsl -lsocket -lsendfile (your's) -lcrypt...

...x1 represents a 'baseline' covariate, whereas x2 represents a 'new' covariate, and my goal is to figure out where x2 adds significant predictive information over x1. Ideally, I could get a p-value for doing this. Originally, I thought of doing some kind of likelihood ratio test (LRT), where i measure the (partial) likelihood of the model with just x1, then with x1 and x2, then it becomes a LRT with 1 degree of freedom. But when i use the summary() function for coxph(), i get the following output (shown at the bottom). I have two questions: 1) What exactly are the p-values i...

...ith -Wchar-subscripts -Wformat=2 >> -Wbad-function-cast -fno-builtin-strftime -Wstrict-aliasing=2 >> -Wl,--as-needed -o .libs/anvil main.o anvil-connection.o anvil-settings.o >> connect-limit.o penalty.o -Wl,--export-dynamic >> ../../src/lib-dovecot/.libs/libdovecot.so -ldl -lrt -Wl,-rpath >> -Wl,/usr/local/lib/dovecot >> ../../src/lib-dovecot/.libs/libdovecot.so: undefined reference to >> `cdb_read' >> ../../src/lib-dovecot/.libs/libdovecot.so: undefined reference to >> `cdb_find' >> ../../src/lib-dovecot/.libs/libdovecot.so: un...

...list(Eth, Sex, Lrn, Age), sum)))) First fit a saturated model and see which term we can drop: > fm <- glm(Count ~ .^4, quine2, family=poisson) > drop1(fm, test="Chisq") Single term deletions Model: Count ~ .^4 Df Deviance AIC LRT Pr(Chi) <none> 6.661e-16 148.916 Eth:Sex:Lrn:Age 2 0.422 145.338 0.422 0.81 Drop the four-way interaction and check which term we can drop next: > fm <- update(fm, . ~ . - Eth:Sex:Lrn:Age) > drop1(fm, test="Chisq")...

In running the following example of a loglinear model for the Titanic data, I was surprised to see NaN reported for the Pearson chisq > loglin(Titanic, margin=list(1:3, 4)) 2 iterations: deviation 2.273737e-13 $lrt [1] 671.9622 $pearson [1] NaN $df [1] 15 $margin $margin[[1]] [1] "Class" "Sex" "Age" $margin[[2]] [1] "Survived" Tracing it back, this occurs because there are zeros in the fitted/expected frequencies for children among the Crew. > # get fitted...

...eem >> to like it still: > > > It looks like the jit doesn't find these because they are located in > librt. Try this (or adapt to whatever system you're using): > > lli -load /usr/lib/librt.so yourprog.bc oh. thanks that worked. > > If you pass '-lrt' when linking your program, it should take care of > this for you. > ah ok. so every library thet gccld can't find as a bytecode lib is added to the shell script then. -- Jakob