similar to: converting an R function into VBA

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, You may try: set.seed(49) m1 = matrix(rnorm(30), nrow = 3) m2 = matrix(rnorm(30), nrow = 3) ?corsP<-vector() ? for(i in 1:3) corsP[i] =? cor(m1[i,], m2[i,]) ?corsP #[1]? 0.58411274 -0.02382329? 0.03760757 diag(cor(t(m1),t(m2))) #[1]? 0.58411274 -0.02382329? 0.03760757 #or mNew<- rbind(m1,m2) ?indx<-rep(seq(nrow(mNew)/2),2) ?sapply(split(seq_len(nrow(mNew)),indx),function(x)

Dear All, It is my first time using this mail list to post a question. And I sincerely hope that this will not bother any subscribers. So far as I know, there are functions like union( ), which can help to combine two sets of discrete data. But what if the data sets are with continuous data. For instance, I got three sets like [2, 4], [5, 6], [3.4, 5.5]. I wonder if there is a quick solution to

Dear all, I am finding difficulty in the following, I would like to create an empty matrix e.g. 10x10 of 0s and sequentially fill this matrix with randomly placed a 1s until it is saturated. Producing 100 matrices of sequentially increasing density., This process needs to be randomized 1000 times., I assume i should run this along the following lines, 1) Create 1000 matrices all zeros, 2) add

Hello R users, I'd like to ask a question about how to add a new column. So, below is my situation. In order to perform the repeated ANOVA, I first imported the following table. score=read.csv("patients_tests.csv"); subject test1 test2 test3 test4 test5 test6 test7 1 ab 0.17687 0.16715 0.17009 0.16480 0.16116 0.24502 0.17975 2 cl

On 12/15/2016 02:44 AM, Sorin Srbu wrote: >> -----Original Message----- >> From: CentOS [mailto:centos-bounces at centos.org] On Behalf Of Sorin Srbu >> Sent: den 15 december 2016 08:16 >> To: CentOS mailing list <centos at centos.org> >> Subject: Re: [CentOS] OT: Blank mails from this list >> >>> -----Original Message----- >>> From:

Easy one. I have win98 clients that login and connect to the shares fine. My problem is the user home directory, each user can see their folders and can even open them but they cannot save to their home directory. Each users home dir is setup as 0700 (-rwx------), as are the files and subdir's, which I believe is Read, Write, Execute for the owner of the file. What am I missing? Chris.

I checked the R procedure HCLUST (hierarchical clustering) but it looks like it requires a full triangular n x n similarity matrix as input, where n = number of observations. The number of variables is 200. My data set has n = 50,000 observations (keywords), and I use ad-hoc similarity measures, not available in R, to measure keyword similarity. Here, the vast majority of the n x n similarities

Dear all, I have discovered that there are many things that I used to do in my GIS which are easily done directly in R, for example calculating interpoint distances using geoR and pick out points inside a polygon using splancs. I now wonder, is there a function to create a line object like a watercourse and then calculate the distances between many points in space and this line? I couldn't

Hello, I have a situation where a host, which is booted off its ''rpool'', need to temporarily import the ''rpool'' of another host, edit some files in it, and export the pool back retaining its original name ''rpool''. Can this be done ? Here is what I am trying to do: # zpool import -R /a rpool temp-rpool # zfs set mountpoint=/mnt

Greetings: I have what might be an oddball question. Is there any way to configure the times at which clients initiate their shutdown sequences? I have a domain controller in my setup that I would have shut down if the UPS have been on backup for any more than a minute, whereas the rest of the servers I would shut down on a low battery condition. I also have a use for such flexibility in our

I've been working on writing a subroutine to page groups of phones at once and I'm having some difficulty. My goal is to have a user call an extension, I record the page they wish to play, I then page out that recorded file to the phones in groups. [sub-masspage] exten => s,1,NoOP same => n,Answer same => n,Set(filename=$PAGE) same => n,Wait(1) same =>

I'm using R 1.7.1 under linux redhat it seems that the eigen values produced by eigen() do not follow a consistant order; I mean either ascending or discending e.g for one system: eigenV<-eigen(V) > print(eigenV$values) [1] -7.706828e+13 -4.702980e+13 -3.267579e+13 -1.701297e+13 -8.041677e+12 [6] -5.707311e+12 -5.053941e+12 -4.774652e+12 -4.280423e+12 -3.798905e+12

Hello: I modified benchmark used in http://www.informatik.uni-frankfurt.de/~stst/benchspl.txt now the attached code will work with R, except for last test for some reason (Sorry ran out of time to play with this). Immediate problem is that R does not show fraction of seconds. So, that kinda makes comparison some what hard. I also looked at memory and CPU usage on NT console. R seems to go

hi! i read some other threads but none could help me. i''m on gentoo trying to boot 2.6.24 and 2.6.25 kernel.org kernels within my 2.6.20-xen kernel but i fail to get a correct .config. thanks yannick more infos: -- dom0 alchemist ~ # xm info host : alchemist release : 2.6.20-xen-r6-yk version : #1 SMP Sat Dec 1 13:47:54 CET 2007

useR's, Consider the variables defined below: yvals <- c(25,30,35) x1 <- c(1,2,3) x2 <- c(3,4,5) x3 <- c(6,7,8) x <- as.data.frame(cbind(x1,x2,x3)) delta <- c(2.5, 1.5, 0.5) h <- delta/2 vars <- 3 xk1 <- seq(min(x1)-0.5, max(x1)+0.5, 0.5) xk2 <- seq(min(x2)-0.5, max(x2)+0.5, 0.5) xk3 <- seq(min(x3)-0.5, max(x3)+0.5, 0.5) xks <- list(xk1,xk2,xk3) xk <-

Dear WizaRds, I am struggling with the use of twophase in package survey. My goal is to compute a simple example in two phase sampling: phase 1: I sample n1=1000 circuit boards and find 80 non functional phase 2: Given the n1=1000 sample I sample n2=100 and find 15 non functional. Let's say, phase 2 shows this result together with phase 1: ...................phase1........

hi, is there any wait function in R. I am running one R script to plot many graphs it is in the for loop. its showing no error but its not plotting well I think i can solve this problem with a wait function. Please help me in this regards. If u need any clarification about programme. u can find the script below. best regards, Deepak.M.R Biocomputing Group University of Bologana. #!/usr/bin/R

Hi everyone! I'm quite a newbie at this Asterisk stuff so please bear with me. We've recently decided to start training in Asterisk via AsteriskNow! Asterisk version is 1.4.18.1 through AsteriskNow! 1.02 The box we have is paired with a Digium TE110P and we've managed to get it to the point where incoming calls via a single DID (from NTT Japan) can be received and answered

There is a series of data contains time in fixed step and energy varying with time, how to test its periodicity?In R, it seems there is no direct tools since I have search the R manual with periodic and I have not found any related topic. Thanks a lot