search for: noninteg

Hello all, I have a fairly simple data manipulation question. Say I have a dataframe like this: dat <- as.data.frame(runif(7, 3, 5)) dat$cat <- factor(c("1","4","13","1","4","13","13A")) dat runif(7, 3, 5) cat 1 3.880020 1 2 4.062800 4 3 4.828950 13 4 4.761850 1 5 4.716962 4 6

...e <- 29 resul <- matrix(rep(0,ncota*nslope*4),ncota*nslope,4) But this doesn't? ncota <- 1 sini <- 0.1; sfin <- 1.5; spaso <- 0.05; nslope <- 1+((sfin-sini)/spaso) resul <- matrix(rep(0,ncota*nslope*4),ncota*nslope,4) I guess the problem is that the division gives a noninteger number. How can I get the second one work? I need to create a zero matrix with its size calculated from a calculation. cheers -- View this message in context: http://r.789695.n4.nabble.com/R-Why-this-deosn-t-work-matrix-rounding-error-tp2968527p2968527.html Sent from the R help mailing list arc...

I see that "case weights" can be optioned in multinom(). I wanted to make sure I understand what weights= is expecting. My weights (not really mine but I'm stuck with them) are noninteger, are not scaled to sum to the sample size, and larger weights are intended to increase influence. The description of various types of weights is a perennial confusion for me; sorry. STS Steven Sullivan, Ph.D. Senior Associate The QED Group, LLC 1250 Eye St. NW, Suite 802 Washingto...

...inate from dpois() which is used to evaluate the deviance of the fitted model. In R versions prior to 1.3 it seemed to work without complaint. So .. (c) what is the recommended way to suppress the warning from dpois() that the responses are not integers ? Ideally I would like dpois() to accept noninteger responses and evaluate using the gamma function.. thanks ---- Adrian Baddeley, Mathematics & Statistics, University of Western Australia <http://maths.uwa.edu.au/~adrian/> -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read ht...

I was looking for an R-package that would reduce matrices to row-echelon form, but Google was not my friend; any leads? If not, I wonder if the problem could be expressed in terms of constraint satisfaction...

Hi, I have count data that have been meddled with enough to make them non integers. Using glm(poisson) returns a "non integer" error but glm(quasipoisson) does not. Just wondering if anyone knows if I am violating the assumptions of a quasipoisson error structure by using these non-integer response data? Thanks! I'd welcome your thoughts and/or references... Mark

...no way to tell what the result=20 should mean for >> non-integer m, n, k. PD> My initial reaction too (and surely it is not a bug that PD> functions behave inconsistently in regions where they PD> are not documented to work at all), but on the other PD> hand, noninteger m,n do appear to give a well-defined PD> distribution, and perhaps there's a way of making sense PD> of it? I wouldn't think it corresponds to noncentral PD> hypergeometric distributions. I'd tend to pretty much agree here. Incidentally (slightly related, but pr...

I have a question about stepAIC and extractAIC and why they can produce different answers. Here's a stepAIC result (slightly edited - I removed the warning about noninteger #successes): stepAIC(glm(formula = (Morbid_70_79/Present_70_79) ~ 1 + Cohort + Cohort2, family = binomial, data = ghs_70_79, subset = ghs_70_full),direction = c("backward")) Start: AIC=3151.41 (Morbid_70_79/Present_70_79) ~ 1 + Cohort + Cohort2 Df Deviance AIC <n...

...s-Jordan elimination). I modified it a bit: rref <- function(A, tol=sqrt(.Machine$double.eps),verbose=FALSE, fractions=FALSE){ ## A: coefficient matrix ## tol: tolerance for checking for 0 pivot ## verbose: if TRUE, print intermediate steps ## fractions: try to express nonintegers as rational numbers ## Written by John Fox if (fractions) { mass <- require(MASS) if (!mass) stop("fractions=TRUE needs MASS package") } if ((!is.matrix(A)) || (!is.numeric(A))) stop("argument must be a numeric matrix") n <- nrow(A) m <- ncol(...

This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. ---1936847065-1111238301-1101309010=:193006 Content-Type: TEXT/PLAIN; charset=iso-8859-1; format=flowed Content-Transfer-Encoding: QUOTED-PRINTABLE On Wed, 24 Nov 2004 Erik.Jorgensen@agrsci.dk wrote: > > dhyper() does not allow non-integer

Hallo Everybody Consider the following vector a=1:10 b=20:30 c=40:50 x=c(a,b,c) I need a function that can tell me that there are three set of continuos sequences and that the first is from 1:10, the second from 20:30 and the third from 40:50. In other words: a,b, and c. regards Christiaan [[alternative HTML version deleted]]

...w I should have started earlier. Kudos to you all. 5. Known issues: i. An additional test suit (book5/book6) corresponding to the not-yet-published extension of my book's appendix of validity tests was added. It found an issue with "stderr of expected survival/Cox model with noninteger case weights/Efron approximation". Not yet addressed, of a 1/n vs 1/(n-1) size. ii. The stderr of terms/survreg/penalized model is different between Splus and R by 50% or more. I haven't yet figured out which formula is correct. iii. The CMD check script still complains about so...

...data so the sum of the absence weights is equal to the sum of presence weights so that the model isn?t swamped by an overwhelming and arbitrary number of background points. I?m trying to do this in R in the standard glm in the stats package and am a bit confused. I understand the issue with noninteger weights in glm and specifying reasonable starting values to ensure convergence to something meaningful and for any individual model fit I can get mostly the same results. For example I have a dataset with 75 presence and 75 absence if instead I had randomly sampled 3150 background points I c...

I don't normally have to go anywhere near this stuff , but it seems to me that this should be a straight-forward process in R. For the purposes of this enquiry I thought I would use something I can work out on my own. So I have my matrix and the right hand results from that matrix tdata <- matrix(c(0,1,0,-1,-1,2,0,0,-5,-6,0,0,3,-5,-6,1,-1,-1,0,0),byrow = T,ncol = 5) sumtd <-