similar to: problems linking c++ code

Hello, I am trying to do some simulation using the rmetasim package and I've run to this problem. --beginning of error msg-- Error in "[<-"(`*tmp*`, slice[l, ], slice[l, ], value = c(0.200000002980232, : number of items to replace is not a multiple of replacement length --end of error msg-- Here is the script I used. --script starts here-- ## load 'rmetasim'

Dear R user: Now I am studying "matrix population models" and Rmetasim package,but I find it's very difficult to understand this model fully and can't find a good teacher.My question is that:Now I do some research about one marine shrimp: Penaeus chinensis, it's life cycle is one year. after it breed it's offsprings, it will died. This shrimp's life cycle is dicrete

Dear listers, I need to export a distance matrix in the following format : a 0.0 1.0 0.2 1.0 1.0 b 1.0 0.0 1.0 1.0 1.0 c 0.2 1.0 0.0 1.0 1.0 d 1.0 1.0 1.0 0.0 1.0 e 1.0 1.0 1.0 1.0 0.0 I tried write.matrix() from the MASS library, which gives: a b c d e 0.0 1.0 0.2 1.0 1.0 1.0 0.0 1.0 1.0 1.0 0.2 1.0 0.0 1.0 1.0 1.0 1.0 1.0 0.0 1.0 1.0 1.0 1.0 1.0 0.0 Does anyone know of an R trick to change

Whether the two variables have the same units does not matter. Moreover, even if there were some way of converting cm to kg the correlation would still be the same because the correlation is invariant under unit conversion as it is invariant under multiplication of its arguments by a constant. As for your second question, the correlation estimator is a continuous function of each of the

Dear Listers, I am trying to match tables that DO NOT have the same length. The tables result from the function "table()" so they look like this: table 1 2 3 4 3 5 7 table 2 1 2 3 6 4 5 I need the following output: (NOTICE THE ZEROS) 1 2 3 4 table1 0 3 5 7 table2 6 4 5 0 Unfortunately, I was not successful using "match()". Previous postings explain how to do

Sorry for the extra info required, but I had digest mode turned on and wouldn't have received the mailing till tomorrow. I have since updated the data manually, but after 'yum update tzdata' was run, the zdump -v for EST5EDT and America/New_York all still showed a date of April 1, instead of March 11. After manually fixing, it is correct. This was mostly a question about why yum

I realized that the script I included in my original post (see below) was written for an older version of rmetasim and will not work with the current version. The only change that need be made to get my script to run with the latest release of rmetasim is to change the command 'simulate.landscape()' to 'sim.landscape()'. My apologies to anyone who tried to use my script and

Dear listers, I am trying to assess variance components for a nested, mixed-effects model. I think I got an answer that make sense from R, but I have a warning message and I wanted to check that what I am looking at is actually what I need: my data are organized as transects within stations, stations within habitats, habitats within lagoons. lagoons: random, habitats: fixed the question is:

Dear Listers, I emailed the list a few days ago about how to bootstrap a community matrix (species by sites) and get a consensus tree with node support. A friend pointed out that a similar question remained unanswered in 2004. I wish to re-word my question: is anyone aware of a package / method to obtain a majority-rule consensus tree from x distance matrices ? Is anyone using R to generate

Hello all, R just ''stopped working'' on my machine (ibook g4, os 10.3.9). between last time I used it and now, I did not do anything but check emails and browse the internet. when I try to start the GUI, it appears and disappears as rapidly. when I tried to run a session from the terminal, I got the following: mkdir: /tmp/Rtmp590-29493: Permission denied mkdir:

I'm trying to use a contributed package (rmetasim) to generate simulated genetic datasets. My scripts work fine when I run them on a Sun workstation running Solaris 7 and when I run them on a ~4 year old laptop PC that I have. However, when I run them on my new laptop (Dell Latitude D410 purchased in August 2005), the simulation will run for a short (variable) period of time, then R

Hi All, I'm trying to specify the formula for a glm based upon the names in the data.frame. I always discard the first two variables in the frame and the third will always be the only response variable. The problem appears because I would like to specify the rest of the variables of the data.frame as independent variables in the formula. I can't think of a way to do this cleanly and

I tested moving a box from Tao 1.0 to CentOS 3.7. All went well, as far as I can tell, with the OS upgrades, but I needed to rpm -ivh --force the desktop-* stuff. This was a minor thing to fix for no real purpose other than to just be doing it. I still get, though, the Tao screens when the graphic login screen appears. I have changed the Application->Settings->Login Screen->

I apologize for this very naïve and stupid question. I have searched the archives and "googled" a number of different queries relating to this question and cannot find the answer. My Division is moving toward running an application server for all statistical applications because the Department of Veterans' Affairs has their computers locked down. They are locked down so much so

Question, I have made numerous attempts at configuring the samba 2.3.3a program to include support for MIT Kerberos 5 v.1.2.3. The configure -help says that the switch to turn it on is"--with-krb5=base-dir". I have tried several different directories as the base-dir with no luck. During the coufigure, completes but the with-krb5 options is set to "no". The only time I

Hi all, I am interested in calculating and displaying the distributions of pairwise comparisons between DNA sequences in populations. The comparisons are the number of nucleotide sites that differ between the two sequences (mismatches). My sequences are stored in a vector of strings. There is an additional vector of the same length that provides the indices to the DNA sequences. Finally, I

Hello I am having an R related problem and can't seem to fix it by looking at any of the online documentation. I have a laptop (using Windows XP home edition) and loaded R 2.2.1. I then designated a mirror (USA-NC). But when I try to load a package, it freezes. I have tried to close R, but it won't close; I end up having to restart the laptop. I have no problems with any other

Hi all, This is as much as statistical/estimation question as an R-specific one, but here goes. I am trying to use logistic regression to predict suitability of habitats for certain plant species. The response variable is a binary one that indicates whether a particular species is found at a site on the landscape. The independent variables represent physical characteristics of the landscape

Hi, I have worked long and hard and looked in the manuals and am having a hard time constructing a diagonal matrix. I can get the diagonals out of a matrix but can't construct the matrix with just the diagonals. I have been on the web site and manuals and I think that it says to use: dsj <- diag (three = 1, nrow, ncol = 7) three is the name of my matrix and

I'm running CentOS 3 and have just installed a Liteon lightscribe DVD/CD unit. I firstly realized I didn't know exactly how to mount it, as /dev/dvd didn't exist, I didn't want to generate anything if it already existed, and I didn't want to mess up my cdrom unit that works just fine. I mounted the device though the X window just to see how it might mount, and find that it