similar to: how do I get to be a member?

Dear List, A student in the Department where I work would like to produce a graphic similar to this one: http://image.guardian.co.uk/sys-files/Guardian/documents/2009/09/16/Public_spending_160909.pdf Does anyone know if the figure in the pdf can be generated in a specific software application for example? Any suggestions would be most gratefully received by the student concerned. Many thanks,

It seems division with numbers bigger than 10 000 000 doesn't work 20000000/21 [1] 952381 > 55555555/23 [1] 2415459 Thank you

Hi all, I am dealing with paleoceanographic data and I have a C14 time serie and one other variable. I would like to perform a spectral analysis (fft or wavelet) and plot it. Unfortunately I don't know the exact script to do this. Does anybody could send me an example to perform my spectral analysis ? I Thank you David Changez de tête et de tenue tous les jours si vous le voulez !

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 --

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

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]]

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 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 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]]

Hi folks, Where can I find document re "how to read anova output"? Google found many of them. But seemingly non of them can explain to me following output:- > tabA = c(5.67, 5.67, 5.55, 5.57) > tabB = c(5.75, 5.47, 5.43, 5.45) > tabC = c(4.74, 4.45, 4.65, 4.94) > tabs = data.frame(tabA, tabB, tabC) > tablets = stack(tabs) > anova(lm(values ~ ind, data =

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)))

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

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 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

Greetings, I have what I hope is a simple question. I would like to change my contour interval on the vis.gam( plot.type="contour") in the mgcv package. Is this a situation where I need to modify the function or is there a default value I can change? Thanks

Hi all, why factorial(150) shows the error out of range in 'gammafn'? I have to calculate the number of subset formed by 150 samples taking 10 at a time. How is this possible? best [[alternative HTML version deleted]]