search for: fem

...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) fem.frek<- c(0.054, 0.099, 0.04, 0.113, 0.234, 0.241, 0.099, 0.182, 0.176, 0.175, 0.126, 0.09) rel.femfec<- c(0.975, 0.713, 1.035, 1.745, 0.77, 1.162, 0.94, 0.636, 2.12, 2.159, 0.428, 1.217) emil.data<- as.data.frame(cbind(pop.size, fem.frek,asfem.frek,rel.femfec)) attach(emil.data) contour(po...

...nalyticaly? 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. The data and the model are based on the examples in zeroinfl(). #aproximate prediction intervals with Poisson regression fm_pois <- glm(art ~ fem, data = bioChemists, family = poisson) newdata <- na.omit(unique(bioChemists[, "fem", drop = FALSE])) prediction <- predict(fm_pois, newdata = newdata, se.fit = TRUE) ci <- data.frame(exp(prediction$fit + matrix(prediction$se.fit, ncol = 1) %*% c(-1.96, 1.96))) newdata$fit <-...

...up MPI for a single windows machine, but I am not sure whether this is also intended to mean that one can not expect it to work with multiple computers. Regards, Samu M?ntyniemi ------------------------------------------ Samu M?ntyniemi Researcher Fisheries and Environmental Management Group (FEM) Department of Biological and Environmental Sciences Biocenter 3, room 4414 Viikinkaari 1 P.O. Box 65 FIN-00014 University of Helsinki Phone: +358 9 191 58710 Fax: +358 9 191 58257 email: samu.mantyniemi helsinki.fi personal webpage: http://www.helsinki.fi/people/samu.mantyniemi/ FEM webpage: ht...

..., therefore I suspect that the problem would be somehow related to the settings of this particular laptop. However, I do not know where to look. Cheers, Samu -- ------------------------------------------ Samu M?ntyniemi University Researcher, PhD Fisheries and Environmental Management Group (FEM) Department of Biological and Environmental Sciences Biocenter 3, room 4414 Viikinkaari 1 P.O. Box 65 FIN-00014 University of Helsinki Phone: +358 9 191 58710 Fax: +358 9 191 58257 email: samu.mantyniemi at helsinki.fi personal webpage: http://www.helsinki.fi/people/samu.mantyniemi/ FEM webpage:...

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. Best regards, P. Olsson [[alternative HTML version deleted]]

...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 names(water)[3:10]){ fname<-paste("Henni/GFX/fem",i,".png",sep="") mname<-paste("Henni/GFX/mal",i,".png",sep="") png(fname,1000,1000) xyplot(N_female~eval(parse(text=i)) |group,xlab=i,ylab="Abundance") graphics.off() png(mname,1000,1000) xyplot(N_mal...

...< *3000 4 *if*3000 *· wrwage*1 *< *4000 5 *ifwrwage*1 *¸ *4000. Estimate an order logit model using the new created discrete variable as the dependent variable with the following control variables: black = 1 if black hispanic = 1 if hispanic otherrace = 1 if black = 0 and hispanic = 0 female = 1 if female wba = claimant's UI weekly bene¯t amount (*=week*) age = age of claimant when claim is ¯led (years) bp1000 = base period wages (thousand $) inuidur1 = duration of unemployment (weeks) I have written the program as: > pennw <- read.table("G:/Econometrics II/penn...

...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 update: usr/src/cmd/volmgt/vold/dev_rmscsi.c update: usr/src/uts/common/sys/fem.h

...In particular, I see a lot of zero-inflated poisson examples that use the at.level(trait, x):variableName syntax. Specifically, the MCMCglmm course notes, available on the package's CRAN page, uses the following example on page 102: m5d.1 <- MCMCglmm(art ~ trait - 1 + at.level(trait, 1):fem + at.level(trait, 1):mar + at.level(trait, 1):kid5 + at.level(trait, 1):phd + at.level(trait, 1):ment, rcov = ~idh(trait):units, data = bioChemists, prior = prior.m5d.1, family = "zipoisson", verbose = FALSE) I have been unable to find an answer to the following questions and would ap...

...dows, trying to install CentOS 4 in a 32 bit's PC. I've downloaded CentOS-4.5-i386-bin.......iso. Thank you for any help! ~ J. Labaki Computational Structural Mechanics Labs Dept. of Computational Mechanics Faculty of Mechanical Engineering Campinas State University/SP - Brazil http://www.fem.unicamp.br/~labaki "Time is the best teacher; unfortunately, it kills all its students."

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,

...with EXTENT_FLAG_LOGGING, because the log code holds one reference. The end result is the warning, but the truth is that we set the flag EXTENT_FLAG_LOGGING only during fsync. So clear flag EXTENT_FLAG_LOGGING for extent maps split from a large one. Reported-by: Johannes Hirte <johannes.hirte@fem.tu-ilmenau.de> Reported-by: Darrick J. Wong <darrick.wong@oracle.com> Signed-off-by: Liu Bo <bo.li.liu@oracle.com> --- fs/btrfs/file.c | 1 + 1 files changed, 1 insertions(+), 0 deletions(-) diff --git a/fs/btrfs/file.c b/fs/btrfs/file.c index 83c790d..7bdb47f 100644 --- a/fs/bt...

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:

...reign") R> fulldat <- read.dta("~/Desktop/psidextract.dta") R> library("plm") R> R> fulldat.plm = plm.data(fulldat,index=c("id","t")) R> R> earn_plm <- plm(lwage~ occ+ south+ smsa+ ind+ exp+ exp2+ wks+ + ms+ union+ fem+ blk+ ed | exp+ exp2+ wks+ ms+ union+ ed, + data = fulldat.plm,model="ht") Error in names(result) <- nf : attempt to set an attribute on NULL I have tried several variations and some other data sets (not so easily reproducible for others as this one) but have yet to o...

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 various syntaxes will now be included into one combined patch and moved forward to be integrated into the Asterisk CVS. If you have code for another language, this is the right time to co...

...xed assert-crash when syncing mailbox deletion. pop3 keeps crashing with a sigsegv in that version. Gotta go back. Feb 12 10:02:52 postamt dovecot: pop3(mhippler) LJZ31jHyygAKK4EU: Fatal: master: service(pop3): child 14635 killed with signal 11 (core dumped) Feb 12 10:02:54 postamt dovecot: pop3(fem-qm-sfb) HlCP1jHyGgCNKtkP: Fatal: master: service(pop3): child 14692 killed with signal 11 (core dumped) -------------- next part -------------- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 246 bytes Desc: OpenPGP digital signature URL: <http://...