[R] how to specify the multinomial distribution in R

p", "a", "u", "tau.u", "tau.s", "phi", "t1.p"), n.chains = 3, n.iter = 1000, OpenBUGS.pgm = "C:/Program Files (x86)/OpenBUGS/OpenBUGS322/OpenBUGS.exe")) Any guys have ever met this problem or know how to f

Re: [R] Join points with arrows based a TIME variable

Hi Uwe, I mistook the arguments for the function of arrows. I understand now. Thanks again. 2009/11/5 Uwe Ligges > > > zhijie zhang wrote: > >> Hi Uwe, >> It works. Why does the following argument generate different results? >> with(data, arrows(x[-length(x

Re: [R] Join points with arrows based a TIME variable

with(data, arrows(x[1], y[1], x[length(x)], y[length(x)])) 2009/11/4 Uwe Ligges > > > zhijie zhang wrote: > >> Hi, >> I have a data set with three variables,X Y and Time. X and Y are the >> coordinates of points, i want to join these points according to the Time >

[R] Join points with arrows based a TIME variable

Hi, I have a data set with three variables,X Y and Time. X and Y are the coordinates of points, i want to join these points according to the Time sequence using arrows? Demo Example data: > x<-c(1:6) > y<-c(1:6) > time<-c(6:1) > data<-cbind(x,y,time) > data x y time [1,] 1 16 [2,] 2 2

Re: [R] Creating a Clustered-Stacked Column Chart

ay be bad. Anyway, thanks a lot. 2009/10/11 Peter Ehlers > > > zhijie zhang wrote: > >> Thanks for your ideas and suggestions. I need to point out that most of >> us >> will create the Clustered-Stacked Column Chart in the matrix layout as >> David >> ga

Re: [R] Creating a Clustered-Stacked Column Chart

ived"),layout=c(1,2)) #missed the children's graph Thanks a lot. 2009/10/10 Peter Ehlers > I think you're missing the point. David _did_ show you how > to create a graph showing 4 clusters of stacked barcharts. > If you want them side-by-side instead of in the ma

Re: [R] Creating a Clustered-Stacked Column Chart

t(title = "Survived"), horizontal=FALSE) > > There are adjustments available to the space between bars. > > barchart(Titanic, scales = list(x = "free"), > auto.key = list(title = "Survived"), horizontal=FALSE, > box.ratio=100) > > > > O

Re: [R] Creating a Clustered-Stacked Column Chart

gt; alternative method of displaying the data using lattice. > > --- On Fri, 10/9/09, zhijie zhang wrote: > > > From: zhijie zhang > > Subject: [R] Creating a Clustered-Stacked Column Chart > > To: r-h...@stat.math.ethz.ch > > Received: Friday, October 9, 2009, 5:3

[R] Creating a Clustered-Stacked Column Chart

Hi all, In R, is there some functions or ways to create a Clustered-Stacked Column Chart as the example in the following page http://peltiertech.com/Excel/ChartsHowTo/ClusterStack.html? I have browsed the R Graph Gallery (http://addictedtor.free.fr/graphiques/) and searched the R site, and di

[R] cannot allocate a vector with 1920165909 length

Dear Rusers, The error for the following was that it cannot allocate the vector of length 1920165909. a <- expand.grid(se1=0:100/100, sp1=0:100/100, se2=0:100/100, sp2=0:100/100, DR=0:100/100) How to solve it? Maybe setwd(dir) can, i am not very sure about it. Any ideas about it? [[alt

Re: [R] How should i change the SAS Codes into R Codes?

ht find Bob Meunchen's book "Muenchen, > R. A. (2008). R for SAS and SPSS Users (1st ed.). Springer." useful. A > shorter version is available as a pdf at http://rforsasandspssusers.com/. > The pdf is a very useful guide for anyone beginning to use R. > > --- On Mo

[R] How should i change the SAS Codes into R Codes?

Dear R users, I have a SAS codes with several loops in it, and i hope to use R to do the same task. The SAS codes are as follows, /*to generate the dataset*/ DATA Single_Simulation; DO se=0 to 1 by 0.01; DO sp=0 to 1 by 0.01; DO DR=0 to 1 by 0.01; TR=(DR+sp-1)/(se+sp-1+1.0e-12); Adj

Re: [R] Does R have function or method to finish my task?

Healthcare > greg.s...@imail.org > 801.408.8111 > > > > -Original Message- > > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > > project.org] On Behalf Of zhijie zhang > > Sent: Thursday, March 05, 2009 8:56 AM > > To: r-h...@stat.mat

[R] Does R have function or method to finish my task?

Dear all, Does R have function or method to finish the following task. Say i have a dataset. The response variable is Y, and the indepedent variables are X1, X2, X3, and YEAR. See an example. Y X1X2X3 X4YEAR 13.4 2.8 3.5

Re: [R] Where to find the source codes for the internal function in stats package

> > - Original Message - From: "Duncan Murdoch" > To: "zhijie zhang" > Cc: > Sent: Saturday, January 17, 2009 6:00 AM > Subject: Re: [R] Where to find the source codes for the internal function > in stats package > > > On 17/01/2009

Re: [R] glmer documentation

r-sig-mixed-mod...@r-project.org maybe more appropriate for your question. You may find answers there. 2009/1/17 Raphaelle > > Hello, > > I am fitting a gmler using poisson, and I was looking for a documentation > to > interpret correctly the output. I'm quite a beginner with these kind of > mod

[R] Where to find the source codes for the internal function in stats package

Dear all, I want to see the source codes for "dchisq(x, df, ncp=0, log = FALSE)", but cannot find it. I input "dchisq" in the R interface, and then enter, the following message return: > dchisq /*/ function (x, df, ncp = 0, log = FALSE) { if

[R] Questions on the results from glmmPQL(MASS)

week>2) -1.61320.478549-3.370.00150.05 -2.5747 -0.6516 -9.35E-7 sigma1.53010.963249 1.590.11860.05 -0.40543.4656 -2.42E-6 -- With Kind Regards, oooO: (..): :\.(:::Oooo:: ::\_)::(..):: :::)./::: ::::::(_/:::: ::

[R] Needs suggestions for choosing appropriate R packages

specified in X1 as above, Which package/functions is possible? Thanks a lot. -- With Kind Regards, oooO: (..): :\.(:::Oooo:: ::\_)::(..):: :::)./::: ::(_/ : [***] ZhiJie Zhang ,PhD Dept.of

Re: [R] Generating series of distributions with the same skewness and different kurtosis or with same kurtosis and different skewness?

t; Renmin University of China, Beijing, 100872, China > > > > On Tue, Sep 23, 2008 at 10:59 PM, zhijie zhang <[EMAIL PROTECTED]> wrote: > > Dear R users, > > I hope to explain the concepts of skewness and kurtosis by generating > > series of distributions with s

[R] Generating series of distributions with the same skewness and different kurtosis or with same kurtosis and different skewness?

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 answ

Re: [R] High resolution figures for a paper?

chive is > available here: > https://stat.ethz.ch/pipermail/r-sig-ecology/2008-May/thread.html > and the thread is very near the bottom. Very detailed recommendations > have been provided. (Also see the past two days, which are not of > course in the May archive.) > > Sarah >

[R] High resolution figures for a paper?

Dear Rusers, My manuscript has been conditionally accepted recently. The problem to generate the high resolution figures in R for the manuscript cannot be solved by me. The journal editor ask me to generate the figures with a minimum resolution of 500 dpi. I have tried the *menu-driven method* t

Re: [R] Great difference for piecewise linear function between R and SAS

-0.3202))3-11.216 2.861 - 3.920 8.86e-05 *** Thanks again. On Mon, Mar 24, 2008 at 8:08 PM, Duncan Murdoch <[EMAIL PROTECTED]> wrote: > On 24/03/2008 7:06 AM, zhijie zhang wrote: > > Dear Murdoch, > >"Compare the predictions, not the coefficients.", this

Re: [R] Great difference for piecewise linear function between R and SAS

so that i can go to find some materials? Thanks a lot. On 3/24/08, Duncan Murdoch <[EMAIL PROTECTED]> wrote: > > On 24/03/2008 5:23 AM, zhijie zhang wrote: > > Dear Rusers, > > I am now using R and SAS to fit the piecewise linear functions, and > what > > sur

[R] Great difference for piecewise linear function between R and SAS

Dear Rusers, I am now using R and SAS to fit the piecewise linear functions, and what surprised me is that they have a great differrent result. See below. #R code--Knots for distance are 16.13 and 24, respectively, and Knots for y are -0.4357 and -0.3202 m.glm<-glm(mark~x+poly(elevation,2)+bs(di

[R] What is the correct model formula for the results of piecewise linear function?

Dear friends, I used the B-spline (degree=1) method to fit the piecewise linear function and the results are listed below. m.glm<-glm(mark~x+poly(elevation,2)+bs(distance,degree=1,knots=c(16.13,24)) +bs(y,degree=1,knots=c(-0.4357,-0.3202 )),family=binomial(logit),data=point) summ

Re: [R] Specify a correct formula in R for Piecewise Linear Functions?

*example*( pmin ) is useful. I will use it from now on. Thanks a lot. On Jan 3, 2008 2:50 AM, Charles C. Berry <[EMAIL PROTECTED]> wrote: > On Thu, 3 Jan 2008, zhijie zhang wrote: > > > Some developments with confusions. I tried the spline method and dummy > > variabl

Re: [R] Specify a correct formula in R for Piecewise Linear Functions?

. <13,13~25, >25. Q3:"fit <- glm( y ~ pmax(x,20)+pmin(x,20), family=binomial)" is good. But if i divide x into three or more parts, how should i specify it in this way? Hope somone can help.Thanks a lot. On Jan 2, 2008 11:58 PM, Thomas Lumley <[EMAIL PROTECTED]> wrot

[R] Specify a correct formula in R for Piecewise Linear Functions?

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 t

Re: [R] help on ROC analysis

Thanks. library(ROCR) was used finally. It also automatically generate a plot beside the value of AUC. On Dec 31, 2007 11:38 PM, Frank E Harrell Jr <[EMAIL PROTECTED]> wrote: > zhijie zhang wrote: > > Dear all, > > Some functions like 'ROC(Epi)' can be used to

[R] help on ROC analysis

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

Re: [R] Two repeated warnings when runing gam(mgcv) to analyze my dataset?

)? > > n Thursday 13 December 2007 17:46, zhijie zhang wrote: > > Dear all, > > I run the GAMs (generalized additive models) in gam(mgcv) using the > > following codes. > > > > m.gam > > > <-gam(mark~s(x)+s(y)+s(lstday2004)+s(ndvi2004)+s(slope)+s

[R] Two repeated warnings when runing gam(mgcv) to analyze my dataset?

Dear all, I run the GAMs (generalized additive models) in gam(mgcv) using the following codes. m.gam <-gam(mark~s(x)+s(y)+s(lstday2004)+s(ndvi2004)+s(slope)+s(elevation)+disbinary,family=binomial(logit),data=point) And two repeated warnings appeared. Warnings$B!'(B 1: In gam.fit(G, family = G

[R] Probelms on using gam(mgcv)

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

[R] Hep on using GAM() in R

Dear friends, I met some problem on using GAM() function in R. Any help or suggestions are greatly appreciated. # My programs and problems are list below# library(splines) library(gam) point<-read.csv("d:/gam.csv",sep=",",header=TRUE) #read the data gam.opt<-gam(mark~lo(x,y,span=0.2)+lo(lst

[R] question on image() function?

Dear friends, My dataset is like the following: xy mcpvalue 0.4603578 0.6247629 1.001 0.4603715 0.62477881.001 0.4603852 0.6247948 1.001 0.4110561 0.5664841 0.995 The x and y variables are unsorted. I use the functi

[R] Anybody has ever met the problem to add a legend to a figure generated by image()?

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 le

[R] How to add legend for image()?

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(

Re: [R] return(x=x,y=y,prob=prob) hasn't been used in R now?

; On Sun, 23 Sep 2007, Peter Dalgaard wrote: > > > zhijie zhang wrote: > >> Dear friends, > >> Now, when i use the argument return(x=x,y=y,prob=prob) , R displays > the > >> waring message: > >> Warning message: > >> The return value for

Re: [R] return(x=x,y=y,prob=prob) hasn't been used in R now?

It works well. Thanks a lot. On 9/23/07, Christian Ritz <[EMAIL PROTECTED]> wrote: > > Hi! > > Use a list structure for all the components you want to have returned > by the function: > > > return(list(x=x, y=y, prob=prob)) > > > > Christian > -- With Kind Regards, oooO: (..):::

[R] return(x=x,y=y,prob=prob) hasn't been used in R now?

Dear friends, Now, when i use the argument return(x=x,y=y,prob=prob) , R displays the waring message: Warning message: The return value for multiple variables wasn't used in: return(x = x, y = gy, prob = prob) I used the methods of "help.search("return")" and "?return" to get some help, but di

[R] How to explain the meaning of mu in the variance function of GLMs?

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 f

Re: [R] What's the corresponding function in R for lo() function inS-PLUS?

lEngineering.com > > > > > > -Original Message- > > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] > On > > Behalf Of zhijie zhang > > Sent: Monday, September 17, 2007 1:19 PM > > To: [EMAIL PROTECTED] > > Subject: [R] What's the

[R] What's the corresponding function in R for lo() function in S-PLUS?

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

[R] Calculate the angles for a point dataset?

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?

[R] handle dates in R?

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.mmdd(), but cannot get it. Anybody can show me how to handle the date data? Thanks very much! My dataset like: