similar to: Multiple linear regression

Hi, I am having a problem with a very slow indexing and sub-sectioning of a 3d array: > dim(arr) [1] 245 175 150 For each point in the array, I am trying to calculate the mean of the values in its surrounding: mean( arr[ (i - radius):(i + radius), (j - radius):(j + radius), (k - radius):(k + radius)] ) Putting that code in 3

Using the input below, can I do something more elegant (and more efficient) than the loop also listed below to pad strings to a width of 5? The true matrix is about 300K rows and 31 columns. ####################### #INPUT ####################### > temp DX1 DX2 DX3 1 13761 8125 49178 2 63371 v75 22237 3 51745 77703 93500 4 64081 32826 v72 5 78477 43828 87645 >

Dear useRs, I have written an ecological model, based on the epidemiology SIR model. I've been trying to simulate it in R. However, I can't simulate it properly. Two guesses: my script isn't right; I'm not setting the parameters properly I have uploaded an image to the model here: http://img24.imageshack.us/img24/743/imagemutr.jpg The script I am using is as it follows:

Hello, I am writing a package with a collection of several models. In order to allow users to play interactively with the models (in contrast to hacking lengthy scripts), I want to put all what is needed to run a particular model into a single list object for each model. Then there will be a collection of functions to run the model or to modify parameters, time steps, integration method ...,

I'm looking for an extension of kappa to measure agreement among multiple raters when there can be more than one response per subject. For example, say a group of doctors assign diseases to patients. Each patient will be assigned one to many diseases, and the number of doctors assigning diseases to any one patient will be two to many. Here's an extremely simple example of the type of

Hello, I have 2d data where x coordinate is not given in usual ascending order (x1,...,x1+l; l>0), and instead in another, regular, but not ascending or descending order (for illustration: x1,-x1,x1+dx1,-x1-dx1,....).y is an array which corresponds to the way x is ordered. I have noticed that giving 'plot(x,y,type='l') produces a plot where the points are connected in a completely

I am quite new to Linux and have moved (almost) from a windoze NT4 environment. My present configuration is running SuSE V 8.0 with KDE3.0.5 desktop on two machines, connecting with Samba to an NT4 PC, and an occasional laptop or other PC that connects locally to the network. After a deal of searching, researching, and seeking advice I have decided to use Shorewall as my firewall.

I am having a problem with the package SSPIR. The code below illustrates it. I keep getting the message: "Error in y - f : non-conformable arrays." I tried to tweak the code below in many different ways, for example, substituting rbind for cbind, and sometimes I get a different error message, but I could not find a variation of this code that would work. Any help will be greatly

A few weeks back I had tried to install CentOS on my old AMD 5x6x (equiv to P75 performance) with plans to use it as my firewall and router. Because there were some problems relating to RPM compile architecture (apparently), I posted and Bryan suggested IPCop, even at peril to his life on these lists ;-) according to him. ;-) ;-) Just wanted to tell him thanks, it's worked out well. Had

Hello again. Today I got following questions: 1.Which version of OpenGL is wine written in? 3.1? 3.0? 2.1? 2.0? Lower maybe? 2.Is wine using custom OpenGL extensions? From which vendors (NVIDIA I guess, ATI, S3, INTEL or others as well?)? 3.What is the minimum OpenGL version to run DX2/3/6/7,DX8,DX8.1 and DX9.0 apps? (I guess for OpenGL apps it doesn't really matter right? You are just

Hi, I've made several attempts to boot syslinux from the harddisk. I have a win 98 startdiskette with syslinux and copybs. I strictly follow the procedure in the syslinux documentation. Yet, the result of the boot is not that linux is starting. Instead, the Microsoft windows screen flashes for a couple of seconds and is then replaced by a Dos prompt. Syslinux is V2.04. The PC is an IBM Aptiva

Full_Name: Albrecht Gebhardt Version: 0.99.0 OS: alpha, osf4.0 Submission from: (NULL) (143.205.180.40) DEC cc doesn't like c++ style comments // like this one Please use standard cc comments or #if 0 ... #endif A quite large patch follows (I hope it will pass mailing with not too much wrapped lines): --- ./src/main/g_her_glyph.c.dec-cc.patch Mon Feb 7 14:48:10 2000 +++

Hi 'tinc users' I'm new whit tinc and this mail-list. I like to build a vpn over the internet between two sites in Holland and Germany. Are there some experience whit the performance (turn-a-round en bandwidth) of a link. I want to use it for a 1,5Mb connection. Greetings, Fred Krom. - Tinc: Discussion list about the tinc VPN daemon Archive:

As this forum proved to be very helpful, I got another question... I'd like to fit data points on which I have an error, dx and dy, on each x and y. What would be the common procedure to fit this data by a linear model taking into account uncertainty on each point? Would weighting each point by 1/sqrt(dx2+dy2) (and taking dx and dy as relative errors) in a lm() fit do the job? I would like to

The nbdkit tests take a long time to run especially on armv7. To try to get an idea of which tests take the most time, and since I happen to have a qemu aarch64 *TCG* guest available, I ran the tests on that guest. A TCG guest is about the worst case for libguestfs, used by some tests, since it uses double-nested TCG, and as you can see below there are tests taking over 6 hours. In fact I

Hello, Are there any packages/functions available for testing the assumptions underlying assumptions for a good GLM model? Like linktest in STATA and smilar. If not, could somebody please describe their work process when they check the validity of a logit/probit model? Regards, Jay

disable bogus touchpad device model, which cause annoying dmesg on 2.6 kernel Signed-off-by: Edwin Zhai <edwin.zhai@intel.com> diff -r 7c951e3eb5ab -r 14d3890843d0 tools/ioemu/hw/pckbd.c --- a/tools/ioemu/hw/pckbd.c Wed Oct 19 11:53:00 2005 +0100 +++ b/tools/ioemu/hw/pckbd.c Thu Oct 20 18:28:18 2005 +0800 @@ -29,6 +29,9 @@ /* debug PC keyboard : only mouse */ //#define DEBUG_MOUSE +/*

hi there here is part of SPECS file: ------------------------------------------------------------------ CFLAGS="$RPM_OPT_FLAGS" ./configure \ --prefix=%{_prefix} \ --localstatedir=/var \ --sysconfdir=/etc \ --with-privatedir=%{_sysconfdir}/samba \ --with-fhs \ --with-quotas \ --with-smbmount \ --with-pam \

Hi, as some of you may know, there were some issues reported recently about the Steam version of this game not starting. Presumably the solution was to compile wine 1.2+ with a recent version of GCC and binutils, so that the copy protection software could hook certain functions. I've done this, and I've posted some details in the bug tracker ( http://bugs.winehq.org/show_bug.cgi?id=23553).

Dear R-users, I need to produce a histogram where for every breaks there are the mean of the data. I tried tu use the function >hist(x, break=20 ... ) but this return the numerosity for every breaks, not the mean. Any hint? Thanks in advance, francesco --------------------------------------- Francesco Nutini CNR-IREA Ist. per il Rilevamento Elettromagnetico dell'Ambiente Via