similar to: "unvector" ?

Hi everybody, recently I had to teach a course on Cox model, of which I am not a specialist, to an audience of medical epidemiologists. Not a good idea you might say.. anyway, someone in the audience was very hostile. At some point, he sayed that Cox model was useless, since all you have to do is count who dies and who survives, divide by the sample sizes and compute a relative risk, and if there

dear R-users, I have to model counts where all counts above some threshold have been censored. In the same dataset I have too many zeroes for a Poisson or even a negative binomial distribution to make sense, so I would need a zero-inflated-censored negative binomial family for use in glmer (or glmmADMB?). That seems not to exist. my question is : how could I add a custom-built family of

Hi, a colleague ran a stepwise discriminant analysis twice in a row and got different results, suggesting some "sochasticity" in the algorithms involved. I looked at her data and found that there was a lot of collinearity, so that I reckoned that maybe "stepclass" (klaR) cannot find a clear winner when trying to include a new variable and makes a random choice. Is that true?

He everybody, I want to add 1 to some elements of a vector: x is a vector u is a vector of idices, that is, integers assumed to be within the range 1..length(x) and I want to add 1 to the elements of x each time their index appears in u x[u]<-x[u]+1 works only when there are no duplicated values in u I found this solution: tu <- table(u) indices <- as.numeric(names(tu)) x[indices]

Hi everybody, I have to run under Ubuntu a programs repeatedly with different arguments and I am using R just to generate the data files and call the external program. basically, in my script I have inside a loop these two lines: command <- paste(<some constant or variable strings>,sep="") system(command,intern=T,wait=T) when I run this script, I get a number of warnings,

Dear R users, I'm trying to use function 'refit' from lme4 and I get this error that I can't understand: > refit(dolo4.model4,cbind(uu,50-uu)) Error in function (classes, fdef, mtable) : unable to find an inherited method for function "refit", for signature "mer", "matrix" if I try: > refit(dolo4.model4,uu) Error in asMethod(object) :

Hi listers, I want to make lattice plots xyplots with the indication of legends inside each panel with only the points and the lines actually ploted inside each given panel according to the group(ing) factor. The code below shows what I have achieved so far and I hope will make clear what I want to have. It seems to me that my solution is a very "dirty hack" and there certainly is

Dear All, After followed the help of nnet, I could get the networks trained and, excitedly, get the prediction for other samples. It is a two classes data set, I used "N" and "P" to label the two. My question is, how do I get the predicted numerical value for each sample? Not just give me the label(either "N" or "P")? Thanks! FYI: The nnet example I

Dear R'helpers, I am reading lines in a .txt file Each line is stocked into a n elements object, as this: [958] " 422 287 339 31 203 602 547 1026 500 366 346 227" [959] " 410 67 11 220 110 451 562 598 732 163 163 220" [960] " 179 513 95 186 102 595 333 1289 804 210 294 459" [961] " 276 153 307 138 126 233 623 739 521 421 209 75" [962] " 64

Hello !! II'm trying to test for my fixed effects using an lmer with quasipoisson errors. Since my lmer model is corrected for overdispersion using this kind of errors, I should use during model simplification in my Anovas *F test *and not *Chi square test* to compare two models. So I write: > anova(model,model2,test="F") but R keeps performing a Chi square instead of the F

I'm looking for a tutorial or notes on the use of contrasts factor in linear model in R, I've found some mails and infos about in various documents about R, but I've probably missed a good review on this subject. -- Robert Espesser Laboratoire Parole et Langage UMR 6057, CNRS 29 Av. Robert Schuman 13621 AIX (FRANCE) Tel: +33 (0)4 42 95 36 26 Fax: +33 (0)4 42 59 50

Hi, I am trying to install pdftools package on R 3.3.1 (Ubuntu 16.04), but the following issue occurs: ------------------------- ANTICONF ERROR --------------------------- Configuration failed because poppler-cpp was not found. Try installing: * deb: libpoppler-cpp-dev (Debian, Ubuntu, etc) * rpm: poppler-cpp-devel (Fedora, CentOS, RHEL) * csw: poppler_dev (Solaris) * brew: poppler (Mac OSX)

Hi everybody, I just tried to run R on one of my projects but it did not want to run: R version 2.12.1 (2010-12-16) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-pc-linux-gnu (64-bit) <some lines deleted> Loading required package: utils Error in loadNamespace(i[[1L]], c(lib.loc, .libPaths())) : there is no package called

Dear all, I would like to solve a trivial problem (I guess it is) but can't find the right way. Maybe someone can help me ? I've got a table with two factors (station = station ID, buffer = buffer size in meters) and a value for each unique combination of those two factors (S = number of habitats within each buffer around each station) like this: TABLE 1 station buffer S Abaia01 200 2

Dear All I have problem with this: > IWJR.temp <-IWJR.missing[sample(1:length(IWJR.missing),1),] Error: subscript out of bounds How I can solved this. Thanks IRD

Hi Francois, Thanks for your quick response. Actually, I had already done that... sudo apt-get install libpoppler-cpp-dev Reading package lists... Done Building dependency tree Reading state information... Done libpoppler-cpp-dev is already the newest version (0.41.0-0ubuntu1.1). 0 upgraded, 0 newly installed, 0 to remove and 107 not upgraded. Therefore, I assume I have this installed. Best

Dear R'helpers, I guess the solution is trivial but despite some hours of testing and looking at the doc, I was unable to fix the following problem. I want to do either: for (i in 1:lat) a[i,which(a[i,]==0] <- NA or: for (i in 1:lat) a[i,which(a[i,]==-9999|a[i,]==9999|a[i,]==9998|a[i,] ==-9998)] <- NA if the word "Time" is or is not present in variable

Greetings, I am having trouble calculating artificial neural network misclassification errors using errorest() from the ipred package. I have had no problems estimating the values with randomForest() or svm(), but can't seem to get it to work with nnet(). I believe this is due to the output of the predict.nnet() function within cv.factor(). Below is a quick example of the problem I'm

Hi! I am trying to install in my R environment the rhdf5 package and library but it seems to have vanished from either the CRAN or BioConductors sites. Can you tell me where it would be possible to find it or any R library (or function) able to read hdf files? Sincerely, Nicolas Degallier UMR 7159 / IRD UR182 Laboratoire d'Oc?anographie et du Climat, Exp?rimentation et Approches

I'm running asterisk 16 on Fedora 31. If I start asterisk as user asterisk, all goes well. But if I start asterisk from systemd: asterisk[1411]: [Jan 21 19:36:47] ERROR[1411]: res_sorcery_config.c:320 sorcery_config_internal_load: Unable to load config file 'pjsip.conf' Jan 21 19:36:47 asterisk.riverside asterisk[1411]: [Jan 21 19:36:47] ERROR[1411]: config_options.c:710