similar to: F test in lmer quasipoisson

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?

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

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 all, I was used to see an htlm documentation when writing the command lines ?plot help.search("plot") The configuration was changed on my computer and help appears now as a R console window, without link with internet. Where can I change the configuration to connect with the html help? Thanks, Sophie Bouchet -- Sophie Bouchet Post-doc INRA UMR de Génétique Végétale équipe

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

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

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

Dear listers, The following command (derived from the example in the ?stars help page) works : data(mtcars) stars(mtcars[, 1:7]) But the following gives an error: stars(mtcars[1, 1:7]) Error in s.y[i, ] : incorrect number of dimensions I was expecting to have the star graph for the first line (Mazda Rx4) The following give an incorrect graph for the first two cars : stars(mtcars[1:2, 1:7])

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

Hi listers, A few days ago I posted a question about the use of the stars function on selected lines of a frame. Thanks to two helpers, a closer look at the scale argument allowed to partially solve the problem. Yet I still have a problem with stars. Allow me to explain what I intend to do (sorry for my poor English and the long post): I want to graph an activity index of a fish during the day

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

Hi, we are having a problem with NFS using RDMA protocol over our FDR10 Infiniband network. I previously wrote to NFS mailing list about this, so you may find our discussion there. I have taken some load off the server which was using NFS for backups, and converted it to use SSH, but we are still having critical problems with NFS clients losing connection to the server, causing the clients

> > Here is a bit of an odd request... II'm in a large company with lots of > politics getting in the way of things. We were able to convince the main IT > team to let us take over a subnet and run our own DHCP server. This lets us > create a really nice PXE environment where we can build servers and ship > them out. However, now we are in need to rebuild servers already

Dear all, I've just downloaded R.12.0 from CRAN (http://cran.r-project.org/bin/windows/base/release.htm). I am unable to find RGUI.exe as well as Rterm.exe. Conversely to what is explained in readme.R.2.12.0: There are two versions of the R executable in R-2.12.0\bin\i386 (32-bit) or R-2.12.0\bin\x64 (64-bit). I can't find directories i386 or x64 in R-2.12.0\bin -- Renaud Lancelot EDEN

Bonjour, I'm running samba 2.0.4 on a Sparc Solaris 2.5.1 and I have some error logs when I'm trying to connect to a sharing directory from *some* NT4 WS (SP4 or SP5): pchp2 (194.xxx.xxx.xxx) connect to service centreAS as user truc (uid=1326, gid=1320) (pid 22939) [1999/08/26 15:52:57, 0] lib/util_sock.c:write_data(415) write_data: write failure. Error = Broken pipe [1999/08/26

I met today a computer crash and our maintenance officer had to reinstall some components of the OS (MS Windows XP Pro) as well as the Internet browser (among other things). Now, I cannot install packages from a distant repository: > utils:::menuInstallPkgs() Error in .readRDS(pfile) : unknown input format > traceback() 5: .readRDS(pfile) 4: .packages(all.available = TRUE) 3: