search for: nrang

Dear all, I wanted to compare Bonferroni vs TukeyHSD correction over a range of groups and group sizes, and wanted to use the function qtukey. In the help documentation it says qtukey(p, nmeans, df, nranges = 1, lower.tail = TRUE, log.p = FALSE) Arguments q vector of quantiles. p vector of probabilities. nmeans sample size for range (same for each group). df degrees of freedom for s (see below). nranges number of groups whose maximum range is considered. log.p logical; if TRUE, probabilit...

...ution of the maximum of r statistics each distributed as the Studentized range of means calculated from c random samples of size n from normal populations. The rc samples are assumed to be mutually independent and a common pooled?within?samplevariance is used throughout." So the connection is nranges == r, and nmeans == c. (n never actually factors in because sqrt(n) is part of the standardization) For the typical application, r is 1 for the usual studentized range distribution. E.g. for two large groups: > qtukey(.95,2,df=Inf) [1] 2.771808 As there is only one difference to consider, t...

...u_notifier); struct hmm_mirror *mirror; - struct hmm_update update; struct hmm_range *range; unsigned long flags; int ret = 0; @@ -173,15 +172,10 @@ static int hmm_invalidate_range_start(struct mmu_notifier *mn, if (!kref_get_unless_zero(&hmm->kref)) return 0; - update.start = nrange->start; - update.end = nrange->end; - update.event = HMM_UPDATE_INVALIDATE; - update.blockable = mmu_notifier_range_blockable(nrange); - spin_lock_irqsave(&hmm->ranges_lock, flags); hmm->notifiers++; list_for_each_entry(range, &hmm->ranges, list) { - if (update.end &l...

...AM PDT, Ursula Garczarek <Ursula.Garczarek at cytel.com> wrote: >Dear all, >I wanted to compare Bonferroni vs TukeyHSD correction over a range of >groups and group sizes, and wanted to use the function qtukey. > >In the help documentation it says > >qtukey(p, nmeans, df, nranges = 1, lower.tail = TRUE, log.p = FALSE) >Arguments >q > >vector of quantiles. > >p > >vector of probabilities. > >nmeans > >sample size for range (same for each group). > >df > >degrees of freedom for s (see below). > >nranges > >number o...

...ole and I am also afraid of making any mistakes while typing the codes, I would like to make this code simple by using for loop. I know that the assign function is useful and another code as follows worked well with it: for(i in 1:N_var){ assign(paste("V", i, sep=""), post_df[Nrange$start[i]:Nrange$end[i],]) } However, I could not find how to apply for loop to 1) - 3) above. 1) >for (i in 1:N_var){ + assign(colnames(paste("V", i, sep=""))[2], par$V2[i]) + } Error in assign(colnames(paste("V", i, sep = ""))[2], par$V2[i]) : inva...

...ne comment on a related cleanup: > > > list_for_each_entry(mirror, &hmm->mirrors, list) { > > int rc; > > > > - rc = mirror->ops->sync_cpu_device_pagetables(mirror, &update); > > + rc = mirror->ops->sync_cpu_device_pagetables(mirror, nrange); > > if (rc) { > > - if (WARN_ON(update.blockable || rc != -EAGAIN)) > > + if (WARN_ON(mmu_notifier_range_blockable(nrange) || > > + rc != -EAGAIN)) > > continue; > > ret = -EAGAIN; > > break; > > This magic handling of er...

...;hmm->ranges, list) { - if (range->valid) - continue; - range->valid = true; - } - wake_up_all(&hmm->wq); - } - spin_unlock_irqrestore(&hmm->ranges_lock, flags); -} - -static int hmm_invalidate_range_start(struct mmu_notifier *mn, - const struct mmu_notifier_range *nrange) -{ - struct hmm *hmm = container_of(mn, struct hmm, mmu_notifier); - struct hmm_mirror *mirror; - struct hmm_range *range; - unsigned long flags; - int ret = 0; - - spin_lock_irqsave(&hmm->ranges_lock, flags); - hmm->notifiers++; - list_for_each_entry(range, &hmm->ranges, list) {...

...llwig <hch at lst.de> One comment on a related cleanup: > list_for_each_entry(mirror, &hmm->mirrors, list) { > int rc; > > - rc = mirror->ops->sync_cpu_device_pagetables(mirror, &update); > + rc = mirror->ops->sync_cpu_device_pagetables(mirror, nrange); > if (rc) { > - if (WARN_ON(update.blockable || rc != -EAGAIN)) > + if (WARN_ON(mmu_notifier_range_blockable(nrange) || > + rc != -EAGAIN)) > continue; > ret = -EAGAIN; > break; This magic handling of error seems odd. I think we should merge rc an...

...p: > > > > > list_for_each_entry(mirror, &hmm->mirrors, list) { > > > int rc; > > > > > > - rc = mirror->ops->sync_cpu_device_pagetables(mirror, &update); > > > + rc = mirror->ops->sync_cpu_device_pagetables(mirror, nrange); > > > if (rc) { > > > - if (WARN_ON(update.blockable || rc != -EAGAIN)) > > > + if (WARN_ON(mmu_notifier_range_blockable(nrange) || > > > + rc != -EAGAIN)) > > > continue; > > > ret = -EAGAIN; > > > break; &...

...(res2) doesn't work and reports NULL res2$df.residual doesn't work and reports NULL Is this an undocumented feature or an error? It sure plays havoc with Dr. Bates' contribution of Tukey.aov which works with res0 and res1 but doesn't work with res2 and reports Error in qtukey(p, nranges, nmeans, df, lower.tail, log.p) : Non-numeric argument to mathematical function A simple but inelegant workaround is to input mse and df.residual as function arguments in addition to the aov object. Does anyone have a better solution? Thank you, Peter B. -- Peter B. Mandeville...

...gt; > list_for_each_entry(mirror, &hmm->mirrors, list) { > > > > int rc; > > > > > > > > - rc = mirror->ops->sync_cpu_device_pagetables(mirror, &update); > > > > + rc = mirror->ops->sync_cpu_device_pagetables(mirror, nrange); > > > > if (rc) { > > > > - if (WARN_ON(update.blockable || rc != -EAGAIN)) > > > > + if (WARN_ON(mmu_notifier_range_blockable(nrange) || > > > > + rc != -EAGAIN)) > > > > continue; > > > > ret = -EAGAI...

Here are seven more patches for things I found to clean up. This was based on top of Christoph's seven patches: "hmm_range_fault related fixes and legacy API removal v3". I assume this will go into Jason's tree since there will likely be more HMM changes in this cycle. Changes from v1 to v2: Added AMD GPU to hmm_update removal. Added 2 patches from Christoph. Added 2 patches as

Dear all, I'm using R to perform multiple comparison testing on agriculture genotype trials. To perform the Duncan's MRT, I use the qtukey() function with the following syntax: qtukey(p = ((1 - 0.05) ^ (pos - 1)), nmeans = pos, df = ni) I experience a strange behaviour when the number of means in the trial and the number of residual degrees of freedom (ni) becomes high

Hi Dan, Jérôme and Jason, below is a series that cleans up the dev_pagemap interface so that it is more easily usable, which removes the need to wrap it in hmm and thus allowing to kill a lot of code Note: this series is on top of Linux 5.2-rc6 and has some minor conflicts with the hmm tree that are easy to resolve. Diffstat summary: 34 files changed, 379 insertions(+), 1016 deletions(-) Git

From: Jason Gunthorpe <jgg at mellanox.com> 8 of the mmu_notifier using drivers (i915_gem, radeon_mn, umem_odp, hfi1, scif_dma, vhost, gntdev, hmm) drivers are using a common pattern where they only use invalidate_range_start/end and immediately check the invalidating range against some driver data structure to tell if the driver is interested. Half of them use an interval_tree, the others

From: Jason Gunthorpe <jgg at mellanox.com> 8 of the mmu_notifier using drivers (i915_gem, radeon_mn, umem_odp, hfi1, scif_dma, vhost, gntdev, hmm) drivers are using a common pattern where they only use invalidate_range_start/end and immediately check the invalidating range against some driver data structure to tell if the driver is interested. Half of them use an interval_tree, the others