similar to: Issues running sensitivity analysis

Dear all mailing listers, Does Anyone have the same problem as mine when using the fast99 (extended-FAST method) to perform SA of model with norm distribution inputs? See the simple example given following. Any suggestion will be greatly appreciated. Thank you! Marino # Simple example # 1. uniform version (It works well) library(sensitivity) Myfun<-function(x){return(rowSums(x))} SA1

Hi! I am trying to implement the Morris One at a Time (OAT) sensitivity analysis technique in R using the package 'sensitivity'. The OAT and similar other techniques in the package requires data frame X (which stores the design or parameter combinations) to be populated by some function in the package such as fast99 or morris. These functions take only the range and step size of the

Hi All, the new release v1.4-0 of the "sensitivity" package is now available on CRAN The "sensitivity" package is devoted to factor screening and global sensitivity analysis of numerical model output. Here are the new features list: New functionalities: * sequential bifurcation method (sb function) [VERSION ALPHA] * 3D plot 3D for the Morris method (function

Hi All, the new release v1.4-0 of the "sensitivity" package is now available on CRAN The "sensitivity" package is devoted to factor screening and global sensitivity analysis of numerical model output. Here are the new features list: New functionalities: * sequential bifurcation method (sb function) [VERSION ALPHA] * 3D plot 3D for the Morris method (function

Hi, I'd like to use Latin Hypercube Sampling (LHC) in the the context of uncertainty / sensitivity analysis of a complex model with approximately 10 input variables. With the LHC approach I'd like to generate parameter combinations for my model input variables. Therefore I came across an simple example here on the mailing list (

Hi, I have a dataset with a continuous response variable and, among other predictors, an ordinal variable. Here is what it could look like treatment <- factor(rep(c("AA", "AC", "AD","AE", "AB"), each = 10)) length <- c(75, 67, 70, 75, 65, 71, 67, 67, 76, 68, 57, 58, 60, 59, 62, 60, 60, 57, 59, 61, 58,

(1) ...need to speed up a monte-carlo sampling...any suggestions about how I can get R to use all 8 cores of a mac pro would be most useful and very appreciated... (2) spent the last few hours trying to get pnmath to compile under os- x 10.5.4... using gcc version 4.2.1 (Apple Inc. build 5553) as downloaded from CRAN, xcode 3.0... ...xcode 3.1 installed over top of above after

Full_Name: Rob Carnell Version: 2.3.0 OS: Windows XP sp 2 Submission from: (NULL) (131.167.72.220) Thank you for all you work on R project! I noticed the following behavior today. It is what I believe to be an inconsistency. > qunif(.5, 0, Inf) [1] Inf > qunif(.5, -Inf, 1) [1] NaN Warning message: NaNs produced in: qunif(p, min, max, lower.tail, log.p) I believe they should yield the

Hi, Is the following function built in somewhere? concat = function(v) { res = "" for (i in 1:length(v)) res = paste(res,v[i],sep="") res } e.g. > concat(c("12","3","45")) [1] "12345" Cheers, Andy -- Andy Fugard, Postgraduate Research Student Psychology (Room F3), The University of Edinburgh, 7 George

Dear All, I am working on a dataset having the dependent variable as ordinal data(discrete data) and multiple independent variables. I need to find the likelihood for the NULL model.i.e the model with only the dependent variable and all other independent variables as zero. Kindly let me know how to find the likelihood for a NULL model in R. Is there any specific function in R that can do

I'm trying something that I thought would be pretty simple, but it's proving quite frustrating... I want to display, for instance, the correlation coefficient "rho" in a graph. I can do something like: text(x, y, paste("rho =", cor)) where cor would be my previously calculated correlation coefficient, and x and y the coordinates. Obviously that displays

The Runuran package includes a simulation for the Triangular Distribution, urtriang. But I didn't find an analogue of the q* functions (rnorm - qnorm, runif - qunif), that would invert the CDF (?uqtriang). Are there any such functions? Alberto Monteiro

I think you should take this conversation private or seek local statistical expertise. This is about strategies for analysis, not about programming in R. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, May 31, 2017 at

Hi, Does anyone know an easy way of retaining rows in a model.matrix where missing values are present in the predictors. Ideally I'd be able to retain these rows as zeros. Thanks, Jarrod -- The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336.

Hi there, Am I correct to believe that tetrachoric correlation is a special case of polychoric correlation when there are only two levels to the ordered factor? Thus it should be okay to use hetcor from the polycor package to build a matrix of correlations for binary variables? If this is true, how can one estimate 95% confidence intervals for the correlations? My guess would be mat =

The (2.2.1) upshid-ups driver reports battery.temperature rather than ups.temperature. Here's a quick patch to have it show up in the status page... --- clients/upsstats.c.DIST 2008-03-04 09:15:40.000000000 +0000 +++ clients/upsstats.c 2008-03-04 09:16:20.000000000 +0000 @@ -787,6 +787,11 @@ return 1; } + if (!strcmp(cmd, "BATTTEMP")) { +

On Fri, Jun 23, 2017 at 11:26 AM, Chris Cooke <cc at inf.ed.ac.uk> wrote: > qemu-kvm-ev-2.6.0-28.el7.10.1 > <https://cbs.centos.org/koji/buildinfo?buildID=17495> has been tagged for > release and will soon be available on CentOS mirrors. > > > Hi Sandro, > > Thanks. That looks like a typo in the package release field? Should be > 28.el7_3.10.1, rather than

Hello, I'm using geom_hline to add a minimum line to my plot (representing the best solution found so far by a search algorithm). I'd like to annotate this line with it's numerical value to save trying to read it off the graph, but I can't see a clear way to do this - any ideas? (Alternatively, if this is against the spirit of the grammar of graphics, is there a better

Thank you very much Rob for your answer. I have some difficulties to understand how to apply my agent-based model to each parameter combination generated by the LHS, in particular when parameters are defined by probability distributions. Indeed, I have multiple parameters in my model: parameters which are defined by a single value (like ?temperature", "pressure?) and parameters which are

Dear list, is there a way to open a .zip folder so that one can extract and modify files inside and then save teh .zip folder again? thanks! Jose -- Dr. Jose I. de las Heras Email: J.delasHeras at ed.ac.uk The Wellcome Trust Centre for Cell Biology Phone: +44 (0)131 6507095 Institute for Cell & Molecular Biology Fax: +44 (0)131 6507360 Swann Building,