similar to: comparison of large data set

Hi everyone, I have received assistance from my previous post (comparison of large data set), but I still need some assistance. This is the revised code (along with some new code that I am requesting assistance with correcting): files <- paste0("MaxFloodDepth_", 1:54, ".txt") fileobs <- "MaxFloodDepth_Observed.txt" readfun <- function(x) { u <-

Hi, how are you? My name is Irucka Embry and I would like assistance on properly subsetting some data obtained from a .csv file. > file <-"info_n.csv" > INFO<- getMetaDataFromFile(file) INFO is 92 observations of 14 variables > INFO<- subset(INFO,site.no==02169570) INFO is 0 observations of 14 variables When the site.no is used I should have 1 observation of 14

Hi, how are you? I have 1000 total columns with 100 rows and I have shown 4 columns of 100 rows below. I want the zeros to be read in as numerical or integer values rather than a character. Right now I have a 100x0 character matrix rather than a 100x1000 double matrix. What do I need to do to read in a .csv spreadsheet with counting the zeros as a numerical or integer value? Thank-you. Irucka

Hi, I have the following function: getDataFromDVFileCustom <- function (file, hasHeader = TRUE, separator = "\t") { DVdatatmp <- as.matrix(read.table(file, sep = "\t", fill = TRUE, comment.char = "#", as.is = TRUE, stringsAsFactors = FALSE, na.strings = "NA")) DVdatatmper <- as.matrix(DVdatatmp[ , c("datetime",

Hi Irucka, I tried it and was able to plot it without any errors.? Here, your code indicates you need two lines. temper[[i]][1] ?temper[[1]][1] # which is the column 1. ? Month 1???? 1 2???? 2 3???? 3 ?temper[[1]][2] #? Data1 #1?? 1.5 #2? 12.3 #3? 11.4 Suppose I use names(temper) instead of seq_along(temper) pdf("irucka.pdf") ?lapply(names(temper),function(i)

HI Irucka, Please check this: temp<- structure(list(`:Bostoncitydata` = structure(list(Month = c(1L, 2L, 3L, NA), Data1 = c(1.5, 12.3, 11.4, NA), Data2 = c(9.1342, 12.31, 3.5, NA)), .Names = c("Month", "Data1", "Data2"), class = "data.frame", row.names = c(NA, -4L)), `:Chicagocitydata` = structure(list(Month = c(1L, 2L, 3L, 4L, 5L, NA), Data1 = c(1.52,

Hi, how are you? I am trying to replicate the binary data f(2) function in the attached document by starting with the simple example found below: observed <- matrix(c(0, 1, 0, 0, 1, 1, 1, 0, 0),3,3,byrow=TRUE) data <- matrix(c(1, 1, 0, 0, 1, 0, 0, 0, 1),3,3,byrow=TRUE) f2 = sum(probability of the matrix element where the matrix element is present in both the observed and the

Hi all, I have the following in my apt sources.list file: deb https:// /bin/linux/ubuntu trusty/ I issued the following command to obtain the public key: sudo apt-key adv --keyserver keyserver.ubuntu.com --recv-keys E084DAB9 I am receiving the following error message though after updating the Synaptic package manager: Failed to fetch

Hi everyone, now I am trying to finish writing the code (I had asked for assistance on subtracting arrays) This is what I what I am running in R: > source("/home/ie/Documents/TTU/GA_Research/GLUE/R-Project/R_GLUE_Example/NSEr.R") NSEr <- function (obs, sim) { {jjh <- (as.vector(obs) - sim)^2 Xjjhs <- apply(Xjjh, 2, sum) Yii <- (obs - mean(obs))^2 Yiis <- apply(Yii, 2,

Hi all, I have a couple of questions. 1) Is there a more effective way to remove the following pattern (patternrm) from the station_id_Timeseries (see below) patternrm <- c(_2000_4_OND, _2001_1_JFM, _2001_2_AMJ, _2001_3_JAS, _2001_4_OND, _2002_1_JFM, _2002_2_AMJ, _2002_3_JAS, _2002_4_OND, _2003_1_JFM, _2003_2_AMJ, _2003_3_JAS, _2003_4_OND, _2004_1_JFM, _2004_2_AMJ, _2004_3_JAS, _2000_4_OND) #

Hi all, I have attached the function file (revisedfunction) and the list of station IDs (StationIDs) in .pdf format. I have written a function to perform various operations on the station IDs in two groups. The first 6 station IDs require one less step than the remaining 68 station IDs. The file referenced at the beginning of the function is the file containing the StationIDs. > file =

Hi all, I have a set of 54 files that I need to convert from ASCII grid format to .shp files to .bnd files for BayesX. I have the following R code to operate on those files: library(maptools) library(Grid2Polygons) library(BayesX) library(BayesXsrc) library(R2BayesX) readfunct <- function(x) { u <- readAsciiGrid(x) } modfilesmore <- paste0("MaxFloodDepth_", 1:54,

Hi everyone, I am having trouble with creating a loop to subtract arrays. In R, this is what I have done: > Vobsr <- read.csv("Observed_Flow.csv", header = TRUE, sep =",") # see data > below > Vsimr <- read.csv("1000Samples_Vsim.csv", header = TRUE, sep =",") # see > data below > Vobsr <- as.matrix(Vobsr[,-1]) # remove column 1 from

Hi, let me correct my previous post ("zeros in double matrix rather than character matrix"). I have included entries from the .csv file that has 4 columns of 1000 columns and 101 rows. The first column with the "Time (day)" will not be removed, but the row headers are the numbers 1 to 1000 and they will remain in the matrix. In the previous post I mentioned that R makes the

I have code like this: - #----------------------------------------------------------------------- ------------------------------------------------------ x=scan() 0 0 0 0 0 1 2 3 4 y=scan() 1 1 1 2 2 1 3 4 5 plot(x,y) identify(0,1,3) #Allows me to select manually to identify co-ordinate (0,1) as being duplicated 3 times identify(0,2,2) #Allows me to select manually to identify co-ordinate

Hi, I've plotted a Kaplan-Meier curve but the curves only range from 0.7 to 1 on the y-axis. Therefore I have used: - ylim=c(0.7,1) [although I think convention dictates that you plot 0.5 to 1 to show the median? A few papers I've read have done this]?? BUT, I would like a symbol like // (but rotated 90 degrees) to indicate that the y-axis has been truncated. I'd need to be able to

Hi all, I'm wanting to be able to find and store the z-score of coxph below: - modz=coxph(Surv(TSURV,STATUS)~RAGE+DAGE+REG_WTIME_M+CLD_ISCH+POLY_VS, data=kidneyT,method="breslow") I know summary(modz) will give me this, but how do i extract the standard error or z-score values in a similar way to obtaining the coefficients by coef(modz) ? I think it must be something to do with

Dear R Help list, My question is regarding extracting the standard error or Z-score from a cph or coxph call. My Cox model is: - modz=cph(Surv(TSURV,STATUS)~RAGE+DAGE+REG_WTIME_M+CLD_ISCH+POLY_VS, data=kidneyT,method="breslow", x=T, y=T) I've used names(modz) but can't see anything that will let me extract the Z scores for each coefficient or the standard errors in the same

Hi Arun, thank-you for your assistance. I'm sorry that I did not provide a reproducible example. I will provide a reproducible example the next time. Using the next to last line of your code I was able to get what I needed without using the for loop. Thanks again. Irucka <-----Original Message-----> >From: arun kirshna [via R] [ml-node+s789695n4659154h50@n4.nabble.com]

Hi Arun, thank you. I will look through that later on today when I get a chance. I have to complete another part of this project now. Irucka <-----Original Message-----> >From: arun kirshna [via R] [ml-node+s789695n4661145h71@n4.nabble.com] >Sent: 3/12/2013 10:35:39 PM >To: iruckaE@mail2world.com >Subject: Re: how do I read certain files from a directory based on number of