similar to: ncp t & Fortran error & power of some tests

Thanks Jeff, I attached a file with the program to my earlier email because the posting guide seemed to imply that non-binary attachments would work. But I see that the file was stripped off. I installed the program file on a web site, but when I downloaded it, the line breaks were stripped out. So I've included the program below: ------------------------------------------------------- #

I am trying to do power calculations for the proportional odds model using the asypow library. The code noncenta90b10<-asypow.noncent(theta.ha=a9010,info.mat=infomatrixa90b10,constraints=constrt) returns Error in max(..., na.rm = na.rm) : invalid "mode" of argument. the various arguments I've used are: a9010 [,1] [1,] -1.7357568 [2,] -0.1928619 specifying the

Per the Posting Guide, why didn't you post the reproducible R code example? On April 24, 2018 8:22:15 PM PDT, Donald Macnaughton <donmac at matstat.com> wrote: >I'm drawing a paneled histogram using the lattice package. I've >succeeded in >adding minor tick marks to the vertical axis, but I can't get the >desired >number of minor tick marks between the major

I'm drawing a paneled histogram using the lattice package. I've succeeded in adding minor tick marks to the vertical axis, but I can't get the desired number of minor tick marks between the major tick marks. I've attached a self-contained program to illustrate the problem. Thanks for your help, Don Macnaughton Here's my sessionInfo: R version 3.4.3 (2017-11-30) Platform:

Craig, Thanks for your follow-up note on using the asypow package. My problem was not only constructing the "constraints" vector but, for my particular situation (Poisson regression, two groups, sample sizes of (1081,3180), I get very different results using asypow package compared to my other (home grown) approaches. library(asypow) pois.mean<-c(0.0065,0.0003) info.pois <-

On Fri, 2020-08-21 at 01:37 +0300, Imre Deak wrote: > On Wed, Aug 19, 2020 at 05:34:15PM -0400, Lyude Paul wrote: > > (adding Ville and Imre to the cc here, they might be interested to know > > about > > this, comments down below) > > > > On Wed, 2020-08-19 at 11:15 -0400, Sean Paul wrote: > > > On Tue, Aug 11, 2020 at 04:04:50PM -0400, Lyude Paul wrote:

I would like to create a binary package for Windows which when installed will give the user an option of whether to install one version or another. Here is the problem, I have a package 'FAME' that can run in two different modes on Windows. An old mode that uses R to connect to a remote linux machine and a new way that uses a locally installed database. I would allow this decision to be

?pt states that ncp=0 is default, so this is correct > pt(3,3) [1] 0.9711656 > pt(3,3, ncp=0) [1] 0.9711656 but this is not > pt(-Inf,3) [1] 0 > pt(-Inf,3,ncp=0) [1] 0.5 is it? Torsten --please do not edit the information below-- Version: platform = i686-pc-linux-gnu arch = i686 os = linux-gnu system = i686, linux-gnu status = Under development (unstable) major = 1 minor

Does anybody has a function or idea how to calculate qt with a noncentrality parameter? I'm porting a SAS-macro for equivalence (Wellek, 1994 TT2ST) to R and it seems that I need this to calculate tinv (SAS-command, I hope that this is R's qt,...I get the SAS-macro manual next week). thanks Peter -- P.Malewski, Limmerstr.47, 30451 Hannover, 0511-2135008 At work:

Hello all, Here's my query: Running R 1.6.2 on FreeBSD 5.0, and on WinXP, and I find that the following hangs the process: dchisq(alpha=0.01, df=1, ncp=295) but it does work for ncp < 294.92. Is this general? Best wishes to all, Andrew Andrew Robinson Ph: 208 885 7115 Department of Forest Resources Fa: 208 885 6226 University of Idaho E : andrewr at uidaho.edu PO

Hi Folks, I'm looking for a neat procedure for the following: Given t0 such that pt(t0,df,ncp=0) = alpha (given) find ncp0 such that for given beta pt(t0,df,ncp=ncp0) = (1 - beta) (In other words, what's the ncp such that you get power (1-beta) to detect it, using a 1-sided test with size alpha when ncp = 0?) In the past I've done the groping by hand, but this time it needs

Dear R-users, I am wondering whether R implements a function returning the non central parameter of a t distribution (equivalent of the TnoncT function from SAS), given /x/ a value from a t distribution, /df /the degrees of freedom and /p/ the probability of x under this distribution. Thanks a lot, Matthieu -- Matthieu Dubois, /Ph.D/. /Student/ Cognitive Neuroscience Unit, UCL,

Hi everyone, I tried to use yum to search for ncpfs-2.2.1-1.i386.rpm, but it couldn't find it. I googled and found this package is included with redhat,but i'm not sure if it's included with CentOS 4.2. Any idea where I can find and install it using yum? Thanks JC

help(qt) states that: "ncp non-centrality parameter delta; currently except for rt(), only for abs(ncp) <= 37.62" so I would expect that calling qt with non-centrality parameter exceeding 37.62 should fail, instead e.g. calling > mapply(function(x) qt(p = 0.9, df = 55, ncp = x),35:45) gives: [1] 40.21448 41.35293 42.49164 43.68862 44.82945 45.97048 47.11170 48.25310 [9]

On 03/01/18 15:45, Danny Smit wrote: > Hi everyone, > > On CentOS 7 I'm running into an issue with the latest nvidia driver > from elrepo: kmod-nvidia-384.98-1.el7_4.elrepo.x86_64 > This driver version seem to introduce issue in detecting video modes > when a monitor is connected using DVI. As soon as the machine attempts > to start X, nothing happens and the monitor goes

Hi everyone, On CentOS 7 I'm running into an issue with the latest nvidia driver from elrepo: kmod-nvidia-384.98-1.el7_4.elrepo.x86_64 This driver version seem to introduce issue in detecting video modes when a monitor is connected using DVI. As soon as the machine attempts to start X, nothing happens and the monitor goes into sleep mode reporting that it has 'no signal'. It is

I've been working on this issue for quite a while. Back in January of 2015 I posted about how well external monitors worked with the out-of-the-box CentOS 7 ATI radeon drivers and my Dell M6500 with an ATI FirePro M7820. At the time the support was seamless and 'Just Worked' without any issue. A few months later, some update came down the pipe and this broke completely, with the

On Wed, Aug 19, 2020 at 05:34:15PM -0400, Lyude Paul wrote: > (adding Ville and Imre to the cc here, they might be interested to know about > this, comments down below) > > On Wed, 2020-08-19 at 11:15 -0400, Sean Paul wrote: > > On Tue, Aug 11, 2020 at 04:04:50PM -0400, Lyude Paul wrote: > > > We're going to be doing the same probing process in nouveau for > >

I have numerous replicated goodness-of-fit experiments (observed compared to expected counts in categories) and these replicates are nested within a factor. The expected counts in each cell are external (from a scientific model being tested). The calculations I need within each level of the nesting factor are a heterogeneity G test, with the total G and the pooled G across replicates. Then I

From: "Richard W.M. Jones" <rjones@redhat.com> Use the macro like this to create temporary variables which are automatically cleaned up when the scope is exited: { CLEANUP_FREE (char *, foo, strdup (bar)); /* char *foo = strdup (bar) */ ... // no need to call free (foo)! } On GCC and LLVM, this is implemented using __attribute__((cleanup(...))). On other