search for: prcc

Hi everybody, I'm trying to create a numerical data frame on which to perform PRCC. So far I have created a data frame that consists of function/vector output that displays in numerical form, but when I try and run PRCC (from epiR package) I get the following error message: "Error in solve.default(C) : Lapack routine dgesv: system is exactly singular" It appears th...

...ess = FALSE, simple = TRUE) and it returned a single series of 500 samples, rather than a series of 500 samples per function. I ultimately need a sample of each variable function that I can run through my model, putting each individual variable function as a constant instead, and then performing PRCC. Is there anyone who can advise on how to do this, or failing that, where I should look for sample code? Thank you for any help you are able to give, Aimee.

...o learn how to use R, but are not supposed to know/understand what maths and stats are under the bonnet of the test (I don't run the course, no comment). So the 'plan' is to give them a lecture about models in epidemiology and sensitivity, and then the punch line should be: 'and the PRCC can be asily done in R using the function ???'. Cheers, Federico [1] S.M. Blower and H. Dowlatabadi, 1994. SENSITIVITY AND UNCERTAINTY ANALYSIS OF COMPLEX-MODELS OF DISEASE TRANSMISSION - AN HIV MODEL, AS AN EXAMPLE. International Statistical Review 62(2) 229-243. -- Federico C. F. Calb...

...4.876022801 4.787737 19873.134 18660.02 9 39.99782 0.4109272317 0.6606016486 0.191627831 5.625588 11086.803 13569.30 I have produced a column (r1) where each cell is a vector of the previous cells in that row. I tried binding them with zz<-cbind(samplevalues, r1) and then running epi<-epi.prcc(zz, sided.test=2) which is what produced the original error message of "Error in solve.default(C) : Lapack routine dgesv: system is exactly singular" I then tried 'sensitivity' with the syntax zz<-pcc(samplevalues, r1). I didn't get any error messages, but print(zz) onl...

Dear R users, I have been trying to estimate global sensitivity indices such as the sobol 1st and 2nd order indices. I managed to obtain the PRCC. The example presented in the sensitivity package on sobol2002 seems to work well for linear models: for example: calculate y for given x values. However, when trying to apply this technique to dynamic models (SIR type), the error messages just keep spinning from one angle to another. Even the d...