search for: ntimes

...aining.dat, Income>10) cs.training.dat_sub2 <- subset(cs.training.dat_sub1, Dept.Ratio<=1) cs.training.dat_sub3 <- subset(cs.training.dat_sub2, Credit.Limit.Ratio<=1) for (i in 1:500){ y.2 <- cs.training.dat_sub3$y[1+[i]:300+[i]] y.1 <- cs.training.dat_sub3$y[-(1+[i]:300+[i])] NTimes.60DaysLate.2 <- (cs.training.dat_sub3$NTimes.60DaysLate[1+[i]:300+[i]]) NTimes.60DaysLate.1 <- (cs.training.dat_sub3$NTimes.60DaysLate[-(1+[i]:300+[i])]) Credit.Limit.Ratio.2 <- (cs.training.dat_sub3$Credit.Limit.Ratio[1+[i]:300+[i]]) Credit.Limit.Ratio.1 <- (cs.training.dat_sub3$Credit...

...in() are much slower for the variables of type 'dates'. The difference gets bigger the larger the 'id' var is - but see this sample code: dts <- dates(c("02/27/92", "02/27/92", "01/14/92", "02/28/92", "02/01/92")) ntimes <- 700000 dts <- data.frame(rep(c(1:40), ntimes/8), chron(rep(dts, ntimes), format = c(dates = "m/d/y")), rep(c(0.123, 0.245, 0.423, 0.634, 0.256), ntimes)) names(dts) <- c("id", "date", "tbs") date() dat.1st &l...

...invalid number of copies in "rep". This happens no matter what values that I give the parameters when I call the function. Any suggestions on clearing this up would be appreciated! The function is printed below. Many Thanks, Monnie McGee mcmcmd <- function(d, w, mean, phi2, varwt, ntimes) { # A function to run a multi dimensional MCMC for a bivariate # normal distribution # d = dimension of problem # w = weight assigned to larger of two modes in target distribution # mu = mean of second mode # phi2 = multiplier of variance factor 2.38 # varwt = weight (0 < varweight <= 1)...

...get [1] 1423.001 note also that if v starts getting large, evaluation of F also gives the error. e.g. F(1, 713.5) any help? seems like looping works nicely here F2=function(z,v) { #check that the v is correct! ie -.5, .5, 3/2, 5/2, ..... if (v==-0.5) { return(1) }else { ntimes = v+1.5 R=array(0,dim=ntimes) R[1]=1 for (i in 2:ntimes) { R[i]= 2*v/z + 1/R[i-1] } return(R[ntimes]) } } Allan Clark ======== Lecturer in Statistical Sciences Department University of Cape Town 7701 Rondebosch South Africa TEL (Office): +27-21-650-3228 F...

I have a combined date and time. I would like to separate them out into two columns so I can do things such as take the mean by time across all dates. meas<-runif(435) nTime<-seq(1303975800, 1304757000, 1800) nDateT<-as.POSIXct(nTime, origin="1970-01-01") mat1<-cbind(nDateT,meas) means1<- aggregate(mat1$meas, list(nDateT), mean) This doesn't do anything as each day

...get [1] 1423.001 note also that if v starts getting large, evaluation of F also gives the error. e.g. F(1, 713.5) any help? seems like looping works nicely here F2=function(z,v) { #check that the v is correct! ie -.5, .5, 3/2, 5/2, ..... if (v==-0.5) { return(1) }else { ntimes = v+1.5 R=array(0,dim=ntimes) R[1]=1 for (i in 2:ntimes) { R[i]= 2*v/z + 1/R[i-1] } return(R[ntimes]) } } Allan Clark ======== Lecturer in Statistical Sciences Department University of Cape Town 7701 Rondebosch South Africa TEL (Office): +27-21-650-3228 F...

Apologies for any obtuseness in the following. We have been working on Version 2.0 of the randomSurvivalForest CRAN package and we're encountering a perplexing 'memory not mapped' segfault that we believe is "influenced" by GC. We essentially have two R functions, rsf.default(..), and predict.rsf(..) and two corresponding entry points, rsfGrow(...), and rsfPredict(...),

Hi I'm encountering an error with some of my Windows clients when using CTDB and samba to access a GPFS hosted share. We see the following error all the time. [2011/03/14 12:31:16.549084, 1] modules/vfs_gpfs.c:1099(vfs_gpfs_ntimes) vfs_gpfs_ntimes: set GPFS ntimes failed -1 The copy continues but is slow. All the timestamps seem reasonable for the data copied. The GPFS version is 3.3.0-9 and samba is 3.5.6 compiled from source with config line. $ ./configure --with-ctdb=/usr/local/src/ctdb --with-cluster-support --...

Hi List, I am building a package for a client to help them create and perform analyses against netcdf files which contain 'a temporal stack' of grids. For my examples and test cases, I create an example dataset in code (as this is a lot more space efficient than providing raw data). The code creates a netcdf file in tempdir() and an object of class 'ncdf' in the global namespace.

...rmat_t old_audio_fmt = { 0, 0, 0, 0, 0 }; + int status_interval = 0; audio_format_t new_audio_fmt; audio_reopen_arg_t *reopen_arg; + + decoder_callbacks_t decoder_callbacks; + void *decoder_callbacks_arg; - /* Flags and counters galore */ - int eof = 0, eos = 0, ret; - int nthc = 0, ntimesc = 0; - int next_status = 0; - static int status_interval = 0; - - /* Reset all of the signal flags */ - sig_request.cancel = 0; - sig_request.skipfile = 0; - sig_request.exit = 0; - sig_request.pause = 0; - + char *convbuffer; + int convsize = AUDIO_CHUNK_SIZE; + + convbuffer...

...distribute any loop. TSVC <https://github.com/shantanuatiith/TSVC_>: There are 151 loops coded in plain ‘C’, none of them got distributed. TSVC has particularly candidates valid for distribution like the one below. The inner loop in this example can be distributed. for (int nl = 0; nl < ntimes/2; nl++) { for (int i = 1; i < LEN; i++) { a[i] += c[i] * d[i]; b[i] = b[i - 1] + a[i] + d[i]; } dummy(a, b, c, d, e, aa, bb, cc, 0.); } MiBench <http://vhosts.eecs.umich.edu/mibench//>: There are 6539 loops in MiBench. *None* of the loops were distributed by the loop-distribu...

Dear All, since years I am struggling with Surv objects in data.frames. The following seems to have to do with it. See below the modified example from the help page of survSplit. The original works, as expected. If, however, a Surv object is added to the data.frame, each record gets doubled. Is there some solution other than avoiding Surv objects in data.frames? Thanks, Heinz

Hola a todos, Me gustaria saber si en un modelo mixto se puede usar el tamaño de muestra como peso o es incorrecto hacerlo. Por ejemplo en el comando 'lmer' del paquete lme4 hay una opcion 'weights' (igual que en 'lm' de stats). Si tengo datos de una medida (por ejemplo el peso) de 10 especies de aves (rango, de 50 a 3000 datos segun la especie) ¿puedo usar esta n en

...man page for more info.\n" diff -Naur vorbis-tools-1.0.orig/ogg123/ogg123.c vorbis-tools-1.0/ogg123/ogg123.c --- vorbis-tools-1.0.orig/ogg123/ogg123.c Sat Jul 6 21:12:18 2002 +++ vorbis-tools-1.0/ogg123/ogg123.c Sun Mar 9 16:59:39 2003 @@ -140,6 +140,7 @@ opts->nth = 1; opts->ntimes = 1; opts->seekpos = 0.0; + opts->endpos = 0.0; opts->buffer_size = 128 * 1024; opts->prebuffer = 0.0f; opts->input_buffer_size = 64 * 1024; @@ -250,6 +251,18 @@ print_statistics_action(NULL, pstats_arg); } +double currenttime (stat_format_t *stat_format, + da...

I have a large data frame, 77 rows, with 10 columns. Each row represents a unique individual with 10 characteristics, some of which are categorical factors and some continuous numeric variables. Each of the ten variables is important (the 10 columns obviously correspond to the individuals of interest). Importantly, this data set represents a population (not sample) of people with a certain medical

...; > TSVC <https://github.com/shantanuatiith/TSVC_>: > There are 151 loops coded in plain ‘C’, none of them got distributed. TSVC has particularly candidates valid for distribution like the one below. The inner loop in this example can be distributed. > > for (int nl = 0; nl < ntimes/2; nl++) { > for (int i = 1; i < LEN; i++) { > a[i] += c[i] * d[i]; > b[i] = b[i - 1] + a[i] + d[i]; > } > dummy(a, b, c, d, e, aa, bb, cc, 0.); > } > > MiBench <http://vhosts.eecs.umich.edu/mibench//>: > There are 6539 loops in MiBench. None of the loops...

Hi, I have around 300 files with the x, y and z coordinates of a DEM that I would like to join in a single data frame object. I know how to automate the import of the 300 data frames. in Bash ls > names in R names<-scan(names...) With rep() and data.frame() construct a series of read.table() commands, which you write to a file and execute them with source(). I do not know however how

Hi There: I was wondering if there is a way to read FORTRAN list data (similar to IDL's readf function). I often use FORTRAN for most of my number crunching, and use something like IDL to visualize and perform statistical analysis on that data. Since the each file is rather large (>100 Mb), formatting the output into columns or tables is impractical, hence the "list"

http://bugzilla.mindrot.org/show_bug.cgi?id=1255 Summary: Solaris contract support has a problem Product: Portable OpenSSH Version: 4.4p1 Platform: All OS/Version: Solaris Status: NEW Severity: normal Priority: P2 Component: sshd AssignedTo: bitbucket at mindrot.org ReportedBy:

...models to data, generated from a negative binomial distribution : count ~ Neg.Bin (mu, theta). Here is the code: ============================================================================== #Generate the data frame set.seed(2153) mydf<-data.frame(subjs=seq(1:nsubjids), counts=rnbinom(nsubjids*ntimes, size=0.5, mu=1.8)) #Model require(MASS); require(lme4) summary(glm.nb(counts~1, data=mydf)) summary(lmer(counts~1+(1|subjs), family=negative.binomial(theta=fixed.nb0$theta), data=mydf)) ============================================================================== The glm.nb output gives : mu=ex...