similar to: Problem in applying conditional looping

Respected All, I am writing a program in R and facing some problem with applying "if statment". Program first draw random numbers from bivariate normal distribution then compute variable say Pi and Pij from that sample and then further computation ..... ..... ..... In some samples Pij is appearing with negitive sign and ultimately resulted in an negative variance estimator.

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 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 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

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

I am looking for R code to be able to fit a linear-linear piecewise model with person-specific changepoint. I have searched the web, but have not been able to locate any code. Below is my attempt at some code: chgpt = function(a1,a2,a3,gam,wave){ yht=numeric(10) y1=(wave <= gam)*(a1+(a2*wave)) y2=(wave > gam)*((a1+(a2-a3)*gam)+a3*wave) yhat=y1+y2 return(yht) } nl.dat <- nlme(y ~

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 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

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.

Hi Atif make is a Unix's command which uses Makefile file for package's compilation. So after installing the complete development package from distribution disk, launch make. Ciao mauro

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

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