similar to: modifynig some elements of a vector

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?

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 ev'rybody, is there a way to pass a vector to a function expecting separate arguments? more specifically, I have a character vector, say u and I want a single string, but >paste(u) doesn't work, so I would like something like >paste(unvector(u)). I am interested in a solution to the general problem too, as the only one I found is maintaining two versions of the functions I

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

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,

Hi everyone, I determined the presence of three types parasites in a passerine bird over two years. I would like to create a bar chart that shows the proportion infected on the y and year/parasite on the x such that each type of parasite is grouped together (single label) and a bar for each year . This would show if there have been changes in the prevalence of a the parasite over two years.

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 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 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, I would like some extra information on the 'cgh' package in R. I noticed that there isn't much activity regarding this package on the R and BioC mailing list (I googled it). I started using this package and I have few questions: 1/ As I have a custom tiling like array @8um features resolution (affy), I have a lot of probes to work with. I'm assuming it is correct to

Hello, I have a blade server with SAS HDD's of 12TB in total. 3 HDD's of 4TB each. Is it possible to install Cent OS 7.8.2003 on 12TB disk space? I will be installing Cent OS on the bare metal HW. I referred = https://wiki.centos.org/About/Product But slightly confused with the 'maximum file size' row for ext4 FS. Thanks & Regards, Amey.

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

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 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

I have a simple script: #!/bin/bash # index=0 total=0 names=() ip=() while read -r LINE do NODENAME=` echo $LINE | cut -f 1 -d ','` IP=` echo $LINE | cut -f 2 -d ','` names[index]="$NODENAME" ip[index]="$IP" index=`expr index+1` total=`expr total+1` done <<< $(cat list.txt) simple file: more list.txt name1,ip1 name2,ip2 name3,ip3 output when

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