similar to: generating clustered data

Hello all: I am interested in computing what the multilevel modeling literature calls a multiple membership model. More specifically, I am working with a data set involving clients and providers. The clients are the lower-level units who are nested within providers (higher-level). However, this is not nesting in the usual sense, as clients can belong to multple providers, which I understand

Dear R users, I am estimating Poisson mixed models using glmmPQL (MASS) and lmer (lme4). We know that glmmPQL do not provide the correct loglikelihood for such models (it gives the loglike of a 'pseudo' or working linear mixed model). I would like to know how the loglike is calculated by lmer. A minor question is: why do glmmPQL and lmer give different degrees-of-freedom for the same

I have read some very good reviews comparing R (or Splus) to SAS. Does anyone know if there are any reviews comparing R (or Splus) to Stata? I am trying to get others to try R in my department, and I have never used Stata. Regards, -Cody Cody Hamilton, Ph.D Institute for Health Care Research and Improvement Baylor Health Care System (214) 265-3618 This e-mail, facsimile, or letter

Where is everyone located? I am in Washington DC. -------------- next part -------------- An HTML attachment was scrubbed... URL: http://lists.digium.com/pipermail/asterisk-users/attachments/20091106/7c73847d/attachment.htm

Hello R users: I am fairly new to R and am trying to figure out how to compute an intraclass correlation (ICC) and/or design effect for binary data? More specifically, I am trying to determine the amount of clustering in a data set - that is, whether certain treatment programs tend to work with more or less severe clients. The outcome variable is dichotomous (low severity / high severity)

Here is a very very basic config. But, not working (: I simply want to dial the DID that is registered with the SIP provider. then, as you can see the call should dial the 703111 number Hints please? sip.conf ;register => 908366554:396444 at carrier.jazzey.com register => 908366554:396444 at sip.jazzey.com [jazzey] type=friend host=sip.jazzey.com username=908366554 secret=396444

For an application, I have formulas defined in terms of a left Kronecker product of matrices, A,B, meaning A \otimes_L B = {A * B[i,j]} -- matrix on the left multiplies each element on the right. The standard kronecker() function is the right Kronecker product, A \otimes_R B = {A[i,j] * B} -- matrix on the right multiplies each element on the left. The example below shows the result of

Hi I'm using kronecker() with a matrix and a vector. I'm interested in the column names that kronecker() returns: > a <- matrix(1:9,3,3) > rownames(a) <- letters[1:3] > colnames(a) <- LETTERS[1:3] > b <- c(x=1,y=2) > kronecker(a,b,make.dimnames=TRUE) A: B: C: a:x 1 4 7 a:y 2 8 14 b:x 2 5 8 b:y 4 10 16 c:x 3 6 9 c:y 6 12 18 > The

I'm trying to calculate the row-wise kronecker product A \Box B of two sparse matrices A and B, and am struggling to find a quick way to do this that takes advantage of sparseness. I thought a good idea would be to use "rep" to construct 2 matrices of the same dimension of the end product, and multiply these two together: library(Matrix) A<-Matrix(c(1,0,0,0,0,1,2,0), 2, 4)

kronecker, with make.dimnames=TRUE uses a hardwired sep=":" in the line tmp <- outer(dnx[[i]], dny[[i]], FUN = "paste", sep = ":") For an application in which dimnames arise from an n-way array, where different dimensions have different roles, and I would like to be able to use kronecker in the form kronecker(A, B, make.dimnames=TRUE,

why does this not work? i simply want to hear the recorded message exten => s,1,Answer() ;exten => s,n,Record(zipcodegutter1.gsm) ;zcg1 exten => s,n,Playback(zipcodegutter1) exten => s,n,Dial(SIP/c000001s/12222222259,120,A,(demo-thanks))

Dear All I am kind of stuck up with a code a part of which seems to be causing a problem, or at least I think so. May be the community can help me. It’s simple but I suppose I am missing something. I generate a data matrix X, say of order n*p, where n represents independent row-vectors and p correlated col vectors. Let the row representation be X = (X’_1, . . ., X’_n)’. I generate the

Does anyone have a Dial script or a hint on how I can dial 10000 numbers in sequence? When the calls are answered, I play a .gsm or .wav. Then, if user presses a defined digit, the call gets bridged to me.

Dear R users, I have a function (below) which encompasses several tests. However, when I run it, only the output of the last test is displayed. How can I ensure that the function root(var) will run and display the output from all tests, and not just the last one? Thank you, b. root <- function(var) { #---Phillips-Perron PP.test(var, lshort = TRUE) PP.test(var, lshort = FALSE)

I want to create a script for IVR that compiles responses, aggregates them to a total number. Then, run an equation based on the result. Press 1 for X (X is a positive number 500) Press 2 for Y (Y is a positive number 200) Press 3 for Z (Z is a positive number 300) Press 20 to calculate the results = 500+200+300 =1000 then, exten => s,n,Read(NUMBER,,1000) exten => s,n,SayDigits(${NUMBER})

Hello all, I just started working with the lme4 package to estimate a multilevel logistic regression and am planning to use this package for a cross-classification / multiple-membership model. I haven't found many worked examples and am trying to figure out how to add variables to the higher-level part of the model. Consider the following example: test.1 <- lmer(y ~ b1 + b2 + b3 +

Dear All, Could someone please suggest how to find the Kronecker sum of two 2x2 matrices, i.e. given two matrices: -A A a -a and -B B b -b I need: -A-B A B 0 a -a-B 0 B b 0 -A-b A 0 b a -a-b Many thanks, Lara [[alternative HTML version deleted]]

So, https://bugs.llvm.org/show_bug.cgi?id=32056 is an example showing our current TBAA tree for union generation is definitely irretrievably broken. I'll be honest here. I'm pretty sure your proposal doesn't go far enough. But truthfully, I would rather see us come closer to a representation we know works, which is GCC's. Let me try to simplify what you are suggesting, and what we

Hello all, A colleague at work set me the challenge to convert some MATLAB code into R, to see which is faster. We'd seen that benchmark comparing MATLAB 6.5 to R1.90 (and others), and so I thought that I should be able to get roughly comparable speeds. The code has lots of multiplications of matrixes, transposes, and MATLAB's "repmat". I did the code conversion, and R was about

Anyone have a cookbook on configuring sendmail to work with Asterisk? Or,a few config examples.