similar to: return(x=x,y=y,prob=prob) hasn't been used in R now?

Dear All, The function I wrote can run well with the small data, but with the large data, the function runs very very slowly. How can I correct it? Thank you very much. My function as below: a<-c(1:240) b<-c(1:240) l=function(a,b){ v=0 u=0 uv=0 v[1]=0 u[1]=0 uv[1]=0 for (i in 1:(length(s)-1)){ v[i]<-((gx[[i]][b,(gx[[i]][a,1]+1)])-(gx[[i]][a,gx[[i]][a,1]+1]))/(gx[[i]][a,gx[[i]][a,1]+1])

Dear all, I have a big data frame: str(data1) 'data.frame': 18272 obs. of 11 variables: $ tag : int 100001 100002 100003 100005 100007 100008 100009 100011 100012 100014 ... $ sp : Factor w/ 18 levels "acassp","acocar",..: 13 5 7 14 14 18 3 11 13 10 ... $ gx : num 20 10 35 68 88 63 123 115 137 136 ... $ gy : num 30 25 24 1 10 40 45 25 23 45 ...

Hello, I have written a function in order to analyse gaze paths. It works with the test data but when I try to apply the function to a data frame that stores "the real data" in columns I receive the error message that the " In if (pp > 1) { : condition has length > 1 only the first element will be used " I interpret this error message as saying that only the first

Dear all, Following the help from gam(mgcv) help page, i tried to analyze my dataset with all the default arguments. Unfortunately, it can't be run successfully. I list the errors below. #m.gam<-gam(mark~s(x,y)+s(lstday2004)+s(slope)+s(ndvi2004)+s(elevation)+s(disbinary),family=binomial(logit),data=point)

Dear Rusers, I have some data in .csv file like "2004-8-1" and "2004-10-1", and i need to convert them into days from the origin (January 1, 1960). I have tried the function date.mmddyyyy(), but cannot get it. Anybody can show me how to handle the date data? Thanks very much! My dataset like: time 2004-8-1 2004-10-1 2001-9-1 2002-9-1 -- With Kind Regards,

Estimate region of highest probabilty density Dear R-community I have data consisting of x and y. To each pair (x,y) a z value (weight) is assigned. With kde2d I can estimate the densities on a regular grid and based on this make a contour plot (not considering the z-values). According to an earlier post in the list I adjusted the kde2d to kde2d.weighted (see code below) to estimate the

Dear friends, The following is an example to explain my question. I want to get a legend which will show the z-values according to different colors in image() function. x<-sort(runif(10)) #x-coordinates y<-sort(runif(10)) #y-coordinates z<-matrix(runif(100),nrow=10) #attributes values image(x,y,z,col=gray((6:3)/6)) # legend(x,y,legend=z,col=gray((6:3)/6)) #error. the colors should

Dear R users, I hope to explain the concepts of skewness and kurtosis by generating series of distributions with same skewness and different kurtosis or with same kurtosis and different skewness, but it seems that i cannot find the right functions. I have searched the mailing list, but no answers were found. Is it possible to do that in R? Which function could be used? Thanks a lot. --

Dear friends, My dataset is like the following: x y mcpvalue 0.4603578 0.6247629 1.001 0.4603715 0.6247788 1.001 0.4603852 0.6247948 1.001 0.4110561 0.5664841 0.995 The x and y variables are unsorted. I use the function image(x,y,mcpvalue) to generate a plot, but the error is that "increasing 'x' and

Dear friends, In S-PLUS, we can use the following argument, but not in R. mode12 <- gam(score1 ~ lo(latitude) + lo(longitude)) I searched the help in S-PLUS, it says lo() Allows the user to specify a Loess fit in a GAM formula, but i didn't find the correponding function in R. Anybody knows how to do the similar task in R? Thanks very much. -- With Kind Regards, oooO:::::::::

Dear R users, We have some programs for the specific task in our research, but they were very commonly used. We want to make some functions for them, anybody can recommend any good tools to facilitate us to create R functions even without going deep into the theories of R functions ? Any suggestions or help are greatly appreciated. -- With Kind Regards, oooO::::::::: (..):::::::::

Dear all, In an attempt to draw fill patterns in grid graphics, I have encountered a behavior of grobX that I cannot understand from the documentation. Consider this, library(grid) ## gTree g1 <- gTree(children=gList( rectGrob(0.5,0.5, width=unit(0.8,"npc"), height=unit(2,"cm")), circleGrob(r=0.3)), vp=viewport(0.5,0.5)) ##

Dear all, I have two variables, y and x. It seems that the relationship between them is Piecewise Linear Functions. The cutpoint is 20. That is, when x<20, there is a linear relationship between y and x; while x>=20, there is another different linear relationship between them. How can i specify their relationships in R correctly? # glm(y~I(x<20)+I(x>=20),family = binomial, data =

Dear friends, Anybody has ever met the problem to add a legend to a figure generated by image()? I have three variables,x,y and z. x and y are the coordinates, and z is the third values. we can use image(x, y, z,...) to generate a figure according to the z-values, but the problem is the figure legend. How can the legend be added to a figure generated by image()? Note that filled.contour() can

Dear list, Following a recent enquiry, I've been playing with the idea of creating a colour gradient for a polygon, using the Grid package. The idea is to draw a number of stripes of different colours, using the grid.clip function. Below is my current attempt at this, library(grid) rotate.polygon <- function(g, angle=0){ # utility function, works fine matR <- matrix(c(cos(angle),

Dear all, Some functions like 'ROC(Epi)' can be used to perform ROC analyssi, but it needs us to specify the fitting model in the argument. Now i have got the predicted p-values (0,1) for the 0/1 response variable using some other approach, see the following example dataset: id mark predict.pvalue 1 1 0.927 2 0 0.928 3 1 0.928 ..................

Dear friends, I need the get the 2.5% and 97.5% quantile from each row of a matrix, how should i get it? BTW, i can get the min/max value from each row of a matrix, using the following programs, is there an easy function to do it? simmin<-matrix(NA,nrow=47,ncol=1) for (i in 1:47) { simmin[i,]<-min(datas[i,]) } Thanks for your help. -- With Kind Regards, oooO:::::::::

Dear Rusers, When i start R, there always the following work to do first, how should i cancel it? *--- Please select a CRAN mirror for use in this session ---* I don't know why it does so, maybe i have done something unintentionally. Thanks. -- With Kind Regards, oooO::::::::: (..)::::::::: :\.(:::Oooo:: ::\_)::(..):: :::::::)./::: ::::::(_/:::: :::::::::::::

Dear Rusers, I'd like to take the cases of cancer of the larynx in chorley(spatstat) to explain my question. I want to join the points of cancer of the larynx with the disused industrial incinerator to generate lines, and then calculate the angles of these line comparing the horizontal line? How can i get it? It seems to be difficult to get the angles Thanks a lot. -- With Kind

Dear R friends, When fitting GLMs in R, we may need to specify the variance function to do our analysis. I had thought it's the mean value, but it seems not. Could anybody expain the correct meaning of *mu* in the variance function of GLMs? The following content is from the R-hlep. variance for all families other than quasi, the variance function is determined by the family. The quasi