similar to: statistics

FWIW the link from ?dbinom to the Loader paper on Binomials is broken but the paper seems to be available here: https://octave.1599824.n4.nabble.com/attachment/3829107/0/loader2000Fast.pdf Roger Koenker r.koenker at ucl.ac.uk<mailto:r.koenker at ucl.ac.uk> Honorary Professor of Economics Department of Economics, UCL Emeritus Professor of Economics and Statistics, UIUC [[alternative

Hi all, I have some problems with the mle2 function > RogersIIbinom <- function(N0,attackR3_B,Th3_B) {N0-lambertW(attackR3_B*Th3_B*N0*exp(-attackR3_B*(24-Th3_B*N0)))/(attackR3_B*Th3_B)} > RogersII_B <- mle2(FR~dbinom(size=N0,prob=RogersIIbinom(N0,attackR3_B,Th3_B)/N0),start=list(attackR3_B=1.5,Th3_B=0.04),method="Nelder-Mead",data=dat) Error in dbinom(x, size, prob, log)

Hi there, I am trying to do multiple pairwise Wilcoxon signed rank tests in a manner similar to: a <- c(runif(1000, min=1,max=50), rnorm(1000, 50), rnorm(1000, 49.9, 0.5), rgeom(1000, 0.5)) b <- c(rep("group_a", 1000), rep("group_b", 1000), rep("group_c", 1000), rep("group_d", 1000)) pairwise.wilcox.test(a, b, alternative="two.sided",

Hi, I'm new to R, and I'm a bit confused with the "convolve()" function. If I do: x<-c(1, 2, 3) convolve(x, rev(x), TRUE, "open") = 9 12 10 4 1 But I expected: 3 8 14 8 3 (like in Octave/MATLAB - conv(x, reverse(x)) ) 3 2 1 x 1 2 3 = 3 2 1 0 6 4 2 0 0 9 6 3 = 3 8 14 8 3 The thing is, that "convolve(x, x, TRUE, "open")" works. For me

Full_Name: Chanseok Park Version: R 2.2.1 OS: RedHat EL4 Submission from: (NULL) (130.127.112.89) pbinom(any negative value, size, prob) should be zero. But I got the following results. I mean, if a negative value is close to zero, then pbinom() calculate pbinom(0, size, prob). dbinom() also behaves similarly. > pbinom( -2.220446e-22, 3,.1) [1] 0.729 > pbinom( -2.220446e-8, 3,.1)

I need to make weighted boxplots. I found that ggplot makes them. I would however like to label them with the boxplot statistics (the median, q1 and q3). In the boxplot function in r-base, I could output boxplot statistics and then write a text on the plot to place the labels. How would one do it with ggplot? Vikas

Am 16.07.2015 um 14:08 schrieb Martin Maechler: >>>>>> G?bor Cs?rdi <csardi.gabor at gmail.com> >>>>>> on Thu, 16 Jul 2015 07:18:00 -0400 writes: > > > ... > > making array.d from array.c > > array.c:33:23: fatal error: duplicate.h: No such file or directory > > compilation terminated. > > ... >

On 23/08/2017 6:25 PM, Bert Gunter wrote: > Doesn't sort by size of subgroup. I interpret the phrase I asterisked as: You were fooled by Peter's tricky single negative. Duncan Murdoch > > Your code does the following: > > First subsets of size 1 are given. > Then all subsets of size 2. > Then all subsets of size 3. > etc. > > Your code does not do this

Hi, The following code, from the examples in ?TkWidgets , immediately crashes R 4.0.1 for Windows: --------------------- snip -------------------- library("tcltk") tt <- tktoplevel() label.widget <- tklabel(tt, text = "Hello, World!") button.widget <- tkbutton(tt, text = "Push", command = function()cat("OW!\n")) tkpack(label.widget,

There seems to be another outage of the Vienna CRAN server [1] though from what I can from afar not related to DNS outage [2]. Could anyone confirm and possibly alert the local team in Vienna? Dirk [1] https://twitter.com/d_olivaw/status/1147161883797086214 [2] The DNS resolves on a few sites but eg simple pings stall -- http://dirk.eddelbuettel.com | @eddelbuettel | edd at debian.org

Hello, # Full disclosure. I am not sure if my problem is a bug(s) in the code, or a fundamental misunderstanding on my part about what I am trying to do with these statistics. I am not familiar with maximum likelihood tests. # I currently have two vectors Aequipecten<-c(0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0,

I am trying to use if...else loop and have included a code snippet which I might like to expand. Maybe you could steer me in the right direction. library(stats) library(prob) { a <- sample ( 1:4,100, replace=T,prob=c(0.1,0.2,0.5,0.3)) b<-sample(3:6,100,replace=T,prob=c(0.2,0.2,0.2,0.4)) } dd <- data.frame(a,b) if (subset finds a vector) ( print that vector) (else

Is there a fuction that will allow me to retrun the filename for a script from within that script. fir instance If I have a script "myscript.r": FileName<-unknown.fucntion() print(FileName) and run it source("myscript.r") will return "myscript.r" Thanks Nevil Amos [[alternative HTML version deleted]]

I'm solving the following problem: Create a variable (column) in the ?wf? dataframe named ?Zone? that takes value of ?tropic? if Latitude is less than or equal to 30, or ?non-tropic? for Latitude greater than 30. Show you Zone variable. Latitude is a column of my dataframe. I don't know the sintaxis and the process to create the column. [[alternative HTML version deleted]]

Hi all, I should have given a better explanation of my problem. Here it is. I extracted from my code the bit that gives the error. Place this in a file called test.c #include <math.h> #include <R.h> #include <Rmath.h> #include <float.h> #include <R_ext/Print.h> int main(){ double prob[3] = {0.0, 0.0, 0.0}; double prob_tot = 0.; prob[0] = 0.3*dnorm(2, 0,

Hi Guy's I was wondering if someone could point me in the right direction. dbinom(10,1,0.25) I am using dbinom(10,1,0.25) to calculate the probabilty of 10 judges choosing a certain brand x times. I was wondering how I would go about simulating 1000 trials of each x value ? regards Brendan -- View this message in context:

Installing packages from binaries on Windows seems broken, when using mirrors that are up to date with CRAN install.packages( 'httr', type = 'binary', repos = "https://cran.rstudio.com/" ) Changing repos to the Kansas CRAN mirror installs the package as expected, but that could be because the KS mirror has not yet synced. Someone pointed out that the PACKAGES.gz

switch(i, ...) extracts 'i'-th argument in '...'. It is like eval(as.name(paste0("..", i))) . Just mentioning other things: - For 'n', n <- nargs() can be used. - sys.call() can be used in place of match.call() . --------------------------- >>>>> peter dalgaard <pdalgd at gmail.com> >>>>> on Mon, 15 May 2017 16:28:42

*Hi all, I got problems installing RSTAR, MSVAR, and MSVECM packages. * > install.packages("RSTAR")Installing package(s) into ‘C:/Program Files/R/R-2.15.1/library’ (as ‘lib’ is unspecified)Warning in install.packages : package ‘RSTAR’ is not available (for R version 2.15.1) > install.packages("MSVAR") Installing package(s) into ‘C:/Program Files/R/R-2.15.1/library’

I think Herv?'s idea was just that if switch can evaluate arguments selectively, so can stopifnot(). But switch() is .Primitive, so does it from C. I think it is almost a no-brainer to implement a sequential stopifnot if dropping to C code is allowed. In R it gets trickier, but how about this: Stopifnot <- function(...) { n <- length(match.call()) - 1 for (i in 1:n) { nm