similar to: Can anyone help me reproduce this SAS Mixed output??

To those who might be interested -- following is the solution to my previous post regarding reproducing output from SAS Proc Mixed for a two-factor crossed random effects ANOVA model. I am graciously endebted to the kind replys from two statisticians for this solution whose names I will refrain from mentioning for the sake of privacy. I hope this helps someone?! -- Phil Turk > hw7 <-

> On Jun 8, 2018, at 11:52 AM, Hadley Wickham <h.wickham at gmail.com> wrote: > > On Fri, Jun 8, 2018 at 11:38 AM, Berry, Charles <ccberry at ucsd.edu> wrote: >> >> >>> On Jun 8, 2018, at 10:37 AM, Herv? Pag?s <hpages at fredhutch.org> wrote: >>> >>> Also the TRUEs cause problems if some dimensions are 0: >>>

Hmmm, yes, there must be some special case in the C code to avoid recycling a length-1 logical vector: dims <- c(4, 4, 4, 1e5) arr <- array(rnorm(prod(dims)), dims) dim(arr) #> [1] 4 4 4 100000 i <- c(1, 3) bench::mark( arr[i, TRUE, TRUE, TRUE], arr[i, , , ] )[c("expression", "min", "mean", "max")] #> # A tibble: 2 x 4

Estimada comunidad tengo un problema del que no encuentro datos que me ayuden mucho en la web. Estoy haciendo calculos por grupos con data,table. Tengo un archivo (zp.res) con tres columnas que clasifican los datos (sol, con, dia) y una columna de datos numericos (media), de la siguiente forma: sol con dia media 1: con 0 1 -22.6 2: con 0 1 -36.6 3: con 0 1 -35.6 y

Hi, R-users, data is a matrix like this AMR BS GE HR MO UK SP500 1974 -0.3505 -0.1154 -0.4246 -0.2107 -0.0758 0.2331 -0.2647 1975 0.7083 0.2472 0.3719 0.2227 0.0213 0.3569 0.3720 1976 0.7329 0.3665 0.2550 0.5815 0.1276 0.0781 0.2384 1977 -0.2034 -0.4271 -0.0490 -0.0938 0.0712 -0.2721 -0.0718 1978 0.1663 -0.0452 -0.0573 0.2751 0.1372 -0.1346

> On Jun 8, 2018, at 10:37 AM, Herv? Pag?s <hpages at fredhutch.org> wrote: > > Also the TRUEs cause problems if some dimensions are 0: > > > matrix(raw(0), nrow=5, ncol=0)[1:3 , TRUE] > Error in matrix(raw(0), nrow = 5, ncol = 0)[1:3, TRUE] : > (subscript) logical subscript too long OK. But this is easy enough to handle. > > H. > > On

Hi How would I analyse time series with - different lengths (i.e. one has 9 entries and the other has 14 entries) - different frequency (i.e. dates are random - no repeated length) - multiple values for the same time entry (e.g. 2009-10-23 below) i.e. my data takes the form: 1st time series 2009-10-07 0.009378 2009-10-19 0.014790 2009-10-23 -0.005946 2009-10-23 0.009096 2009-11-08 0.004189

Hi, How can I extract the label only from a given data frame. Fore example from this data frame. > print(dataf) V1 V2 V3 V4 V5 V6 V7 V8 V9 11145 14.3 17.1 31.2 41.7 45.8 49.8 68.6 70.6 72.9 3545 10.2 15.6 20.9 23.2 31.4 31.7 36.2 48.4 51.9 8951 15.2 17.5 20.0 21.4 32.4

Hi, I''ve adapted my ferret_killer1.rb script to reproduce the segfaults I''m still having in my app with Ferret 0.11.3. ferret_killer1 didn''t reproduce these - I had to add code to retrieve the documents from the index after searching (and fixed the index creation bug Dave spotted). http://johnleach.co.uk/downloads/ruby/ferret/ferret_killer3.rb To use, in one terminal

Hi, I run into a consistent crash in the GOG.com version of ToTSC (latest patched version). Specifically when entering house 6 as shown on http://www.forgottenwars.com/bg1/maps/x_0700.jpg Config: Ubuntu 10.10 x64, nVidia driver 270.41.03 kernel 2.6.38, wine 1.3.18 Can someone please check if they're also affected (or not)? I can provide a save game if necessary. Thanks in advance, FD

Hi Juergen, I've had a trawl through the SuSE update but it's not clear which patch might impact this issue. However SuSE are usually very good about upstreaming their fixes so I think there is a very good chance this will be fixed in the 4.1 packages in Wheezy. Are you able to give it a try? Thanks, Ian.

Dear All, when I try to reproduce an example of survexp, taken from the help page of survdiff, I receive the error message "Error in floor(temp) : Non-numeric argument to mathematical function" . It seems to come from match.ratetable. I think, it has to do with character variables in a ratetable. I would be interested to know, if it works for others. With an older version of

Dear R Users, I'm having trouble reproducing the results in Section 5.1 of Culp, M., Johnson, K., Michailidis, G. (2006). ada: an R Package for Stochastic Boosting Journal of Statistical Software, 16 They build and display a boosting model with the code: library("ada") n <- 12000 p <- 10 set.seed(100) x <- matrix(rnorm(n*p), ncol=p) y <-

The R code to reproduce the solutions to the example problems in Design of Experiments with MINITAB by Paul Mathews (ASQ Quality Press, 2005) is now posted at http://www.mmbstatistical.com/DOEwithMINITAB.html. Enjoy! Paul Mathews -- Paul Mathews Mathews Malnar and Bailey, Inc. 217 Third Street, Fairport Harbor, OH 44077 Phone: 440-350-0911 Fax: 440-350-7210 E-mail: Paul: pmathews at

On 1 Mar 2014 at 5:42, Ady wrote: From: Ady <ady-sf at hotmail.com> To: syslinux at zytor.com Date sent: Sat, 1 Mar 2014 05:42:25 +0200 Priority: normal Subject: Re: [syslinux] VirtualBox 4.1.x can reproduce bug 45 Send reply to: For discussion of Syslinux and tftp-hpa <syslinux at zytor.com> > > > > > Short summary:

Hello all, did anybody try to use the application under Wine? I need help in reproducing and researching the bug. The program requires modem and I don't have one. For details see: http://bugs.codeweavers.com/show_bug.cgi?id=411 The page has description how to reproduce and gather information about the issue. Please, add your comments to the bug. In general - you are welcome to submit

>>>>> Winston Chang <winstonchang1 at gmail.com> >>>>> on Tue, 4 Apr 2017 15:29:40 -0500 writes: > I've done some more investigation into the problem, and it is very > difficult to pin down. What it looks like is happening is roughly like this: > - `p` is an environment and `p$e` is also an environment. > - There is a loop.

On 5 April 2017 at 15:46, Winston Chang wrote: | On Wed, Apr 5, 2017 at 2:24 PM, Robert McGehee <rmcgehee at walleyetrading.net> | wrote: | | > Winston, | > I had a similar experience to you tracking down an insanely difficult bug | > in my R code that "disappeared" whenever slight changes were made to the | > script (e.g. like adding cat() statements). In my case, it

> On 05 Apr 2017, at 20:40 , Winston Chang <winstonchang1 at gmail.com> wrote: > > I think there's a good chance that this is due to a bug in R. I have > been trying to track down the cause of the problem but haven't been > able find it. > > -Winston Apologies in advance if this is just stating the obvious, but let me try and put some general ideas on the

On 05.04.2017 23:54, peter dalgaard wrote: > >> On 05 Apr 2017, at 20:40 , Winston Chang <winstonchang1 at gmail.com> wrote: >> >> I think there's a good chance that this is due to a bug in R. I have >> been trying to track down the cause of the problem but haven't been >> able find it. >> >> -Winston > > Apologies in advance if this