search for: sdat

glm(S ~ -1 + Mdif, family=quasipoisson(link=identity), start=strt, sdat) gives error: > Error in glm.fit(x = X, y = Y, weights = weights, start = start, etastart > = > etastart, : > cannot find valid starting values: please specify some strt is set to be the coefficient for a similar fit glm(S ~ -1 + I(Mdif + 1),... i.e. (Mdif + 1) is a vector si...

...e mean of yvar, diff and the p-value using tapply & by compared to a for loop. I cannot see what I am doing wrong. Can anyone help? > # Simulate some data - I'll do 2 simulations... > > xvar = rnorm(40, 20, 5) > yvar = rnorm(40, 22, 2) > num = factor(rep(1:2, each=20)) > sdat = data.frame(cbind(num, xvar, yvar)) > > # Define a function to do a simple t test and return some values... > > kindtest = function(varx, vary){ + res = t.test(varx, vary) + x.mn = res$estimate[1] + y.mn = res$estimate[2] + diff = y.mn-x.mn + pval = res$p.value + ca...

...the grand sum of the entries this function then calculates a 2D-discrete-ECDF. I keep thinking I am missing the obvious, but I did try searching. Here is my effort at creating that functionality: ecdf.tbl <- function (.dat) { .dat <- data.matrix(.dat) #speeds up calculations .sdat <- matrix(0, nrow(.dat), ncol(.dat) ) .sdat[] <- sapply(1:ncol(.dat), function(x) sapply(1:nrow(.dat), function(y) sum(.dat[1:y, 1:x], na.rm=TRUE ) ) ) return(.sdat) } > tst <- read.table(textConnection("NA 5 6 4 5...

Newbie here. Many apologies in advance for using the incorrect lingo. I'm new to statistics and VERY new to R. I have a "nx2" matrix , I want to sort the values based on the average of 2 columns and put k lowest (or highest) values in bin1, second k high/low values in bin2, and so on (bins would be of the same dimensions). I should also know what the first index (or position) of

...etter 'y'. For example, ? BA1y1 would equal 1, ?C3A2r3y1 would equal 1 and ?MA3r3y1r3 would equal 1. I currently have my code ?cycling through by the 5th character but my rows have variable length and the y can occur in several different locations in the row name. Here is my current code sdat <- read.csv (paste(data.dir,"/summary.data.csv", sep="")) year <-c("1", "2", "3") for (y in year) { sdat2 <- sdat[sapply(strsplit(as.character(sdat$GCM), ""), function(zzz)zzz[5] == y),] Thanks Daisy Englert Duursma Bioclimati...

...#39;m trying to fit nonlinear mixed effects model using nlme function but getting an error message. Here is what I have: fitted_model = nlme(scores~spline(b1,b2,b3,kt,time), fixed = list(b1~1, b2~1, b3~1, kt~1), random = b1+b2+b3~1, groups= ~id, data = sdat, start = c(b1=3.5,b2=2,b3=.60,kt=3.5),verbose=T) Error: Error in spline(b1, b2, b3, kt, time) : invalid interpolation method In addition: Warning message: In if (is.na(method)) stop("invalid interpolation method") : the condition has length > 1 and only the first elemen...

I used split() to split a variable by 3 years, but am wondering how to call up that split data and use it in further analyses. Can I make separate columns for the 3 resulting year groups? [[alternative HTML version deleted]]

...i/index.php?title=Export_maps_to_GE The lines that I can run are: data(meuse.grid) coordinates(meuse.grid) <- ~x+y gridded(meuse.grid) <- TRUE proj4string(meuse.grid) = CRS("+init=epsg:28992") # raster to polygon conversion; writeGDAL(meuse.grid["soil"], "meuse_soil.sdat", "SAGA", mvFlag=-99999) rsaga.geoprocessor(lib="shapes_grid", module=6, param=list(GRID="meuse_soil.sgrd", POLYGONS="meuse_soil.shp", CLASS_ALL=1)) and it crashes there with the following error message: Error in setwd(env$worksp...

...weights, jd, vp, ne, nx, nlam, flmin, ulam, thresh, isd, lmu = integer(1), a0 = double (nlam), ca = double(nx * nlam), ia = integer(nx), nin = integer(nlam), rsq = d ouble(nlam), alm = double(nlam), nlp = integer(1), jerr = integer(1), PACK AGE = "glmnet") 2: glmnet(xstd$sdat, ystd$sdat, family = c("gaussian"), lambda = lambda, stan dardize = FALSE, type = c("naive")) 3: sl_cv.glmnet(Xnew, ynew, K = 10, fit_all$lambda) aborting ... Thanks again! Dan Kamalic Director of Systems Analysis / Administration Boston University Department of Biomedic...

This is a bit updated set of chain.c32 changes that simplifies a few things (and in partiter part), fixes few minor issues and adds a few new features. Details are in the following commits, below is the summary and pull details at the end. Shao - any chance to peek over them ? Most of those are relatively simple changes and well tested, though of course something might have slipped my attention.

...f [ 41.333569] sd 44:0:0:0: [sdas] Mode Sense: 00 3a 00 00 [ 41.333569] sd 44:0:0:0: [sdas] Write cache: disabled, read cache: enabled, doesn't support DPO or FUA [ 41.333569] sdas: unknown partition table [ 41.346492] sd 44:0:0:0: [sdas] Attached SCSI disk [ 41.346492] sd 45:0:0:0: [sdat] 976773168 512-byte hardware sectors (500108 MB) [ 41.346492] sd 45:0:0:0: [sdat] Write Protect is off [ 41.346492] sd 45:0:0:0: [sdat] Mode Sense: 00 3a 00 00 [ 41.346492] sd 45:0:0:0: [sdat] Write cache: disabled, read cache: enabled, doesn't support DPO or FUA [ 41.346492] sd 45:0:0:...

Hi, I currently manage an environment where my Unix and Windows usernames are defined, but different. So, I must use a username map. But I have several Samba servers (one for each NFS package) and I must maintain all these username map files. As we want everybody to map the same way on every Samba servers in a NIS domain, we added the 'Windows user' information to the NIS

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

...BI.ySBAR 6000: 0a 05 00 0a 10 5b 81 3a 53 4d 42 49 01 48 53 54 .....[.:SMBI.HST 6010: 53 08 00 08 48 43 4f 4e 08 48 43 4f 4d 08 54 58 S...HCON.HCOM.TX 6020: 53 41 08 44 41 54 30 08 44 41 54 31 08 48 42 44 SA.DAT0.DAT1.HBD 6030: 52 08 50 45 43 52 08 52 58 53 41 08 53 44 41 54 R.PECR.RXSA.SDAT 6040: 10 14 44 04 53 53 58 42 0a a0 07 53 54 52 54 a4 ..D.SSXB...STRT. 6050: 00 70 00 49 32 43 45 70 0a bf 48 53 54 53 70 68 .p.I2CEp..HSTSph 6060: 54 58 53 41 70 69 48 43 4f 4d 70 0a 48 48 43 4f TXSApiHCOMp.HHCO 6070: 4e a0 12 43 4f 4d 50 7d 48 53 54 53 0a ff 48 53 N..COMP}HSTS..HS 6...