search for: fems

Dear all, I would like to know if there is a function for sorting a data.frame based on one column. sort() does it for a single vector but not for the whole data. order() does it if there are relations between the columns but in my case it is simply a data.frame. Further, I would like to plot other data, pop.size<- c(800, 800, 1500, 1000, 158, 300, 740, 250, 2000, 1500, 250, 700)

Dear all, I'm using zeroinfl() from the pscl-package for zero inflated Poisson regression. I would like to calculate (aproximate) prediction intervals for the fitted values. The package itself does not provide them. Can this be calculated analyticaly? Or do I have to use bootstrap? What I tried until now is to use bootstrap to estimate these intervals. Any comments on the code are welcome.

R-users, My question is related to earlier posts about benefits of quadcore over dualcore computers; I am trying to setup a cluster of windows xp computers so that eventually I could make use of 10-20 cpu:s, but for learning how to do this, I am playing around with two laptops. I thought that the package snow would come handy in this situation, but to use snow, I would probably need to install

I have this problem with R 2.6.1 and 2.7.1 (have not tried others) running on Windows XP: Html help has wrong links to help files. The index.html is correctly found from file:///C:/PROGRA~1/R/R-26~1.1/doc/html/index.html as well as the search engine file:///C:/PROGRA~1/R/R-26~1.1/doc/html/search/SearchEngine.html However, when I search for, say, "paste", I will be pointed to

Dear R helpers, running the following code of a glm model of the family poisson, gives predicted values < 0. Why? library(MASS) library(stats) library(mvtnorm) library(pscl) data(bioChemists) poisson_glm <- glm(art ~ fem + mar + kid5 + phd + ment, data = bioChemists, family = poisson) predicted.values = predict(poisson_glm) range(predicted.values) Thank you in advance for any hints.

hello list! I'm sorry, I just stumbled over this strange behaviour (at least I am not able to explain the behaviour, therefore I assume it to be a strange behaviour): attach(water) # I know, this is not recommended names(water[3:10]) [1] "temp" "pH" "DO" "BOD" "COD" "no3" "no2" "po4" for (i in

Hi, I am new to R. My problem is with the ordered logistic model. Here is my question: Generate an order discrete variable using the variable wrwage1 = wages in first full calendar quarter after benefit application in the following way: * wage*1*Ordered *= 1 *if*0 *· wrwage*1 *< *1000 2 *if*1000 *· wrwage*1 *< *2000 3 *if*2000 *· wrwage*1 *< *3000 4 *if*3000 *· wrwage*1 *<

...udia x <- c(1989,1990,1991,1992,1993,1994,1995,1996,1997,1998,1999,2000,2001, 2002,2003,2004,2005,2006,2007,2008,2009) y <- c(987,937,810,749,1087,807,1050,1294,1455,1022,927,403,698,1191,1078, 1176,1125,936,1263,647,868) dat <- data.frame(x,y) f <- function(x,a,b){a*exp(b*x)} fems <- nls(y ~ f(x,a,b), data = dat, start = c(a=100, b=0)) plot(y~x, xlim=c(1989,2010), ylab="Nesting Females", xlab="Year") curve(f(x, a=co[1], b=co[2]), add=T) [[alternative HTML version deleted]]

Is there a way to force certain formula parameters to be nonnegative? What I want to do is to estimate student capacity over time, namely by > capacity ~ Student + Student:Day I add this formula to a glm call and obtain negative learning slope estimates (Student:Day) in some cases. However, I don't want to allow for that. In such a case, glm should solve > capacity ~ Student and

Hello, I have the TCl/Tk command "tkmessageBox(titel="",message="x",icon="question",type="okcancel")" in my R script. Now I want to perform some operation in relation to the user's choice, something like "if (okpressed) xxx else yyy" What values does this command give and how are they used? Thank you, Richard -- Richard M?ller -

Author: petede Repository: /hg/zfs-crypto/gate Revision: ef5d30f0abff39691966b9805b895f3c757b5a9f Log message: 6269165 misleading comments in usr/src/cmd/stat/iostat/iostat.c 6286482 remove the only occurrence of "shit" in OpenSolaris 6231501 Typo in <sys/fem.h> Contributed by Shawn Walker <binarycrusader at gmail.com>. Files: update: usr/src/cmd/stat/iostat/iostat.c

I have a question for users of MCMCglmm that have experience implementing the zero-inflated poisson model. I find that the documentation, and previous questions, do not offer a lot of clear guidance on specifying and interpreting the zipoisson model. In particular, I see a lot of zero-inflated poisson examples that use the at.level(trait, x):variableName syntax. Specifically, the MCMCglmm

Hello! I just joined this mailing list a couple of minutes ago. I'll start to use CentOS for academic purposes. We'll try to build a cluster based in machines with this OS. First of all, I'd like to beg you for patience, because I'm comple- tly new in Linux. My first question is: every tutorial on installation talks about inserting an installation CD, but I'm not sure about

If I combine elements into a list b <- c(22.4, 12.2, 10.9, 8.5, 9.2) my.c <- sample.int(round(2*mean(b)), 5) my.list <- list(b, my.c) the names of the elements seems to get lost in the process: > str(my.list) List of 2 $ : num [1:5] 22.4 12.2 10.9 8.5 9.2 $ : int [1:5] 11 8 6 9 20 If I explicitly name the elements at list-creation, I get what I want: my.list <- list(b=b,

Users report that an extent map''s list is still linked when it''s actually going to be freed from cache. The story is that a) when we''re going to drop an extent map and may split this large one into smaller ems, and if this large one is flagged as EXTENT_FLAG_LOGGING which means that it''s on the list to be logged, then the smaller ems split from it will also be

Hi useRs, In trying to take R to engineering undergraduate students, I have been looking for context that would make R more accessible to the said audience. Though R is primarily a statistical tool, I would want to demonstrate the use of R for certain engineering courses (Design of Machine Elements - gear design, ball bearings, etc.) which would generate interest in it and provide students a way

Hello, I have a simple question but I don't know which method is best to use for my problem. I have the following strings: str1 <- "My_name_is_peter" str2 <- "what_is_your_surname_peter" I would like to apply predefined abbreviations for peter=p and name=n to both strings so that the new strings look like the followings: str1: "My_n_is_p" str2:

Dear all- I am have trouble in using the model="ht" option in function plm from the plm library. I am using Package: plm Version: 1.1-1; R version 2.8.1 (2008-12-22) running on a FC-8 linux machine. Here is what I am trying to do: ##---------------------------------------------------------------------------- R> ###Prob 6 Chapter 3 Use R! Applied Econometrics with R (Kleiber

http://bugs.digium.com/bug_view_page.php?bug_id=0001429 Saying numbers is not always easy, especially if you want one software to be able to do it in many different languages with different syntaxes for how to construct numbers like "one-hundred-twenty-four" or "fem-hundra-tjugo-?tta" or "quatre-vingt-dix-neuf". All the code that we had on the bug tracker for

why are people starting now every minute a new thread about the same problem already reported as reply of the release announce? -------- Original-Nachricht -------- Betreff: Re: [Dovecot] v2.2.11 released Datum: Wed, 12 Feb 2014 10:05:17 +0100 Von: Ralf Hildebrandt <Ralf.Hildebrandt at charite.de> Antwort an: dovecot at dovecot.org An: Dovecot Mailing List <dovecot at dovecot.org> *