similar to: OT: Software for specific visualisation of data...ideas?

Hello there, I'm trying to plot vectors with p<0.1 in a NMDS ordination plot using p.max. Below the scripts I'm using. I guess I'm missing something! could you please give me a hand? species<-metaMDS(species_matrix)ef<-envfit(species,environmentaldata_file,permu=999,na.rm=TRUE)efplot(species, dis="sites")plot(ef,p.max=0.1) Error in plot.envfit(ef, p.max = 0.1) :

Hi I am trying to run Factor Analysis using R...I am using the syntax factanal(m1, factors=3) but it's giving me an message Error in cov.wt(z) : 'x' must contain finite values only ...I am using a data set which is having only numeric variables and have some NA's also in it..What should I do next..Someone please help me out with the syntax..Thanks in advance Cheers Arup --

Sorry to disturb you, but can you help me to solve one little problem? I want to draw a graphic after another with R but I cannot find the first one after that. Do you know the command to hold the graphic with R? I remember with Matlab you may use "hold on". Thanks. -- Baohua Yang Email:yangbaohua@gmail.com [[alternative HTML version deleted]]

Dear R, the condition: identical(length(x),1) returns FALSE but print(length(x)) returns 1 and: is.vector(x) is TRUE. is.integer(length(x)) is TRUE length(x) ==1 is TRUE I am puzzled. Regards -- Corrado Topi Global Climate Change & Biodiversity Indicators Area 18,Department of Biology University of York, York, YO10 5YW, UK Phone: + 44 (0) 1904 328645, E-mail: ct529 at york.ac.uk

Hi, there, I cannot get accurate value for calculation. for example: ld<-sqrt(1*0.05*0.95*0.05*0.95) 0.05*0.95-ld=-6.938894e-18 0.05*0.95-ld==0 is False. I met this problem in my program, how can I handle it. Thanks. xj.

Dear List, I am writing a formula method for a function in a package I maintain. I want the method to return a data.frame that potentially only contains some of the variables in 'data', as specified by the formula. The problem I am having is in writing the function and wrapping it around model.frame. Consider the following data frame: dat <- data.frame(A = runif(10), B = runif(10), C

Dear all, How can I run a constrained correspondence analysis with the following data: 15 animals were measured repeatedly month-wise (over to 2 years) according to ther diet composition (8 food categories). our data.frame looks like this: food 1 2 ... 8 sex season year animal freq 12 8 ... 1 0 summer 2011 1 freq 0 7 ... 0 1 winter 2011 1 ... freq 0 7 ... 0 1 spring 2011 15 We

Hi list I loaded the content of a file dureetasks.txt to variable a. This file contains an interger per line. when I print a vector, it displays correctly. however, when I try to print the histogram, I get this error message > a=read.table("dureetasks.txt") > hist(a) Error in hist.default(a) : 'x' must be numeric Can you help please? regards -- PhD candidate in

Dear all I try to format labels on axis as standard ones does not look well, but I am not able to make axis.POSIXct to work. here is an example x<-seq(as.Date("2008-1-1"), as.Date("2009-9-6"), by="2 months") y<-rnorm(11) plot(x,y) #not very sophisticated axis plot(x,y, axes=F) axis(2) axis.POSIXct(1, at=x, format="%m/%Y") # no axis The same is on

There is a R function to calculate weighted mean : weighted.mean() under stats package. Is there any direct R function for calculating weighted variance as well? [[alternative HTML version deleted]]

Dear r-Community, Step1: I would like to calculate a NMDS (package vegan, function metaMDS) with species data. Step2: Then I want to plot environmental variables over it, using function envfit. The Problem: One of these environmental variables is cos(EXPOSURE). But for flat releves there is no exposure. The value is missing and I can't call it 0 as 0 stands for east and west. Therefore I

Dear List, I have a data set stored in the following format: > head(dat, n = 10) id sppcode abundance 1 10307 10000000 1 2 10307 16220602 2 3 10307 20000000 5 4 10307 20110000 2 5 10307 24000000 1 6 10307 40210000 83 7 10307 40210102 45 8 10307 45140000 1 9 10307 45630000 1 10 10307 45630600 41 >

Dear all, Are there any functions in R to reverse the order of the string. smth like reverse("abc") to get "cba"? Thanks a lot. [[alternative HTML version deleted]]

Dear list, Following advice posted to this list a while back by Prof Ripley [1], I have been trying to draw a per mille character [2] in an axis label. This should give the correct character: plot(1:10, ylab = "\u2030") but all I get is '"S'. I'm running linux (FC5) and have fonts installed that have the correct character (viewed in the Gnome character map at least).

Hi Sir I have made a data set having 23 stations of rainfall. when I use the attach function to approach indevidual stations then following error occurr. *>attach(data)* *>S.Sharif #S.Sharif is the station name which has 50 data values* *Error: object "S.Sharif" not found* Now how to solve this problem. Thank You Regards -- AMINA SHAHZADI Department of Statistics GC

Hi, I find it a bit annoying that aggregate.default forces the returned object to loose the 'name' of the variable aggregated, replacing it with 'x'. A brief example: > dat <- data.frame(A = runif(100), B = rnorm(100), + Group = gl(4, 25)) > with(dat, aggregate(A, by = list(Group = Group), FUN = mean)) Group x 1 1 0.6523228 2 2

I have source a script running with no errors that ends in these lines: ============== for (i in 2:3) win.graph() dev.set(2) levelplot(avg ~ procs * size, rrt, drape=T, colorkey=T,main="ops/sec",scales=list(x=list(tick.number=15),y=list(log=10))) dev.set(3) levelplot(avg ~ procs * size, ort, drape=T, colorkey=T,main="ops/sec",scales=list(x=list(tick.number=15),y=list(log=10)))

Dear list, Is it intentional that example() opens a new device and leaves it in a mode where it asks (prompts) the user to Press return to see the next plot for *all* subsequent plots on that device. For example; with an already opened device, example() works as I would expect: > par("ask") [1] FALSE > example(lm) .... OUTPUT CLIPPED .... lm> plot(lm.D9, las = 1) #

Dear list, I have been using and testing the alpha/betas of R 2.4.0, sources for which I got from svn.r-project.org/R/branches/R-2-4-branch, and which I compiled on my Fedora Core 5 machine. No problems there. I have just updated the source tree today and re-compiled - which also works just fine. When run, this gives: R version 2.4.0 Patched (2006-10-03 r39576) So I assume from that, that I am

Dear List, Apologies for this off-topic post but it is R-related in the sense that I am trying to understand what R is telling me with the data to hand. ROC curves have recently been used to determine a dissimilarity threshold for identifying whether two samples are from the same "type" or not. Given the bashing that ROC curves get whenever anyone asks about them on this list (and