similar to: circular filter

Hi, I have a dataframe containing data from individuals 1, ..., 12 (grouping variable "g" in the data frame below), which belong either to "A" or "B" (grouping variable "f"): set.seed(1) tmp <- data.frame(

Hi R-helpers I have a dataframe DF (lets say with the variables, y, x1, x2, x3, ..., clust) containing relatively many NAs. When I fit an ordinal regression model with the function lrm from the Design library: model.lrm <- lrm(y ~ x1 + x2, data=DF, x=TRUE, y=TRUE) it will by default delete missing values in the variables y, x1, x2. Based on model.lrm, I want to apply the robust covariance

I have a spatial data set with ordinal response variable containing four levels. I would like to know if and how spatial autocorrelation can be taken into account when ordinal logistic regression is used (e.g. the function lrm from the Design package). Thanks for your help! Christof

Dear R-users, has anybody combined the glm function with the Newey-West estimator of variance, similar as in Stata 7.0? I'd like to estimate corrected standard errors within a logistic regression model, taking into account the auto-correlated binary observations within individuals. I use R1.5.1 on Mac OS X (10.2). Thanks, Christof

Dear R users, is there a way to extract directly rows instead of column vectors using the function sample()? I have already tried to transpose and back-transpose my dataframe, but I couldn't manage to get rid off the levels in the resulting dataframe. Macintosh G4 (Mac OS9.1), R Version 1.3.0 Any help is welcome! Christof -- Christof Bigler Mountain Forest Ecology ETH-Zentrum HG G21.4

Is there any way to fill the area within a plotted line and a horizontal line, e.g. in the following example: x <- rnorm(100) plot(x,type="l") abline(h=0) Thanks for any help. Christof

I want to assign different colors to different strings using mtext(...). However, when I use something like mtext(text=c("string1","string2","string3"), col=c("black","blue","red"), side=3, line=0) string2 and string3 are printed over string1. When I use paste("string1","string2","string3"), the series

Dear R users, can anybody explain me why the function glmmPQL(.) behaves in different ways, depending on the number of measurements/individuals you use? To show you this, I generated two examples. The first one includes 20 indivduals with each 100 repeated measurements (binary response), the second one includes 40 individuals. The 'individuals' differ only in different x values. I

Le 25/08/2019 ? 23:07, Stephen John Smoogen a ?crit?: > I have had several people explain to me that my comments were > out-of-place and out-of-turn. When I wrote them originally I thought I > was going from experience of previous releases and thought I was being > matter of fact. In rereading I realized I was crass and grumpy. My > un-reserved apologies for that. I should have

Dear R-users, Is there a procedure to identify neighbors in a regular lattice using either a "rook" or a "queen" criterium? To be more specific, suppose: My lattice: 1 2 3 4 5 6 7 8 9 "Rook" Neighbors 1 has neighbors 2,4 2 has neighbors 1,3,5 ... 5 has neighbors 2,4,6,8 "Queen Neighbors" 1 has neighbors 2,4,5 2 has neighbors 1,3,4,5,6 For each case,

raf wrote: > Warlich, Christof wrote: > > ... > > I want to be able to ssh to all internal hosts that live in the internal.sub.domain.net, > > i.e. that are only accessible through the internal.sub.domain.net jumphost without > > having to list each of these hosts somewhere, as they may frequently be added or > > removed from the internal domain and without being

Dear R experts, I am plotting the population of students who live in a city, and in successive circular bands made of the contiguous districts that surround the city. This is a stylized figure, where I specify the area of each successive circle based on the cumulative population of students. I want to compare two sets of concentric circles across different populations - such as 'All

When the world was young, "M. Tyson Bigler" <bigler@shellus.com> carved some runes like this: > I am having difficult getting the hosts allow/deny to work with IP ranges > or netmasks. Ideally I would prefer to use the IP/Netmask method, e.g. > 111.22.33.0/255.255.255.0, because we have a flat ATM network with > multiple class C's (netmask is really

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi, As a result of not being able to control the shares, (invisible network), i wanted to enable swat, but the files are not were they supposed to be.... Can anyone help me? - -- Have a nice day, M9. Now, is the only time that exists. OS: Linux 2.6.18.2-34-default x86_64 Huidige gebruiker: monkey9@tribal-sfn2 Systeem:

I have had a RH5.2 server running Samba up for a long time. I upgraded to Mandrake 6.0. It broke Samba, so I reinstalled using the RPM on the CD. Samba works just fine, but SWAT will not! If I try to start it manually (/usr/sbin/swat) it just hangs. In /etc/services I have a line like this: swat 901/tcp In /etc/inetd.conf I have: swat stream tcp nowait.400 root

Hello,everybody. My name is Chuan Zun Liang. I come from Malaysia. I am just a beginner for R. Kindly to ask favor about median filter. The problem I facing as below: > x<-matrix(sample(1:30,25),5,5) > x [,1] [,2] [,3] [,4] [,5] [1,] 7 8 30 29 13 [2,] 4 6 12 5 9 [3,] 25 3 22 14 24 [4,] 2 15 26 23 19 [5,] 28 18 10 11 20 This

Hello, I have a little question and this regarding the compilation of the latest release of OpenSSH on a Linux Slackware version 8 box. We are currently using Kerberos 5 for user authentification and I saw that in SSH there is only an option to configure called: --with-kerberos4, so my question is: what do I need to do to get Kerberos 5 support into OpenSSH ? I am using the MIT kerberos version

Hi I want to do something like a <- c(10,20,15,43,76,41,25,46) tab <- data.frame(a) name <- "a" for (v in unique(tab[[name]])) { r <- subset(tab, name==v) # this does not work ... } i.e. a "string" on the left side of the select expression (subset). How could I solve this? thx Christof

Hi xyplot(y ~ x | subject) plots a separate graph of y against x for each level of subject. But I would like to have an own function for each level. Something like xyplot(y ~ x | subject, panel = function(x,y) { panel.xyplot(x,y) panel.curve(x,y) { # something that dependents on the current subject ... } }) How I get the current

require(circular) c <- circular(rep(0, 20), zero = pi/2, rotation = 'clock') plot(c, stack = TRUE, shrink = 1.5) Can anyone tell me why the stack is offset from 0? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick