similar to: Changing selected elements of an array

Hi Guys, I have this part of a program: library(survival) Gastric <- cbind.data.frame(Gp=c(rep(1,45),rep(0,45)), ### 2nd gp 0 time=c(1,63,105,129,182,216,250,262,301,301,342,354,356,358, 380,383, 383,388,394,408,460,489,499,523,524,535,562,569,675,676, 748,778,786,797,955,968,1000,1245,1271,1420,1551,1694,2363,2754,2950,

Hi, I have a financial (zoo) time series with prices and volumes (although I can get the coredata as a matrix). Due to the data-source some indices have multiple observations. I want to aggregate these according to a weighted average. 11:00:01 34 1000 11:00:01 35 500 11:00:01 35 1000 11:00:02 34 500 11:00:02 35 500 should become 11:00:01 34.6 2500 11:00:02 34.5 1000 I currently do this

Dear R-help group, I have a dataset with daily closing prices from a stock exchange (consecutive 5 trading days) from a firm trading a specific commodity. The date variable looks like: quote_date 20080411 With the format; yyyymmdd. Moreover, I have another data set with a (average) weekly price of the underlying commodity. The date variables in this dataset are only year and a week number. I

Dear R-help list, I have grouped data, looking like this: cases <- c(23,12,56,81) total <- c(123,234,248,390) x1 <- c(0,0,1,1) x2 <- c(0,1,0,1) Data <- as.data.frame(cbind(cases,total,x1,x2)) Data I would like to run a logistic regression with group weights on these, where cases and (total-cases) are equal to the group weights (w). My final data would look

Hi, I'd like to blow up portions of my graph and put it in boxes beside the graph.Is there an addon to do this? -- Rajesh.J [[alternative HTML version deleted]]

Hello :) I'm just starting out with R and would appreciate your help with a couple of problems I am running into. I have used Sys.glob to get a list of all filenames having a particular file extension (in my case, *.txt) I would now like to use this list in the following manner: I would like to use each filename from the list and open that file into a tab separated matrix and proceed. How can

Dear R users, A have a dataframe (matrix) with two collumns (plot, and diameter (d)). I want all diameters values for different combination of plots. For example I want all d values for all posible combination, 100C2 (all d values for plot 1 with all d values in the plot 2.......with all d values from plot 1 with all d values from plot 100, ...... with all d values from plot 99 with all d values

Hello All: Using the below dataset how can I make a barplot with Date(X) and NumEggs(Y) by Site. Then plot Temp(lineplot) It seems really simple, but I am having a hard time trying to do it by Site. Thanks Date NumEggs Site Temp 1 2008-04-22 0 Massacre Flat (RK424.5) 51.20 2 2008-04-23 0 Massacre Flat (RK424.5) 50.80 3 2008-04-24

How I make the below matrix | a^0 0 0 ... 0 | | a^1 a^0 0 ... 0 | | a^2 a^1 a^0 ... 0 | | . | | . | | . | | a^n a^(n-1) a^(n-2) ... a^0 | with no loops, where "a" is a constant? ======================================================================== C?zar de Freitas Depto. de Estat?stica - UFPE Recife - Pernambuco - Brasil

Hello, I'm cc'ing R-Help. Sorry but your question was asked 3.5 years ago, I really don't remember it. Can you please post a question to R-Help, with a reproducible example that describes your problem? Rui Barradas ? Citando catalin roibu <catalinroibu at gmail.com>: > Dear Rui, > ? > I helped me some time ago with a code..... regarding aggregated data >

Drug US1 US2 Aptecha Celebrex 235.54 269.99 121.02 Detrol LA 157.99 190.99 55.3 Flomax 166.00 190.99 93.45 Lipitor 174.99 200.99 137.7 Novaldex 108.6 129.99 22.48 Norvasc 186.66 203.99 161.93 Plavix 107.99 106.99 64.53 Prevacid 117.39 134.99 59.83 Prilosec 115.99 126.99 57.75 Zyrtec 181.1 200.99 58.79 US1=c(235.54,157.99,166,174.99,108.60,186.66,107.99,117.39,115.99,181.10)

Dear R users, I want to aggregate all *d *data from all combination of n *plots* taken by k. Thank very much! My data is like that: plot d 1 14 1 13 1 12 1 14 1 18 1 20 1 21 1 43 1 108 1 43 2 41 2 61 2 83 2 61 2 84 2 45 2 21 2 12 2 11 ... 100 10 100 12 -- --- Catalin-Constantin ROIBU Forestry engineer, PhD Forestry Faculty of Suceava

> x<-as.POSIXct("1970-1-1", tz="UTC")-.5 > y<-as.POSIXct("1970-1-1", tz="UTC")+.5 > x==y [1] FALSE # of course but x and y "appear" to be the same when formatted, even with extra precision: > format(x, format="%Y-%m-%d %H:%M:%OS2") [1] "1970-01-01 00:00:00.50" > format(y, format="%Y-%m-%d

Hello, dear R-community. This is a question about TukeyHSD between factor combinations of a Three-Way ANOVA, which is - since it is a multi measure ANOVA - not a simple ANOVA but a Generalised Linear Mixed Model (GLMM), calculated with "lmer". > growth <- groupedData(length~meas|box_id,outer=~spec*comp*water,data=all.spec) > model <-

dum is a simple data frame transferred to Splus using the dump() command in Splus and the source() in R. All fields are numeric. There are no missing data. The data frame looks like it is should: > apply(dum,2,mode) yrcl sland s02 s234 "numeric" "numeric" "numeric" "numeric" > apply(dum,2,is.vector) yrcl sland s02 s234

All of our Kafka clusters are fairly write-heavy. The cluster in question is our second-heaviest ? we haven?t yet upgraded the heaviest, due to the issues we?ve been experiencing in this one. Here is an iostat example from a host within the same cluster, but without the RAID check running: [root at r2k1 ~] # iostat -xdmc 1 10 Linux 3.10.0-327.13.1.el7.x86_64 (r2k1) 05/27/16 _x86_64_ (32 CPU)

Hi, I am trying to plot the Case-Shiller index found at: http://www2.standardandpoors.com/spf/pdf/index/CSHomePrice_History_022445.xls The way I'm importing it into R is as follows: library(gdata) W <- read.xls("http://www2.standardandpoors.com/spf/pdf/index/CSHomePrice_History_022445.xls ", header=TRUE) attach(W) To give you and idea of what the data looks like: >

is there a way to use the actual index value for plotting zoo objects this is the way that the index is set up and a sample range of what I would like to plot 01/01/06 00:00:00 - 01/01/06 23:45:00 { library(zoo) # chron library(chron) fmt.chron <- function(x) { chron(sub(" .*", "", x), gsub(".* (.*)", "\\1:00", x)) }} x <- structure(c(15.57, 15.5,

Hello everybody! I have 3 gluster servers (*gluster 3.12.6, Centos 7.2*; those are actually virtual machines located on 3 separate physical XenServer7.1 servers) They are all connected via infiniband network. Iperf3 shows around *23 Gbit/s network bandwidth *between each 2 of them. Each server has 3 HDD put into a *stripe*3 thin pool (LVM2) *with logical volume created on top of it, formatted

Guess you went through user lists and tried something like this already http://lists.gluster.org/pipermail/gluster-users/2018-April/033811.html I have a same exact setup and below is as far as it went after months of trail and error. We all have somewhat same setup and same issue with this - you can find same post as yours on the daily basis. On Wed, Apr 11, 2018 at 3:03 PM, Anastasia Belyaeva