similar to: Generating Data using Formulas

Hi, after looking around I was not able to get info about these two questions: 1. Is there a function to have a "jackknifed corrected var/cov estimate" (as described in MacKinnon and White 1985) in a standard OLS regression? 2. Does R possess a LAD (Least Absolute Deviation) regression function? Any help? Thanks -- ======================================================== Angelo

Does anyone have an R function for the estimation of Tobit regression models (for censored responses), as described, e.g., in Davidson, R. and MacKinnon, J.G., 1993, Section 15.7, Estimation and Inference in Econometrics, Oxford University Press? Thanks. FC. -- Francisco Cribari-Neto voice: +55-81-32718420 Departamento de Estatistica fax: +55-81-32747425 Universidade

Full_Name: James MacKinnon Version: Version 1.0.0 OS: Debian 2.1 Submission from: (NULL) (24.141.113.138) I am a novice at R, so I am not totally sure that this is a bug, but I think it must be. When I type help.start() or ?help, Netscape starts up, as it should, but it displays my normal home page, not any help for R. The relevant documents exist in /usr/lib/R/doc/html, but Netscape does not

Has any one run sage line 50 with samba As i have setup a samba server to hold the data directory for sage was fine for bout week then slowed right down a samba restart seemed to help but sage still runs rather slow We are only talking about 5 users at a time the data dir is about 2gb in size Any help would be handy or tweaks for samba Samba version 2.2.8a Freebsd 4.9 2.4 cpu with 512 ram 100

Dear R-community, I am just conducting various mediation analyses with data from individual participants who are nested in teams. I found the ?sobel.lme? function to be pretty helpful when doing so. (A big thank you to Paul Bliese for his multilevel-package and the awesome manual!) Here?s the example from the Help-file: data(bh1996) library(nlme) #A small but significant indirect effect indicates

Is there an option for running SUR and 2SLS regressions with weighting (I am analysing mortality in towns, hence want to weight by population size) Many thanks Jon Anson -- Yonathan (Jon) Anson Department of Social Work Ben Gurion University of the Negev 84105 Be'er Sheva, Israel. Tel: +972 8 647 93 14(w) +972 8 6489286 (h) 067 233279 (m) Fax: +972 8 647 29 33

Okay, I've read as much as I can, and I think i've followed instructions, but I'm still having problems with * and firefly... I can get outgoing to other freshtel working, but not incoming (I get the "not available" voicemail), or outgoing to landline. I'm using the debian asterisk package (0.9.1-RC1-4) My iax.conf has in general (under my FWD register, which

Hi All, Is the Levenberg-Marquardt algorithm available in R. This method combines the steepest descent algorithm and Newton's method. Thanks in Advance, Dermot MacSweeney. ************************************************************** Dermot MacSweeney NMRC, Email: dsweeney at nmrc.ucc.ie Lee Maltings, Tel: +353 21 904178 Prospect Row, Fax: +353 21 270271 Cork, WWW:

Full_Name: Petr Savicky Version: 2.4.0 OS: Fedora Core release 2 Submission from: (NULL) (62.24.91.47) the error is > binom.test(0.56*10000,10000) Error in binom.test(0.56 * 10000, 10000) : 'x' must be nonnegative and integer while > binom.test(5600,10000) yields correct result. The same error occurrs for > binom.test(0.57*10000,10000)

Hi guys. Quite happy with mocha for testing. It''s been a bit of an eye opener for me, and I expect to begin using it heavily. One other thing I''d like to do is stub out particular class methods in development. Is there a way I can do this with mocha? Jason

Consider the following little "benchmark" > require(Matrix) > tmp <- Matrix(c(rep(1,1000),rep(0,9000)),ncol=1) > ind <- sample(1:10000,10000) > system.time(tmp[ind,]) user system elapsed 0.004 0.001 0.005 > ind <- sample(1:1000,10000,replace=TRUE) > system.time(tmp[ind,]) user system elapsed 0.654 0.006 0.703 >

I would like to compare every column in my matrix with every other column and get the r-squared. I tried using the following formula and looping through every column: > summary(lm(matrix[,x]~matrix[,y]))$r.squared If I have 10,000 columns, the loops (10,000 * 10,000) take forever even if there is no formula inside. Then, I attempted to vectorize my code: > cor(matrix)^2 With 10,000

Dear readers we are using Dovecot 2.2.7 and all of our users are using Thunderbird as their mail client. Some of them additionally use their iPad/iPhone and a very few an Android Mail-Client. Now one user noticed that two of his mail folders disappeared. He first believed that he accidentally deleted those folders but then he realized that they are still visible from his iPad. I checked this

ogg123 seems to have a issue when streaming for long peroids of time. If I have ogg123 up connected to a streaming server for more than 3.5 hours it begins to stutter or add blocks of static. I've verifed this with 2 different soundcards in both alsa oss emulation and regular oss. It's 100% repeatable for me. It does not happen with xmms using libvorbis as the output plugin so it

Hi ! When trying to change users password from NT or XP workstations (using CTRl-ALT-DEL), I get an error under Windows and this in my samba logs: Sep 10 16:15:08 testbox kernel: pid 2709 (smbd), uid 10000: exited on signal 6 Sep 10 16:15:08 testbox smbd[2709]: [2003/09/10 16:15:08, 0] lib/util_sec.c:assert_gid(114) Sep 10 16:15:08 testbox smbd[2709]: Failed to set gid privileges to

Hi there Trying to solve an issue with samba and windows 7 permissions dialog. Problem is that sometimes windows 7 permissions dialog is lacking ldap users and groups. Looks like my problem is related to this one: https://forums.freenas.org/index.php?threads/users-and-groups-not-showing-up-in-windows-7.46023/ Sadly there is no solution in that thread. Consider the following setup: linux debian

## This is an example where ``-quoted non-syntactic names work in ## simple formulas, but not in formulas with an Error() term. ## Is this intentional or an oversight when ``-quoted names were added? tmp <- data.frame(y=1:6, `x^1`=factor(c(1,2,1,2,1,2)), `x^2`=factor(c(1,1,2,2,3,3)), check.names=FALSE) tmp aov(y ~ `x^1` + `x^2`, data=tmp)

Dear useRs, version 1.0-0 of the "Formula" package has just been released on CRAN http://CRAN.R-project.org/package=Formula accompanied by an article in the Journal of Statistical Software http://www.jstatsoft.org/v34/i01/ It provides simple infrastructre for processing formulas like y ~ x1 + x2 | z1 + z2 where there are two groups of explanatory variables, e.g.,

Dear useRs, version 1.0-0 of the "Formula" package has just been released on CRAN http://CRAN.R-project.org/package=Formula accompanied by an article in the Journal of Statistical Software http://www.jstatsoft.org/v34/i01/ It provides simple infrastructre for processing formulas like y ~ x1 + x2 | z1 + z2 where there are two groups of explanatory variables, e.g.,

> From: Peter Dalgaard > > Hallgeir.Grinde at elkem.no writes: > > > Dette er en melding med flere deler i MIME-format. > > --=_alternative 004613C000257091_= > > Content-Type: text/plain; charset="US-ASCII" > > > > And some more informastion I forgot. > > R does not crash if I write out the formula: > > > > set.seed(123)