similar to: reformatting some data

Hi, Try this: #mat1 is the data res<-do.call(cbind,lapply(seq_len(nrow(mat1)),function(i) {new1<-do.call(rbind,lapply(seq_len(nrow(mat1[-i,])),function(j) {x1<-rbind(mat1[i,],mat1[j,]); x2<-(abs(x1[1,1]-x1[2,1])*abs(x1[1,5]-x1[2,5]))+(abs(x1[1,2]-x1[2,2])*abs(x1[1,6]-x1[2,6]))+(abs(x1[1,3]-x1[2,3])*abs(x1[1,7]-x1[2,7]))+(abs(x1[1,4]-x1[2,4])*abs(x1[1,8]-x1[2,8]))}));new1}))

Dear R user I want to create a table (as below) to summarize the attached data (Test.csv, which can be read into R by using 'read.csv(Test.csv, header=F)' ), to indicate the day that there are any data available, e.g.value=1 if there are any data available for that day, otherwise value=0. 28/04 29/04 30/04 01/05 02/05 532703 0 1 1

Dear friends, I have a weekly time series which starts on Jan 4th, 2003 and ends on december 31st, 2016. I set up my ts object as follows: MyTseries <- ts(mydataset, start=2003, end=2016, frequency=52) MyModel <- auto.arima(MyTseries, d=1, D=1) MyModelForecast <- forecast (MyModel, h=12) Since my last observation was on december 31st, 2016 I expected my forecast date to start on

Een ingesloten tekst met niet-gespecificeerde tekenset is van het bericht gescrubt ... Naam: niet beschikbaar Url: https://stat.ethz.ch/pipermail/r-help/attachments/20061106/756382a8/attachment.pl

Hello R-sig-debian and (hopefully) Dirk: On Debian wheezy, I have the R packaging that CRAN (you) provide. I run into a little trouble while trying to fiddle with alternative BLAS. I know you and I went around on this last year and I think perhaps I've found something wrong in the framework, or I've just done something wrong. I installed the packages openblas-base and openblas-dev, and

> On Dec 6, 2017, at 5:07 AM, Paul Bernal <paulbernal07 at gmail.com> wrote: > > Dear friends, > > I have a weekly time series which starts on Jan 4th, 2003 and ends on > december 31st, 2016. > > I set up my ts object as follows: > > MyTseries <- ts(mydataset, start=2003, end=2016, frequency=52) > > MyModel <- auto.arima(MyTseries, d=1, D=1)

Thank you very much David. As a matter of fact, I solved it by doing the following: MyTimeSeriesObj <- ts(MyData, freq=365.25/7, start=decimal_date(mdy("01-04-2003"))) After doing that adjustment, my forecasts dates started from 2017 on. Cheers, Paul 2017-12-06 12:03 GMT-05:00 David Winsemius <dwinsemius at comcast.net>: > > > On Dec 6, 2017, at 5:07 AM, Paul

I can't get my simultaneous equations to work using system fit. Please help. #Reproducible script Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv") View(Empdata) str(Empdata) Empdata$gnipc<-as.numeric(Empdata$gnipc) install.packages("systemfit") library("systemfit") pdata <- plm.data(Empdata,

Am 22.07.2015 um 17:41 schrieb Windsor Dave (AdP/TEF7) <Dave.Windsor at us.bosch.com>: > Sorry for the top post, Outlook defaults strike again..... Outlook forces you to write above ? :-) -- LF

Dear all, As I am new to the R community - although eager to advance- I would like to pose a question to the community. I have an SPSS file which I have imported it in R (with the read.spss command) which conists of scale (continuous) variable "adiponectin" and the corresponding categorical value "death" (0=No, 1=Yes). In all there are 60 observations (among which

I was looking at the manpage for rsyslog.conf, primarily because I need to filter my manager's new fedora 22 logs coming to our loghost, because of the bug that I forwarded (if it gets through). At any rate, I am surprised: under selectors, I see that " The keywords error, warn and panic are deprecated and should not be used anymore." Huh? If I only want warn or more severe,

Hi I have attempted to install Samba on an ageing Motorola running SVR3. Unfortunately it has a limitation of 14 characters for filenames, which means that when I untar the package, many of the files are not created. Does anyone have a solution to this problem, or am I stuck? Many thanks -- Pete ---------------------------------------------------------------- morph@morhp.dircon.co.uk

Just thought I'd pass along an error I received this weekend after upgrading two servers to 5.6. Both are HP Proliant DL380 servers running 64-bit, but one is a G6 model and one is a G7. After the upgrade and the reboot to the 2.6.18-238.9.1.el5 kernel, both servers displayed the following error at boot: "pci_mmcfg_init marking 256MB space uncacheable." Some Googling

Hi Firstly apologies but i have to install an app that uses postgres but i am more a MySQL/Oracle guy so i maybe doing something really wrong here.. I am following instructions for creating a DB user and encountering this # createuser -ADP -U postgres foo Enter password for new user: Enter it again: createuser: could not connect to database template1: FATAL: IDENT authentication failed for

Hello, I have current data from a nortek ADP, which is basically current speed and direction data in a 3 dimensional X Y Z format http://www.nortekusa.com/usa/products/current-profilers/aquadopp-profiler-1 The instrument logs data in a complex way and I was wondering if anyone has had any experience using R to at least read the data, and perhaps smooth it as well. If so, are there any

Hi I have been facing a strange connectivity problem from CentOS and Fedora boxes. When I telnet to $ telnet adp.eease.com 443 it works fine some times and every 5-10 tries or so it says connection refused. And this does not happen from my windows box. Can anyone from this list try this on their linux installation and see if they get the same behaviour? Thanks Tharun

Dear R users, I created a matrix that tells me the first day of use of a category by id. #Calculate time difference test$tdiff<-as.numeric(difftime(as.Date("2002-09-01"), test$ftime, units = "days")) # obtain the index date per person and dcategory index.date.test<-tapply(test$tdiff, list(test$id, test$rcat), max) Nonetheless, at the moment I think will be

______________________________________________________ Get Your Private, Free Email at http://www.hotmail.com

OK, Let's try this again! Here is the reproducible script; it is long because I had to copy the panel dataset here. My question is related to systemfit; I don't know how to get the result for the entire panel. #Reproducible script Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv") View(Empdata) install.packages("systemfit")

Hello Everyone, I recently installed CentOS 5.5 x86_64 on a brand new ProLiant DL380 G7. I have identical OS software running reock-solid on two other DL380 ProLiant servers, but they are G6 models, not G7. On the G7, the installation went perfectly and the machine ran great for about 2 weeks, when it just seemed to "stop". The system stopped responding on the network, and there was