search for: henric

>>>>> Henric Winell <nilsson.henric at gmail.com> >>>>> on Tue, 7 Feb 2017 13:37:42 +0100 writes: > Hi, On 2017-02-07 13:12, Benjamin Tyner wrote: >> Hello >> >> When loading a package, I'm wondering if it's frowned >> upon for t...

...ecking. Since my data are in binary form I'd like to collapse over the cross-classification of the factors before the model checking. Are there any nice and simple ways doing this? If so, how? If not, I'd be grateful for receiving some hints on how this can be accomplished. Many thanks, Henric --------------------------------------------------------------------------------------- Henric Nilsson, Statistician Statisticon AB, ?stra ?gatan 31, SE-753 22 UPPSALA Phone (Direct): +46 (0)18 18 22 37 Mobile: +46 (0)70 211 68 36 Fax: +46 (0)18 18 22 33 <http://www.statisticon.se>

...predict.prcomp* [13] predict.princomp* predict.smooth.spline* [15] predict.smooth.spline.fit* predict.StructTS* Non-visible functions are asterisked I'm sure there's a good reason for this, but I haven't been able to figure it out. Please enlighten me. Thanks! Henric

...0 > seq(0.7, 0, by = -0.1) [1] 7.000000e-01 6.000000e-01 5.000000e-01 4.000000e-01 3.000000e-01 2.000000e-01 1.000000e-01 -1.110223e-16 Is this really the intended behaviour? I ended up using > seq(0.7, 0, length = 8) [1] 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 which does what I want. //Henric

The last line of the example in droplevels' manual page seems to be incorrect to me. I think it should read: "table(droplevels(aq$Month))". Amazingly (I don't understand) both variants seem to produce the same result (R 3.3.3): --- > aq <- transform(airquality, Month = factor(Month, labels = month.abb[5:9])) > aq <- subset(aq, Month != "Jul") >

Hello, Inline. Em 12-04-2017 16:40, Henric Winell escreveu: > (Let's keep the discussion on-list -- I've added back R-devel.) > > On 2017-04-12 16:39, Ulrich Windl wrote: > >>>>> Henric Winell <nilsson.henric at gmail.com> schrieb am 12.04.2017 >>>>> um 15:35 in >> Nachricht <...

...> > x <- c(1, 1, 2, 3) > > rank2(x, ties.method = "last") > [1] 1 2 4 3 > > That doesn't look right to me -- I had expected > > > rev(sort.list(x, decreasing = TRUE)) > [1] 2 1 3 4 > Indeed, well spotted, that seems to be correct. > > Henric Winell > ------------------------------ In the particular example (of length 4), what is really wanted is the following. ind <- integer(4) ind[sort.list(x, decreasing=TRUE)] <- 4:1 ind The following gives the desired result: sort.list(rev(sort.list(x, decreasing=TRUE)))

...this function always estimates the dispersion using Deviance/df. On the other hand, the cooks.distance function uses the Pearson residuals and extracts the dispersion from the fitted model using summary, which to me seems more sensible for a GLM. Can someone please shed some light on this? //Henric

...kage ,whick library should I include in R ? Best Regards! Ivy Li YMS in Production & Testing Semiconductor Manufactory International(ShangHai) Corporation #18 ZhangJiang Road, PuDong New Area, Shanghai, China Tel: 021-5080-2000 *11754 Email: Ivy_Li at smics.com -----ԭʼÓʼþ----- ·¢¼þÈË: Henric Nilsson [mailto:henric.nilsson at statisticon.se] ·¢ËÍʱ¼ä: 2004Äê10ÔÂ15ÈÕ 17:39 ÊÕ¼þÈË: Ivy_Li ³­ËÍ: r-help at stat.math.ethz.ch Ö÷Ìâ: Re: [R] R plot problems At 16:17 2004-10-15 +0800, you wrote: >[...] I want to rotate the direction of x-coordinates' letter so that it >can show al...

>>>>> Rui Barradas <ruipbarradas at sapo.pt> >>>>> on Wed, 12 Apr 2017 17:07:45 +0100 writes: > Hello, Inline. > Em 12-04-2017 16:40, Henric Winell escreveu: >> (Let's keep the discussion on-list -- I've added back >> R-devel.) >> >> On 2017-04-12 16:39, Ulrich Windl wrote: >> >>>>> Henric Winell <nilsson.henric at gmail.com> schrieb am >> 12.04.2...

...or and x is a continuous covariate. If I want to fit a similar gam model, is it correct to fit y~a+s(x)+s(x,by=a.1)+s(x,by=a.2)+s(x,by=a.3), where a.1--a.3 are dummy variables representing each level of the factor? Or is the s(x) term redundant? Any hints are greatly appreciated. Best wishes, Henric --------------------------------------------------------------------------------------- Henric Nilsson, Statistician Statisticon AB, ?stra ?gatan 31, SE-753 22 UPPSALA Phone (Direct): +46 (0)18 18 22 37 Mobile: +46 (0)70 211 68 36 Fax: +46 (0)18 18 22 33 <http://www.statisticon.se>

(Let's keep the discussion on-list -- I've added back R-devel.) On 2017-04-12 16:39, Ulrich Windl wrote: >>>> Henric Winell <nilsson.henric at gmail.com> schrieb am 12.04.2017 >>>> um 15:35 in > Nachricht <b66fe849-bb8d-f00d-87e5-553f866d57e0 at gmail.com>: >> On 2017-04-12 14:40, Ulrich Windl wrote: >> >>> The last line of the example in droplevels' manual pag...

Hello When loading a package, I'm wondering if it's frowned upon for the package to alter the state of the random number generator? I guess not, since the parallel package does it? > set.seed(6860) > old.seed <- .GlobalEnv$.Random.seed > library(parallel) > new.seed <- .GlobalEnv$.Random.seed > identical(old.seed, new.seed) [1] FALSE I ask

...believing that this is the intended behaviour, but I've been proved wrong before... Try library(lattice) y <- x <- 1:10 par(mfrow = c(2, 2)) # plots on a new page plot(y ~ x) # plots on a new page xyplot(y ~ x) # overplotting occurs plot(y ~ x) for a reproducible example. Best wishes, Henric

In Henric's recent post, he included this output: @BOOK{R:Harrell:2001, AUTHOR = {Frank E. Harrell}, TITLE = {Regression Modeling Strategies, with Applications to Linear Models, Survival Analysis and Logistic Regression}, PUBLISHER = {Springer}, YEAR = 2001,...

...e, 71, 175-179. 4. For each purely parametric covariate a t-test is produced; I'd like to have something like S-plus' anova.gam() to get an overall test. (Perhaps with the addition of a choice between Type I and Type III tests, but I guess that may be controversial). Is it possible? //Henric --------------------------------------------------------------------------------------- Henric Nilsson, Statistician Statisticon AB, ?stra ?gatan 31, SE-753 22 UPPSALA Phone (Direct): +46 (0)18 18 22 37 Mobile: +46 (0)70 211 68 36 Fax: +46 (0)18 18 22 33 <http://www.statisticon.se>

...ized version of my package can no longer reproduce any subsequent results depending on random number generation (unless a call to 'set.seed' was issued *after* attaching my package). I'd be most grateful for any help that you're able to provide here. Many thanks! Kind regards, Henric Winell > sessionInfo() R Under development (unstable) (2014-01-26 r64897) Platform: x86_64-redhat-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=sv_SE.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8 [7] LC_...

Hello R community members, I'm trying to install the 'rqpd' package which is developed by Roger Koenker and Stefan Bache. When I try to install the package using the command 'install.packages("rqpd",repos="http://R-Forge.R-project.org")' I'm getting the following two messages: i) package ?rqpd? is available as a source package but not as a binary

On Tue, Oct 20, 2015 at 10:26 AM, Henric Winell <nilsson.henric at gmail.com> wrote: > Den 2015-10-09 kl. 12:14, skrev Martin Maechler: > I think so: the code above doesn't seem to do the right thing. Consider > the following example: > > > x <- c(1, 1, 2, 3) > > rank2(x, ties.method = "last&qu...

...-08 variety 1 3.065 18 39.686 0.211 host:variety 1 6.408 17 33.278 0.070 As expected: > 1-pchisq(6.408/dscale,1) [1] 0.07040576 2. When using summary.glm() on a glm object fitted using family=quasibinomial the reported tests are t-tests. Why? Thanks, Henric