similar to: Problem in conditional looping

Respected All, I hope you are enjoying good health, I am tring to write a program in R but could not be very sucessful. My program draws random sample form bivariate normal distribution and then compute a variable PIJ. For certian samples some entries of variable PIJ is apearing as negative, which result in negative variance estimator. I want to introduce a loop in my program that verify the each

Respected All, I need some efficient program or package to draw all possible samples of size two without replacement. I am using "combinat" package to list all possible samples but it hangs my computer for larger populations say 10,000 (i.e. 49995000 all possible samples). I wish to even work for larger populations then this and replicate this procedure for many times. Kindly can anyone

hello R: we have a two-parameter IRT simulation code. The goal is to generate a response matrix.But the "for" part doesn't run. we don't know what is wrong with it. Thanks so much~~~ I <- 10 J <- 5 response <- matrix(0, 10, 5) pij <- function(a,b,theta) { a <- rnorm(J, 0.8, 0.04) a b <- rnorm(J, 0, 1) b theta <- rnorm(I, 0,1) theta for( i in 1:I ) { for(

Hi, I have network data in the form of a couple of edgelists containing weights in the format "x,y,weight" whereby x represents row header and y represents column header. All edgelists are based on links among 634 nodes and I need to convert them into a 634*634 weighted matrix. I searched for online help using possible keywords I know, but I could not find a clue how to do this in

Please suggest a way out to the following problem. I have a T by n data matrix (say Y) where coulmns are time series of length T. To do some analysis in WinBUGS I need to construct my data as follows. yy<-rep(Y,k) ## this will be a vector Yk<-array(yy,dim=c(T,n,k)) ## data array Here the definition of dim indices is first index: T rows second index: n columns third index: for kth T by n

Hi Friends, I'm trying to model the consumer decisions (Click-Through Rate and Conversion) in Search Engine Advertising using a hierarchical Bayesian binary logit. The input data is the weekly CTRs and Avg. Position for each search keyword. CTR is modeled as (for each keyword i and week j): Pij = exp(C + Bi x Positionij + A1 x Lengthi + A2 x Brandi + A3 x ProductSpecifici) / [1 + exp(C +

I am trying to solve one of my graduate level physics problems with ruby... Here is what I have so far... a6=0.0 for n1 in -10..10 for n2 in -10..10 for n3 in -10..10 if n1!=0 and n2!=0 and n3!=0 p=Math.sqrt(n1**2+n2**2+n3**2+n1*n2/1.414+n1*n3/1.41+n2*n3/1.414) a6+=(1/p)**6 end end end end puts a6 What I''ve got here is a 10x10x10 face-centered

Hi, I am facing problem in sync time within DomU with Dom0. After gone through mailing lists, I have set /proc/sys/xen/independent_wallclock to 0. But still DomU are not syncing their time :-S Machine using kernel-xen-2.6.19-1.2911.fc6 xen-3.0.3-3.fc6 xen-libs-3.0.3-3.fc6 Please advice to solve this weird problem, Thanks _______________________________________________ Xen-users mailing list

Dear UseRs! We are pleased to announce the release of our new package 'PVAClone'. The 'PVAClone' package implements Population Viability Analysis (PVA) methodology using data cloning. The data cloning algorithm by Lele et al. (2007, 2010) is employed to compute maximum likelihood estimates of the state-space model parameters and the corresponding standard errors, heavily

Dear UseRs! We are pleased to announce the release of our new package 'PVAClone'. The 'PVAClone' package implements Population Viability Analysis (PVA) methodology using data cloning. The data cloning algorithm by Lele et al. (2007, 2010) is employed to compute maximum likelihood estimates of the state-space model parameters and the corresponding standard errors, heavily

Dear samba orgnisation, i have installed my samba in my linux server, in my network there are 5 machines, 1 is linux and other are win2k, i have configured the file /etc/samba/smb.conf, i change workgroup and etc... and i make a "testparm", it's all right, and i restart my samba, it's right. i can see my linux name from windows network but when i click it the message "The

Hi there, I'm new to R but I thought this is the most likely place I could get advice or hints w.r.t the following problem: I have a series of measurements xi with associated uncertainties dxi. I would like to construct the probability density histogram of this data where each density estimate has an associated error that is derived from the dxi. In other words, for large dxi the

Dear all, Is there a function in R that checks transitivity and acyclicity of a given nXn matrix with entries representing a decision-maker's comparisons of n objects? Like 0 1 0 1 1 1 0 0 0 1 0 0 1 0 0 0 1 1 0 0 1 0 0 0 0 1 0 1 0 1 0 1 0 1 0 0 1 represents xPy and 0 represents ~xPy. Is there a vectorized solution to this? n can be quite large. Thanks in advance, Alex

Hi, I am trying to install asterisk on Suse 9. I downloaded asterisk-1.2.1.tarand untar this package. I am following the README and the installation instruction to run "make" ans "make install". But I can not find any "make" or "make install" in the directory asterisk-1.2.1. Can any one please help me how can I install asterisk-1.2.1 on Suse? What am I

Dear R Community, I am happy to introduce the next release of the *dclone* package (see R Journal, 2(2): 29-37, 2010) which now includes a suite of functions for parallel computations: - see functions parJagsModel, parUpdate, parCodaSamples using 'snow' type clusters; and the - jags.parfit function that can be used with 'snow' type clusters or multiple cores. - iterative model

Dear R Community, I am happy to introduce the next release of the *dclone* package (see R Journal, 2(2): 29-37, 2010) which now includes a suite of functions for parallel computations: - see functions parJagsModel, parUpdate, parCodaSamples using 'snow' type clusters; and the - jags.parfit function that can be used with 'snow' type clusters or multiple cores. - iterative model

Hi, when I am connecting to samba, openingor saving a MS Word document, nearly every time there is a core-file. In the log-file there I found following message: [2002/08/12 15:13:44, 0] lib/util_sec.c:assert_gid(111) Failed to set gid privileges to (-1,-24) now set to (0,0) uid=(0,0) [2002/08/12 15:13:44, 0] lib/util.c:smb_panic(1064) PANIC: failed to set gid What does it mean? The

Yes unfortunately I do not have much luck recieving help from the SAMBA mailing list. After much playing I found something weird that was seemingly causing incorrect comments. Most printers would recieve the comment - "no entries" which originally I assumed to be the default comment given by Win95 when it could not get a proper comment. However I later found that this is not the case.

Also, a lot of measurement/research on deployment of OpenSSH rely on version advertising for their statistics. It's going to be harder to know impact of deprecation of certain legacy features without statistics. I also agree with Mark here. On Wed, Feb 20, 2019 at 10:57 AM Mark D. Baushke <mdb at juniper.net> wrote: > Nagesh writes: > > > Cyber security team has

I think most of us would expect prod(0:1000) to return 0, and ... ... it does. However, many of us also expect prod(x1, x2) to be equivalent to prod(c(x1,x2)) the same as we can expect that for min(), max(), sum() and such members of the "Summary" group. Consequently, prod(0, 1:1000) should also return 0, but as you see, it gives NaN which may be a bit puzzling... The