similar to: R enquire

Hi, there I am a graduate student new to coding in S who is hitting a bit of a wall at present using an "optim" function. I am running into some troubles, and was hoping someone might be able to recognize where I am going wrong. As background: I have constructed a loop that carries out a 365-day calculation for a mass-balance model. Basically, the model depends on 2 variables (p,

Dear All, I have several sets of data such as this: year jday avg_m3s 1 1960 1 4.262307 2 1960 2 4.242308 3 1960 3 4.216923 4 1960 4 4.185385 5 1960 5 4.151538 6 1960 6 4.133846 ... There is a value for each day of multiple years. In this particular data set it goes up to 1974. I am am looking to obtain the minimum and maximum values for each year, but also know on which

hello all, this is my reproducible example data frame test<-structure(list(date = structure(c(1262300400, 1262304000, 1262304000, 1262307600, 1262307600, 1262311200, 1262311200, 1262314800, 1262314800, 1262318400, 1262318400, 1262322000, 1262322000, 1262325600, 1262325600, 1262329200, 1262329200, 1262332800, 1262332800, 1262336400, 1262336400, 1262340000, 1262340000, 1262343600,

Dear R users, Please help me with the following problem. I have a data set of the form ´dat´   dat <- data.frame(id = seq(1:5), trt=c(0,0,1,1,1),tid=c(0,0,0,0,0),                   ntid=c(0,0,0,0,0))   The function auto is used to generate the autoregressive model:   auto <- function(seed,delta,beta,maxt,dat){ set.seed(seed)  sp <- lapply(split(dat,dat$id),function(x){  while(x$tid <

Dear R users, Kindly help me to look into the following problem. I have a dataset like the following dat <- data.frame(id = seq(1:5), trt=c(0,0,1,1,1),tid=c(0,0,0,0,0),                   ntid=c(0,0,0,0,0)) I want to create waiting times using   B <- rexp(1) f1 <- (1/phi)+ ntid; f2 <- abs(1+phi*delta*tid*exp(beta*trt)) f3 <- phi*delta*exp(beta*trt) wait <-

I consider this to be an old bug, which also persists in Splus 7. It is unnecessary, and annoying. ## Section 1: Consider a simple data frame with three possible factors (in levels) d <- data.frame(a=sort(rnorm(10)*10), b=factor(c(rep("A",4), rep("C", 6)), levels=c("A","B","C"))) tapply(d$a, d$b, mean) # returns three results, which I

The issue is not with boxplot, but with split. boxplot.formula() calls boxplot(split(split(mf[[response]], mf[-response]), ...), but look at what split() returns when there are empty levels in the factor: > f <- factor(gl(3, 6), levels=1:5) > y <- rnorm(f) > split(y, f) $"1" [1] 0.4832124 1.1924811 0.3657797 1.7400198 0.5577356 0.9889520 $"2" [1] -1.1296642

Sorry if this is a stupid question but I tried reading, studying and finding solution but am not making it alone :) On my home lan I have a Linux (Debian sid) box (192.168.174.244) running 3.0.10 samba and cups 1.1.23. My objective is to print from this box to the other home PC which is running XP (192.168.174.242) and sharing an Epson R200 printer. Follows the troubleshooting I could do

Hi there I thought this would be of particular interest to people using 'optim' functions and perhaps people involved with R development. I've been beaten down by R trying to get it to perform an optimization on a mass-balance model. I've written the same program in excel, and using the 'solver' function, it comes up with an answer for my variables (p, ACT, which

To All How does one get a list of functions from a loaded package so that one can then get the appropriate help for each of the functions. Currently my method is based on a lot of trial-and-error. Here's an example of what I mean... >From this forum I learn that an interesting package called "multtest" exists on Bioconductor. I then use R Console's "Packages" --

Dear friends of lme, After so many year with lme, I feel ashamed that I cannot get this to work. Maybe it's a syntax problem, but possibly a lack of understanding. We have growth curves of new dental bone that can well be modeled by a linear growth curve, for two different treatments and several subjects as random parameter. By definition, newbone is zero at t=0, so I tried to force the

I just tried setting: performance.parallel-readdir on features.cache-invalidation on features.cache-invalidation-timeout 600 performance.stat-prefetch performance.cache-invalidation performance.md-cache-timeout 600 network.inode-lru-limit 50000 performance.cache-invalidation on and clients could not see their files with ls when accessing via a fuse mount. The files and directories were there,

2017-10-10 8:25 GMT+02:00 Karan Sandha <ksandha at redhat.com>: > Hi Gandalf, > > We have multiple tuning to do for small-files which decrease the time for > negative lookups , meta-data caching, parallel readdir. Bumping the server > and client event threads will help you out in increasing the small file > performance. > > gluster v set <vol-name> group

You can get help on the whole package by > help(package="multtest") which is likely pretty close to what you want. There's the index etc for the package on the web as well. You can also just look in the package's installation directory. If it's loaded you can do an ls(2) say if it loaded in position 2, to get all objects in the package. Reid -----Original Message-----

The following works fine as a standalone ruby program: ********************************************************************************************** t = Time.now starter = Time.local(t.year,t.mon, t.day) + (24 *3600) ender = starter + 7 * 24 *3600 lat = 40.352039 lon = -74.191961 require ''soap/wsdlDriver'' params = {:maxt => false, :mint => false, :temp => true, :dew

Dear R helpers, I have what another member on this forum described as an embarrassingly parallel problem. I am trying to fit models on subsets of some data based on unique combinations of two id factors in the dataset. Total number of combinations is 30^5, and this takes a long time. So, I would like fit models for each of the datasets produced by subsetting on the unique combinations, splitting

Hi everybody I have a questionnaire with a lot of questions that allow for more than one option to be chosen (like a tickbox in a html form). The data captured on a mobile device and is supplied in a format where every option is a separate variable (logical). I want to develop a generic function to process these questions. As part of the analysis I want make a matrix of the p-value from the

hi, we recently upgraded from samba 2.2.8a w/ldap to samba3 w/ldap and having troubles with domain group mappings. first symptom is that net groupmap list returns nothing. [root@lola root]# net groupmap list [root@lola root]# if we try to add a groupmapping, the command fails: net groupmap add sid=<sid-gid> ntgroup="Domain Admins" unixgroup=domadm type=domain --debuglevel=4

Hi List Apologies if this isn''t the correct place for this query (I''ve tried a search of the mail archives but not had much joy). I''m running R (2.14.0) on a Mac (OSX v 10.5.8, 2.66GHz, 4GB memory) and am having a few performance issues with reading data in from a Postres database (using RPostgreSQL). My query / code are as below # -----------------------------

2009/1/22 Spiros Ioannou <sivann at gmail.com> > Cher Arnaud, > Hi Spiros, > I send you a demo that reads info from an Emerson Liebert > GXT2-10000T230 through the RS-232 port. This uses the "ESP-II (ESP 2) > Protocol". > Current supported data includes: measurements (voltage, temperature, > etc) , ID strings (serial number, etc), status flags (inverter on,