similar to: transition matrix problem

Hello R-help Readers I'm currently in the process of trying to write a permutation test procedure for looking at differences between groups with a multivariate data set (something equivalent to ANOSIM - analaysis of similarities for those familiar with this test). As with other permutation tests, for cases where there are large numbers of possible permutations, randomly sampling from these

Dear all, I've taken a subset of data from a data frame using crb<-subset(all.raw, creek %in% c("CR") & year %in% c(2000,2001) & substrate %in% ("b")) this works fine, except that all of the original factor levels are maintained. This results in NA's for these empty levels when I try to do summaries based on factors using by(). Is there a simple way to drop

Hi I have a 3 x 2 contingency table: 10 20 30 40 50 60 I want to update the frequencies to new marginal totals: 100 130 40 80 110 I want to use the ipf (iterative proportional fitting) function which is apparently in the cat package. Can somebody please advice me how to input this data and invoke ipf in R to obtain an updated contingency table? Thanks. By the way I am quite new to R. -- Dr

How do I remove a "$" character from a string sub() and gsub() with "$" or "\$" as pattern do not work. > sub("$","","ABC$DEF") [1] "ABC$DEF" > sub("\$","","ABC$DEF") Error: '\$' is an unrecognized escape in character string starting "\$" >

I have some square matrices with na values in corresponding rows and columns. M<-matrix(1:2,10,10) M[6,1:2]<-NA M[10,9]<-NA M<-as.matrix(as.dist(M)) print (M) 1 2 3 4 5 6 7 8 9 10 1 0 2 1 2 1 NA 1 2 1 2 2 2 0 1 2 1 NA 1 2 1 2 3 1 1 0 2 1 2 1 2 1 2 4 2 2 2 0 1 2 1 2 1 2 5 1 1 1 1 0 2 1 2 1 2 6 NA NA 2 2 2 0 1 2 1 2 7 1 1 1 1 1 1 0 2 1 2 8

Hi all, I'm hoping someone might help with a query about conditionally applying formulas to a dataframe. In essence I have 3 lookup tables (Table A, B & C) and a dataframe with a variable Type.Code, which identifies the Lookup Table to which each record belongs. The lookup tables reference different sensor types for which I need apply a different formula to values in Column3 in each row

Can I get some advice on what is required to claim a GPL for a package? A package I maintain (pcurve) uses some S code from statlib. I have an email from the author of the S code saying that he would be happy to issue a GPL for the code. Is the email from the author enough, or do I need something more formal? Chris -- Chris Walsh CRC for Freshwater Ecology Water Studies Centre PO Box 23

We have used the Syntax TotalaNET product (TAS) for some years. This provides similar funtionality to Samba. However the current version (5.2) does not support Win98 and the upgrade cost is very high. In looking at Samba as an alternative, I have two questions: 1. Does Samba fully support Win98? 2. Our users move easily between PCs, X-terminals and workstations, depending on their

I am trying to get a dll compiled for use with dyn.load. I use R.1.7.1 under Windows. I have tried the following trivial example based on the "Writing R extensions" manual. rtest.h -------- class X { public: X (); ~X (); void Give7(double*); }; class Y { public: Y (); ~Y (); }; rtest.cpp --------- #include <iostream.h> #include "rtest.h" static Y y;

Hi all, I have two samba machines running v1.9.18p4. When I try and list the shares on the second machine (calvin) fromn the first machine (hobbes) I get this: Session request failed (131,129) with myname=HOBBES destname=CALVIN Not listening for calling name Try to connect as another name (instead of HOBBES) You may find the -n option useful for this Any idea what this is trying to tell me?

Hi, everyone I've been trying to make a scatter plot of POSIXct format data to see the difference between two start time. Like X: [1] "1995-05-26 19:00:00" "1995-05-27 01:00:00" "1995-05-27 07:00:00" [4] "1995-06-29 01:00:00" "1995-06-29 07:00:00" "1995-06-30 13:00:00" [7] "1995-06-30 19:00:00" "1995-07-01

I am getting the following erro rmessage in ordistep. I have a number of similarly structured datasets using ordistep in a loop, and the message only occurs for some of the datasets. I cannot include a reproducible sample - the specific datasets where this is occur ing are fairly large and there are several pcnm's in the rhs of the formula. thanks for any pointers that may allow me to

Hello, I am the author of the Bonk audio compression program... i've just been looking at your comparrison table, and i noticed bonk gets marginally better compression than Flac on some files (actually i was rather surprised to see bonk on the list at all, it's not exactly high profile :-) ). Bonk in lossless mode is a pretty naive implementation of a predictive coder, so the main

Hi, everyone I got some problems when trying to subscript the value of vector and list, by using calculated indices. Here is the vector I am generated lon<-rep(0,886);lat<-rep(0,691) for (i in 1:886){ lon[i]<-112+0.05*(i-1) } for (i in 691:1){ lat[i]<--44.5+0.05*(691-i) } For a given location of lon(xp) and lat(yp), I would like to calculate the position of them,

I'm trying to fit a mixed-effects model of the form Y = a + bX + cZ + e where X are fixed effects and Z are random. i.e., c is a vector of random coefficients with mean 0. There is no "grouping" variable such as would be used with longitudinal data. Can the nlme package handle this sort of thing? Although it is a simpler model than would be used with grouped data, I can't seem

Hi all: I am trying to use the BiodiversityR package in Rcommander since I need to box-cox transform some data. However, even when trying with the example files, I get this error just after selecting the box cox transformation in the Menu, and Rcommander and Rstudio just shut down. What is it? What could I do? Thanks a lot The instruction at "0x00000001" referenced memory at

The generic stats::median method is defined as median <- function (x, na.rm = FALSE) {UseMethod("median")} I suggest that this should become median <- function (x, na.rm = FALSE, ...) {UseMethod("median")} This would allow additional S3 methods to be developed with additional arguments. Currently I have to over-ride this generic definition in the demography

I am receiving the above error ( full r session output below) the script runs OK in windows. and "genotypes" and "ploidy" are both correct arguments any suggestions would be most welcome Nevil Amos MERG/ACB Monash University School of Biological Sciences > library(Geneland) Loading required package: RandomFields Loading required package: fields Loading required

>>>>> Martin Maechler <maechler at stat.math.ethz.ch> >>>>> on Mon, 27 Feb 2017 10:42:19 +0100 writes: >>>>> Rob J Hyndman <Rob.Hyndman at monash.edu> >>>>> on Wed, 15 Feb 2017 21:48:56 +1100 writes: >> The generic stats::median method is defined as median <- >> function (x, na.rm = FALSE)

Hi, can you advise me is there any ROC(Receiver Operating Characteristic)analysis program in R? Thanks, Dechao ===== Dechao Wang Tel: (44) 01223 719718 Mob: (44) 07729 411134 __________________________________________________ Everything you'll ever need on one web page from News and Sport to Email and Music Charts