search for: 2km

Dear all, I have a set of data that is 2 columns wide and 35,000 rows long (see extract below). The first column contains codes for tetrads (2km x 2km squares within NW England) and the second column has the names of moss species present within each tetrad. I wish to convert this into a matrix where rows are labelled by 'Tetrad', columns are labelled by 'Moss species' and each cell contains a 1 or 0 to denote the presence o...

Dear R support network I have a problem that is driving me crazy. I have a dataframe with about 74000 landscapes which I call "land". A landscape is a 2km -by- 2km square. Land has three columns: land$lat, land$long, and land$description. The last one holds a NON-unique (integer) description of each landscape. There are maybe 100 distinct descriptions. Identifying landscapes that have identical descriptions gives equivalence classes. I am interes...

Dear R gurus, I have some climatic data for a region of the world. They are monthly averages 1950 -2000 of precipitation (12 months), minimum temperature (12 months), maximum temperature (12 months). I have scaled them to 2 km x 2km cells, and I have around 75,000 cells. I need to feed them into a statistical model as co-variates, to use them to predict a response variable. The climatic data are obviously correlated: precipitation for January is correlated to precipitation for February and so on .... even precipitation a...

Hi All, Does anybody know of any existing functions that will calculate distance inland from a coastline? It's possible to test if a lon,lat location is land or sea using map.where(), but I need to add a buffer to this of say 2km, to allow for points that are just on the coast, and below the resolution of the worldHires database. I'm working with a marine mammal satellite telemetry dataset and wish to filter out spurious locations on land. On another issue, does anybody know of any free vector map datasets that are mo...

...card, 2 red modules fitted... Both channels work fine under 1.2/Zaptel. With 1.4/DAHDI both channels still work OK, but only for one line - the 2nd line causes it to refuse to dial-out no matter which port it's plugged into. The Lines are bog-standard BT analogue lines and we're about 2Km from the exchange. Both sound good to me and dial out OK with a test phone connected to them, but only one will dial-out via the PBX. This is what I see: [Jan 1 05:14:14] WARNING[1200]: app_dial.c:1237 dial_exec_full: Unable to create channel of type 'DAHDI' (cause 0 - Unknown) == E...

...n is an object of class "varest" generated by VAR() -exactly the class of object that var.2c is: > class(var.2c) [1] "varest" May it be that I'm using the latest version of R -in Windows- (ie 2.15.3 ), perhaps? Jose Wrap Up and Run 10k is back! Also, new for 2013 ? 2km intergenerational walks at selected venues. So recruit a buddy, dust off the trainers and beat the winter blues by signing up now: http://www.ageuk.org.uk/10k Milton Keynes | Oxford | Sheffield | Crystal Palace | Exeter | Harewood House, Leeds |...

I appreciate the help I've been given so far. The issue I face is that the data I'm working with has 53000 rows, so in calculating distance, finding all recids that fall within 2km and summing the population, etc. - a) takes too long and b) have no sense of progress. Below is a loop that reads each recid one at a time, calculates the distance and identifies the recids that fall within 2 km. It iterates through all records successfully. Where I'm stuck is how to get the...

...g.uk> Twitter @jose.iparraguirre@ageuk Tavis House, 1- 6 Tavistock Square London, WC1H 9NB www.ageuk.org.uk<http://www.ageuk.org.uk> | ageukblog.org.uk<http://ageukblog.org.uk/> | @ageuk<http://twitter.com/ageukpa>campaigns Wrap Up and Run 10k is back! Also, new for 2013 – 2km intergenerational walks at selected venues. So recruit a buddy, dust off the trainers and beat the winter blues by signing up now: http://www.ageuk.org.uk/10k Milton Keynes | Oxford | Sheffield | Crystal Palace | Exeter | Harewood House, Leeds |...

In R I have 3 variables x,y and z that contain: | x y z [1,] 1 12 122 [2,] 1 13 113 [3,] 1 11 111 [4,] 3 13 123 [5,] 4 14 134 [6,] 5 15 155| I am creating a coplot with the following syntax: |given.depth<- co.intervals(z, number=3, overlap=.1) coplot(x~y|z,pch=19,col="red",cex=1.5,given.values=given.depth)| The resulting plot however has has bars in the

Dear useRs, In a thesis, I found a mention of the use of pairwise deletion in linear regression and GLM (binomial family). The author said that he has used R to do the statistics, but I did not find the option allowing pairwise deletion in both lm and glm functions. Is there somewhere a package allowing that ? Thanks, Arnaud [[alternative HTML version deleted]]

Dear useRs, With the following set of command i managed to create the copula density ================================================ >library(MASS) >library(evd) >X<-matrix(sample(1:3000), ncol=2) >U=cbind(rank(X[,1])/(nrow(X)+1),rank(X[,2])/(nrow(X)+1)) > mat1=kde2d(U[,1],U[,2],n=150) image(mat1$x,mat1$y,mat1$z,col= rev(heat.colors(1500)),xlab="",ylab="")

I have a data. frame to which you want to change the names to some of their columns. For example: > seba <- data.frame ('constant' = 3, 'numbers' = 1: 10, 'letters' = LETTERS [1:10], otros = 2:11) List their names: > names (Seba) [1] "constant" "numbers" "letters" I want to rename c the column called constant and b the column

Hello, I have a problem regarding calculation of Cumulative Incidence Function. The event of interest is failure of bone-marrow transplantation, which may occur due to relapse or death in remission. The data set that I have consists of- lifetime variable, two indicator variables-one for relapse and one for death in remission, and the other variables are donor type (having 3 categories), disease

Hi, I am trying to draw a dotted circle by using symbols command: symbols(0, 0, circles = 1, inches = FALSE, fg = 'gray', lty=2) But it gives me a solid line. I was just wondering whether am I making a mistake? Ved [[alternative HTML version deleted]]

Hi all, I want to plot a facet plot with column names as x and column values as y. One plot for each row. here is part of my dataset: Gene T0h T0.25h T0.5h T1h T2h T3h T6h T12h T24h T48h NM_001001130 68 95 56 43 66 62 68 90 63 89 NM_001001144 0 1 4 0 1 1 1 4 1 2 NM_001001152 79 129 52 50 24 45 130 154 112 147 NM_001001160 1 1 2 0 0 1 0 0 0 1 NM_001001176 0 0 0 0 0 0 0 0 0 0 NM_001001177

Dear all, I have a matrix with two columns: "Names" and "Values" In names: there are 4 groups they are, CK113234, CK116296, CK116292 and CK114042 I want to *sort values* (decreasing order) based on each group and take average of the *top two numbers* in each of the groups. > dput(x) structure(list(Names = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 5L, 5L, 5L, 5L, 2L,

Dear R Helpers, I'm struggling with a data preparation problem. I feel that it is a quite easy task but I don't get it done. I hope you can help me with that. I have a data frame looking like this: ID <- c(1,1,1,2,2,3,3,3,3) T <- c(1,2,3,1,4,3,5,6,8) x <- rep(1,9) df <- data.frame(ID,T,x) >df ID T x 1 1 1 1 2 1 1 3 1 2 1 1 2 4 1 3 3 1 3 5 1 3 6 1 3 8 1 I want to

Hello all, this is not really an * question but it is somehow related, i am trying to develop a working proposal for cheap and quick telephony services using Voip running over *. By running a wireless network (over 802.11 a/b/g devices), i plan to be able to reach customers directly with eithe table top or handheld 802.11 sip enabled phones. But the disadvantage is that how do i power each radio

Hello, How can I judge if a string is in a group of string? For example, I would like to have if (subpool in pool){ }else{ } Where > pool = c("s1","s2") > subpool = c("s1") How can I write the "subpool in pool" right in R? Thanks very much! Cheers, Rebecca ---------------------------------------------------------------------- This message,

Hi R bayesians, I need an advise how to resolve the two different estimates applying a traditional glm (TG) and a bayes glm (BG), and different results depending on the data formats of response data and the prior specs using bayesglm in R. I'm not familiar with bayes estimate and my colleague asked me to look into this because the EPA from France reported a quite different estimates for