search for: xseq

Hi all, I have a data frame with sequences of x and y from a map. I would like to know it both ways: 1. How to make a matrix from that; 2. how to make a data frame of all points in a map. Probably it is a silly question, but please tell me where to read about it or tell me how to do it. Miha Staut

Hi Peter, You didn't give a very specific example, but it seems to me that what you wish to do is not really complicated. I suppose you have created a table of sequences vs. say hyprophobicity, charge, etc..., something like... seq hydroph arom b0001 0.104762 0.000000 b0002 0.035122 0.065854 b0003 0.024193 0.070968 b0004 -0.096729 0.084112 b0005 -0.973469 0.091837 b0006

Hi all, I tried to produce some kriged surfaces with geoR (latest version). The size of the grid should be around 900 x 650 cells (what I find is not a very big grid), and the number of points is around 2500. The command krige.conv stopped after arround 5 min saying it can not allocate a vector with around 1.5 billion units. Sounds reasonable. Is there a workaround? How would I partition the

Dear All I would like to use filled.contour to plot something with arrows on it. I did it in this way: filled.contour(1:nx,1:nz,u,col=gray(rev((0:20/20)))) for (j in zseq ) for (i in xseq) arrows(i, j, i+u[i,j],j+w[i,j],length=angleng,angle=angarrow,code=2, lwd=1,col="black") The peoblem is that the arrows use the whole device area, including the key on the right hand side of the filled contour. Does anyone know how to work out this problem? Many thanks for your kind rep...

vector with a sequence starting from 1 to sum of all the values of the raster data and length(the same length as the sum). In another data frame (G) I have got a vector with the environment variables (G$xseq and G$yseq). How can I combine the x and y coordinates with tha raster data to form a data.frame? length(G$xseq)=3099 length(G$yseq)=5032 rastername=rst length(rst)=(3099*5032) Next question: I would like to broaden every unit in rst in the following way: rst2<-c(rep(rst[1:1], times=rst[1:1])...

...? >nomogramCF = nomogram(cph_age6_40avp4, +????????????????????lp.at= seq(-10,10,by =0.5),lp = TRUE, +?????????????????????? +??????????????????????funlabel="5year survival", +??????????????????????fun=surv5.CFdata5,fun.at= c(0.01,seq(0.1,0.9,by=0.2),0.99)) Error inapprox(fu[s], xseq[s], fat, ties = mean) : ??needat least two non-NA values to interpolate ? ? I have fit similar nomograms based on cox ph models usingsimilar code, so I?m not sure what I?m doing wrong this time. ? I heard there might be some an issue because of a softwareupdate? ? Thanks, Annie [[alter...

....pressure, cholesterol, sex) > options(datadist='ddist') > > fit <- lrm(y ~ blood.pressure + sex * (age + rcs(cholesterol,4)), + x=TRUE, y=TRUE) > > par(mfrow=c(2,2)) > plot(fit) # Plot effects of all 4 predictors Error in value.chk(at, ix, xseq, if (plot.type == "curves") 100 else 40, : variable blood.pressure does not have limits defined by datadist I guess somthing goes wrong with function datadist > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mi...

...,lp) at.surv=c(0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9) at.med=c(120,80,60,40,30,20,15,10,8,6,4,2,0) nomogram(f,lp=F, fun=list(srv120, med),funlabel=c("120-mo Survival","Median Survival"),fun.at=list(at.surv, at.med)) I get a the following warning: Warning message: In approx(fu[s], xseq[s], fat) : collapsing to unique 'x' values However, a great nomogram is constructed. But then I try to stratify... Stratified Nomogram: f...

Hi to all, I'm trying to run Battery3 with Wine, but I'm experiencing some problems. 1. When I launch Battery (command line) I see these messages: --- fixme:win:RegisterDeviceNotificationW (hwnd=0x50062, filter=0x32fcb0,flags=0x00000004) returns a fake device notification handle! JackActivationCount::Signal value = 0 ref = 5 JackActivationCount::Signal value = 0 ref = 5

...ent to position dodge so it's now possible to accurately dodge things with different widths to their physical widths * added median summary * new examples: * logistic regression example in stat_smooth * bugs fixed: * evaluation of arguments to layer is no longer delayed * can use categorical xseq with stat_smooth * x and y axes named incorrectly (thanks to Dieter Menne for spotting this) * can now pass position objects to qplot * y jitter calculated correctly, and jittered data rescales axis now * removed silly legend from quantile plot * extra arguments not being passed on to geoms/st...

...ent to position dodge so it's now possible to accurately dodge things with different widths to their physical widths * added median summary * new examples: * logistic regression example in stat_smooth * bugs fixed: * evaluation of arguments to layer is no longer delayed * can use categorical xseq with stat_smooth * x and y axes named incorrectly (thanks to Dieter Menne for spotting this) * can now pass position objects to qplot * y jitter calculated correctly, and jittered data rescales axis now * removed silly legend from quantile plot * extra arguments not being passed on to geoms/st...

Hello, I have a question about how to plot a series of data. The folloqing is my data matrix of n > n 25p 5p 2.5p 0.5p 16B-E06.g 45379 4383 5123 45 16B-E06.g 45138 4028 6249 52 16B-E06.g 48457 4267 5470 54 16B-E06.g 47740 4676 6769 48 37B-B02.g 42860 6152 19276 72 35B-A02.g 48325 12863 38274 143 35B-A02.g 48410 12806 39013 175 35B-A02.g 48417 9057 40923

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

Dear R users, Does anyone have a function for putting a legend on an image plot? I couldn't locate an R equivalent of image.legend....has anyone written such a thing? kind regards andy --------------------------------------------------------------------- J. Andy Royle, U.S. Fish and Wildlife Service - Office of Migratory Bird Management; 11510 American Holly Drive , Laurel, MD

Full_Name: Roger Bivand Version: 0.99.0 OS: RH Linux 6.1 Submission from: (NULL) (158.37.60.152) I am working on an interface between R and the GRASS geographical information system, written in R, with no dynamically loaded code. I have written full examples, and tested then under R 0.90.1, both by entering example() for each function and R CMD check, both of which worked without problem. Under

...g package: 'Design' The following object(s) are masked from 'package:Hmisc': strgraphwrap The following object(s) are masked from 'package:survival': Surv > ols(Retro~rcs(log(Frequency),4),data=Data)->Data.ols > plot(Data.ols) Error in value.chk(at, ix, xseq, if (plot.type == "curves") 100 else 40, : variable Frequency does not have limits defined by datadist > datadist(Data.ols)->d Error in sort.int(x, na.last = na.last, decreasing = decreasing, ...) : 'x' must be atomic > is.atomic(Data.ols) [1] FALSE > is.recursive...