[R] package flexmix: models for ordinal data?

Hi Klaus, I saw you post “package flexmix: models for ordinal data?” on grokbase. I was wondering did you find the answer? Or what package did you use at the end? I really appreciate if you would like to help me. Best, Xin __ R-help@r-project.org

[R] pgmm error message, Error in terms.default(formula) : no terms component nor attribute

a.frame(data, index=c("id","year_begin")) y<-urban[,7] x1<-urban[,8] x2<-urban[,9] x3<-urban[,10] x4<-urban[,11] z1<-pgmm(formula=y~lag(y,0:1)+x2+x3+x4,data=urban, effect = "twoways", model = "twosteps") Error in terms.default(formula) : n

Re: [R] Help: Error when installing R 2.8.1 in IBM AIX system from source code

Hi Jeff, Thanks for your quick response. Let me describe more about my situation. 1. Patch file - Yes, I found patch didn't change any file, so I manually changed 4 files(R-2.8.1/configure, R-2.8.1/configure.ac, R-2.8.1/Makeconf.in, R-2.8.1/src/main/Makefile.in) by reading the patch file. 2. why

[R] Help: Error when installing R 2.8.1 in IBM AIX system from source code

w how to fix it. ../../../bin/R[206]: "${}": bad substitution make: The error code from the last command is 1. Stop. Thanks and Best Regards. De Xin Luan (æ ¾ å¾

[R] how to put barchart and line chart in the same plot in ggplot2

dear userR: I am trying to plot two dependent variables in the same plot in ggplot2. because these two variables have very different magnitude, I have to use a second Y axis. I hope one variable to be line and the other to be barchart. The x axis is continuous. Yet since I have to make barchart, I

[R] statistical model with censored independent variable

hello: Does anyone know any R function which handles statisitcal model when the independent variable is censored? I know survival package does the analysis for censored dependnent variable. thanks -- View this message in context: http://r.789695.n4.nabble.com/statistical-model-with-censored-inde

Re: [R] how to add row index based a categorical column

thank you everyone. how can I not be aware of the existence of ave()? I try the following: tapply(data$value, data$level, rank). However, I have a very difficult time merging the resulting rank variable back to the original data frame. thanks a lot! -- View this message in context: http://r.78

[R] how to add row index based a categorical column

hello, I have the following data manipulation issue. the following is the sample data: value level 4 A 5 A 2 A 10 B 9 B 34 B 100 C 34 C 101C. I hope to get the following result: value level rank 2 A 1 4 A 2 5 A 3 9

[R] How to extract information from the following dataset?

the measurements only from a specific day and hour, i.e. if I want measurements from Jan 27 2010 16:--:-- then I get 125,125,176,159,142,142,125,125,125,125,125,125. Many thanks!! -- Xin Zhang Ph.D Candidate Department of Statistics University of California, Riverside [[alternative HTML

Re: [R] a question on sqldf's handling of missing value and factor

Dear Gabor: I did not have time to look at this issue these few days. Now I saw your investigation. I am very grateful that you quickly identify the root cause of this. It seems that a little caution needs to be exercised when applying sqldf to text file with large number of blanks (I have no contr

Re: [R] a question on sqldf's handling of missing value and factor

Dear Mr. Grothendieck : thank you so much for your attention. You are the real expert here. the following is a mock text file: a b c aa 23 aaa 34 77 note that both b and c column contain missing value (blank) I save it under my C drive and use bo

Re: [R] a question on sqldf's handling of missing value and factor

I am sorry for posting the wrong source file. the correct source file is as follows: a b c aa 23 aaa 34.6 77.8 They are tab delimited but somehow could not be displayed correctly in browser. -- View this message in context: http://r.789695.n4.nab

[R] a question on sqldf's handling of missing value and factor

Dear subscribers: I am using the following code to read a large number of big text files: library(sqldf) tempd <- file() tempdx <- sqldf("select * from tempd", dbname = tempfile(), file.format = list(header = T, sep="\t", row.names = F)) The problem is: all my numberical variable become facto

[R] problem with installing packages

Dear: I am recnetly trying to install some libraries. However, I found this issue for both my laptop and desktop even I uninstall and install it again. I even can not update the R now. I wonder if you have the similar issue. Thakns! Xin > chooseCRANmirror() Warning message:

[R] Color Alaska in USA map

Hello: I have a problem when I tried to color the USA map with different colors. The following is my data (I only used the second column of data): alaska,1, 2 Hawaii, 0, 0 USA,5, 5 And here is my code: library("latticeExtra") library("mapproj") state<-read.table

Re: [R] Data frame self-join with no duplicates

; On Mon, Sep 27, 2010 at 4:29 AM, Xin Zhang wrote: >> >> Suppose I have the following data frame (df): >> >> Year    Price >> --- >> 2001    10 >> 2002    20 >> 2003    30 >> >> I would

[R] Data frame self-join with no duplicates

Suppose I have the following data frame (df): YearPrice --- 200110 200220 200330 I would like to produce another data frame like this: a.Yeara.Priceb.Yearb.Price 200220200110 200330

Re: [R] why I could not reproduce the Mandelbrot plot demonstrated on R wiki

(programming_language) where the nice Mandelbrot plot and sample codes are provided. i would appreciate your help if you can provide further hint. thanks Peter Dalgaard-2 wrote: > > On 09/24/2010 10:41 PM, xin wei wrote: >> >> I am trying to reproduce the nice loo

[R] why I could not reproduce the Mandelbrot plot demonstrated on R wiki

I am trying to reproduce the nice looking of Mandelbrot demonstrated by R wiki page by the following code: library(caTools)# external package providing write.gif function jet.colors = colorRampPalette(c("#7F", "blue", "#007FFF", "cyan", "#7FFF7F", "yel

Re: [R] how to generate a random data from a empirical distribition

hi, Frank: how can we make sure the randomly sampled data follow the same distribution as the original dataset? i assume each data point has the same prabability to be selected in a simple random sampling scheme. thanks -- View this message in context: http://r.789695.n4.nabble.com/how-to-gene

Re: [R] how to generate a random data from a empirical distribition

this is very insightful. sounds exactly like what I want to do. thanks. Frank. -- View this message in context: http://r.789695.n4.nabble.com/how-to-generate-a-random-data-from-a-empirical-distribition-tp2302716p2304346.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] how to generate a random data from a empirical distribition

good point. It seems to be important to investigate the nature of distribution. I might be too naive to assume that a "empirical probability distribution" would be automatically generated from a cloud of data points. -- View this message in context: http://r.789695.n4.nabble.com/how-to-gener

Re: [R] how to generate a random data from a empirical distribition

hi, Dennis: points well taken. it seems to be important to investigate the nature of distribution. I may be too naive to assume a "empirical probability distribution" would be computed from a could of data points -- View this message in context: http://r.789695.n4.nabble.com/how-to-generate-

[R] how to generate a random data from a empirical distribition

hi, this is more a statistical question than a R question. but I do want to know how to implement this in R. I have 10,000 data points. Is there any way to generate a empirical probablity distribution from it (the problem is that I do not know what exactly this distribution follows, normal, beta?

[R] integrate of function f(x,y,z) respect to x only

Dear: Anyone has experience to use R integrate a function with more than two dimensions, but you only need to integrate the function respect to one variable, say x. Because I also want to maximise the log-likelihood after integration of function respect to x. Thanks! Xin --- Dr.Xin Shi

[R] Connecting to MS Access database

Hi All, Can anyone please suggest me from where should I start to learn about 'how to connect to access db' ? How if someone has some written code and I can go over that to understand and make necessary changes... any help would be highly appreciated, -- Xin Ge. [[altern

Re: [R] help with reshape is needed again!

thank you very much for such a thorough solution! -- View this message in context: http://r.789695.n4.nabble.com/help-with-reshape-is-needed-again-tp2260640p2260943.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org

[R] help with reshape is needed again!

hi, folks: i need to transpose the following data: genetissuepatient1 patient2 patient3. - gene1 breast 10100 1 gene2 breast 20200 4 gene3 breast 3050

Re: [R] help for reshape function

I am afraid that your solution is not solving the problem. it seems that timevar="gene" just create the followings: GENESAMPLE value id 1.11 Kidney 3.69351 1 2.11 Kidney 5.42710 2 3.11 Kidney 5.26883 3 4.11 Kidney 2.88098 4 5.11 Kidney 4.68519

[R] help for reshape function

hi, everyone: i have a question on the reshape function. i have the following dataset : gene tissuepatient1 patient2 patient3. _ gene1 breast 10 20 50 gene2 breast 20 40 60 gene3 br

[R] panel.abline {lattice} help

5) to (5, 5) I'm unable to control the length of these lines, can anyone please help me with this, I would highly appreciate any help. (I tried using panel.curve, but of no use) # Plotting code: library(lattice) bwplot(c ~ b, data = x, panel = function(...) { panel.abline(v=2.5,

[R] how to parse out string separated by special character

hi, everyone: i have a data frame with one single character column as followings: DS_xxx_yyy_24hr_zzz_1 DS_xxx_yyy_24hr_zzz_2 DS_xxx_yyy_24hr_zzz_3 DS_xxx_yyy_24hr_zzz_4 DS_xxx_yyy_24hr_zzz_1 DS_xxx_yyy_24hr_zzz_2 DS_xxx_yyy_24hr_zzz_3 DS_xxx_yyy_24hr_zzz_4 DS_xxx_yyy_24hr_zzz_1 i need to parse o

[R] forest() in {metafor} :: edit labels

> PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Xin Ge. [[alternative HTML version deleted]] __ R-help@r-project.org mailing

Re: [R] Multiple plots; single x(y) labels

) http://support.sas.com/documentation/cdl/en/grstateditug/61951/HTML/default/images/multi-panel-data.gif * In the process I learnt so many other stuff, ~Xin On Fri, May 14, 2010 at 1:29 PM, Shi, Tao wrote: > Gurmeet, > > I think Xin is more interested in the common a

Re: [R] Multiple plots; single x(y) labels

x27;m using) par(mfcol=c(2,2)) plot(1:3, 4:6, type = "s", main = "", ylab="", xlab="") plot(2:5, 5:8, type = "s", main = "", ylab="", xlab="") plot(3:6, 6:9, type = "s", main = "", ylab="&q

Re: [R] Multiple plots; single x(y) labels

clearly before. Any comments? Thanks, Xin On Thu, May 13, 2010 at 12:51 PM, David Winsemius wrote: > > On May 13, 2010, at 12:35 PM, David Winsemius wrote: > > >> ?plot # ylim >> >> and you need to have the data in a form (before plotting) where you can >>

[R] Multiple plots; single x(y) labels

quot;, main = "", ylab="", xlab="") plot(x <- sort(rnorm(47)), type = "s", main = "", ylab="", xlab="") plot(x <- sort(rnorm(67)), type = "s", main = "", ylab="", xlab="") also

[R] what the problem could be if i am suddenly unable to add abline to the scatter plot?

I am doing very regular stuff like the following: attach(wtana) fm<- lm(Body.Wt.on.SD1~Heart.Wt, data=wtana) #fm<- lm(wtana$Body.Wt.on.SD1~wtana$Heart.Wt) lrf<- loess(Body.Wt.on.SD1~Heart.Wt, wtana) #lrf<- loess(wtana$Body.Wt.on.SD1~wtana$Heart.Wt) plot(Body.Wt.on.SD1,Heart.Wt) #plot(wtana$Bod

[R] Help: coxph() in {survival} package

lue, I would like to extract these from R output (for each data set). Can any one help please, thanks in advance, Xin P.S. OR, may be if anyone can suggest me some other function which can help me with this. [[alternative HTML version deleted]] __

Re: [R] finding weekly average...

Thanks Dennis and Daniel for your reply! On Wed, Apr 7, 2010 at 12:24 PM, Daniel Malter wrote: > > Assume you have suitable "week" indicator and your data is stored in a > data.frame named "data." Then you get the weekly averages by > > #simulate data > week <- rep(1:52,each=7) > x <- week+rnorm

[R] finding weekly average...

Hi All, I have a time series data with two continuous variables (say Var1 and Var2) for 4 years (***not continuous, do have some breaks because of missing data***). Something like this: Date Var1 Var2 12/01/2004 7 0 12/01/2004 0 0 12/01/2004 0 7 12/01/2004 7 0 12/01/2004 0 7 12/01/2004 0 7 12/02/

Re: [R] Three most useful R package

- ggplot2 - lattice - reshape Request to R-team: Please make ggplot2 and reshape as default packages in R, they are really helpful! Thanks, Xin On Fri, Mar 5, 2010 at 10:48 PM, kMan wrote: > (1) - nlme, lattice, stats > (2) - a usable large-file/out of memory regression package that abs

Re: [R] how to make R plot under Linux

thank you all your guys. You all show a great deal of tolerance on my ignorance. I installed reflection X on windows and problems solved! BTW, I am not with school anymore unfortunately and none of my local colleagues use either R or Unix.. -- View this message in context: http://n4.nabble

Re: [R] how to make R plot under Linux

hi, Guys: thank you so much for all the suggestion. Now I seem to be able to set up x11 forwarding in PUTTY. however, I still could not get plot and I get the following error msg: Error in function (display = "", width, height, pointsize, gamma, bg, : X11 I/O error while opening X11 connectio

Re: [R] how to make R plot under Linux

hi, Kevin and K.Elo: thank you for the suggestion. Can you be more specific on these? (like how exactly get into x-switch or man ssh). I am totally ignorant about linux and SSH:( Memory limitation forces me to switch from windows to Linux cluster. Xin -- View this message in context

Re: [R] how to make R plot under Linux

thank you for reply. I just type: hist(x) from SSH terminal, expecting a histogram to pop up like what i got under windows.instead I got the following error msg: Error in X11(d$display, d$width, d$height, d$pointsize, d$gamma, d$colortype, : unable to start device X11cairo In addition: Wa

[R] regularised incomplete beta function

Dear: I am trying to work out the regularised incomplete beta function in R. I searched google on this and found UCS library. But I can not find this in R packages. Does anyone have use this before or how to insert UCS in R? Many Thanks! Xin Xin --- Dr.Xin Shi Senior Lecturer in

Re: [R] re-ordering x-lables using barchart()

, scales = list(x = list(abbreviate = TRUE, rot=45, minlength = 5))) On Tue, Dec 8, 2009 at 10:35 PM, Xin Ge wrote: > Hi Gary, yes I'm using barchart() a lot these days... and was facing same > problem of re-ordering... > > > On Tue, Dec 8, 2009 at 10:28 PM, Gary Mil

Re: [R] re-ordering x-lables using barchart()

Thanks David, it worked! On Tue, Dec 8, 2009 at 5:03 PM, David Winsemius wrote: > > On Dec 8, 2009, at 4:42 PM, Gary Miller wrote: > > Hi R Users, >> >> I'm trying to re-order the "site names" ("Waseca", "Morris", ...). I'm >> using >> following code: >> >> libarry(lattice) >> > > # slip this cod

Re: [R] Forest Plot

Tel: +31 (0)43 388-2277 > School for Public Health and Primary Care Office Location: > Maastricht University, P.O. Box 616 Room B2.01 (second floor) > 6200 MD Maastricht, The Netherlands Debyeplein 1 (Randwyck) > > > Original Message > From: Xin Ge

[R] {Lattice} help.

Q20 AXR10 AXS20 AYP20 AYQ5 AYR20 AYS5 BXP 20 BXQ10 BXR20 BXS10 BY P30 BYQ50 BYR30 BYS50 Thanks, ~Xin [[alternative HTML

[R] R ANOVA question

coefficients. My Question is: How to get back the original month coefficients (including April's coef.). So that I can find seasonal effects (for e.g. I'll average Feb, Mar, Apr to get an estimate for "SPRING" season). Any suggestions on how can I do this? Thanks, Xin [

Re: [R] Forest Plot

r and std. deviation? So, a 95% confidence interval would be [estimate +/- 1.96*sqrt(variance of estimate)]. Am I correct in saying this? Thanks again, Xin On Sat, Dec 5, 2009 at 6:21 PM, Viechtbauer Wolfgang (STAT) < wolfgang.viechtba...@stat.unimaas.nl> wrote: > The figure that you linked to

[R] Forest Plot

e following forest plot: http://bm2.genes.nig.ac.jp/RGM2/R_current/library/metafor/man/images/big_plot.rma.uni_001.png Thanks, Xin [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-h

Re: [R] Help: barchart() {lattice}

Thanks David, it worked! On Fri, Dec 4, 2009 at 10:36 PM, David Winsemius wrote: > > On Dec 4, 2009, at 10:17 PM, Xin Ge wrote: > > @ Francisco: Thanks, it worked. >> >> >> @ All: I'm able to change the colors of legend using following code: >> >

Re: [R] Help: barchart() {lattice}

auto.key = list(points = FALSE, rectangles = TRUE, space = "right"), scales = list(x = list(rot = 45))) Thanks, Xin On Fri, Dec 4, 2009 at 8:04 PM, Francisco J. Zagmutt < gerifalt...@hotmail.com> wrote: > Hello Xin, > > Take a look at the examples under ?print.

[R] Help: barchart() {lattice}

c(1,6), stack = TRUE, auto.key = list(points = FALSE, rectangles = TRUE, space = "right"), ylab = "Barley Yield (bushels/acre)", scales = list(x = list(rot = 45))) par(mfrow=c(1,1)) Thanks, Xin [[alternative HTML version deleted]]

Re: [R] (Grouped + Stacked) Barplot

. > Statistical Data Center > Intermountain 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 Xin Ge > > Sent: Wednesday, December 02, 2009 8:2

Re: [R] (Grouped + Stacked) Barplot

Hello David, Can you please explain a little more, I could not exactly got it. Thanks, Xin On Wed, Dec 2, 2009 at 10:50 PM, David Winsemius wrote: > > On Dec 2, 2009, at 10:25 PM, Xin Ge wrote: > > Hi All, >> >> I have googled and tried finding if someone has ever tri

[R] (Grouped + Stacked) Barplot

t "Factor1" on x axis -- each level of this factor should have three (grouped) bars (one for each level of Factor2). Which will further be stacked by each level of Factor3. Can anyone guide me please, thanks, ~Xin [[alternative HTML version deleted]]

[R] pass an array from a perl script to a R script

pass it by argument? Thanks! Xin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible

[R] dereferencing in R

Hi, All, I passed an array reference to the R script and do not know how to do dereferencing in the R script. Anybody has some suggestion? Many thanks __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

Re: [R] R's documentation

use 'mean'. -Original Message- From: Prof Brian Ripley [mailto:rip...@stats.ox.ac.uk] Sent: Friday, May 29, 2009 12:28 PM To: Zheng, Xin (NIH) [C] Cc: r-help@r-project.org Subject: Re: [R] R's documentation Which version of R is this? From the BUGS section in NEWS for R 2.9.0

Re: [R] R's documentation

ch the exact point of its algorithm. Why not say more clearly that "make use of partial sorting and calculating mean of middle values". Maybe I'm too strict. It's not a really bug. Let's stop this thread. From: Kjetil Halvorsen [mailto:kjetilbrinchmannhalvor...@gmail.com] S

[R] R's documentation

coming or drive people away. Xin Zheng [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html a

[R] why warning appears when set both 'scale' and 'breaks' meanwhile in heatmap.2 of gplots

only one or the other. What does 'unpredictable results' mean? Could anyone explain that? Thanks a lot. Xin Zheng [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-he

Re: [R] how to avoid add 'X' before numeric colnames when read.table

substring(header, 2,) works for the purpose perfectly. Xin Zheng [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org

Re: [R] how to avoid add 'X' before numeric colnames when read.table

, May 28, 2009 3:40 PM To: Zheng, Xin (NIH) [C] Cc: r-help@r-project.org Subject: Re: [R] how to avoid add 'X' before numeric colnames when read.table ?read.table check.names=FALSE On Thu, May 28, 2009 at 3:25 PM, Zheng, Xin (NIH) [C] mailto:zheng...@mail.nih.gov>> wrote: Hi all, Is

[R] how to avoid add 'X' before numeric colnames when read.table

Hi all, Is there any way to keep numeric colnames as is? Any hint will be appreicated. Xin Zheng [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] problem with ci for lmer

<-lmer(y~x2+(1|x1)) ci(fit) Does anyone know the solution? In addition, I want to do the daily update for R. I run the code below and it is not working. svn log -v -r HEAD:\{`date +%Y-%m-%d -d'7 days ago'`\} https://svn.r-project.org/R Does anyone use t

[R] how to take away the same varible when I use "merge"

Dear: I am trying to merge two tables by a common variable. However, there are a few same variables which are in both of two tables. How can I take them away when I merge the two tables? Thanks! Xin __ R-help@r-project.org mailing

[R] any other fast method for median calculation

I am wondering if there is another function with better algorithm. Any hint? Thanks, Xin Zheng __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html a

[R] how to do trend.test for subset data

Dear: I try to do the trend.test {pastecs} for a subset data. Anyone have this experience? Many thanks! Xin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

[R] any package for connecting berkeley db in R?

Hi there, Is there any such package? I searched but found none. Thanks in advance. Xin Zheng [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] how to add labels to existing plot for the subset of data

Dear: I am trying to plot x against y for a particular subset of data, say z=1, and labelling data by another variable say, k. My data are: y<-c(69.7, 82.3, 66.3, 107.3, 90.1, 63.7, 82, 74.4, 61.7, 93.4, 73.4, NA, NA, 70.7, 67.7, NA, NA) x<-c(71.2, 82.6, 67.4, 107.1, 90.5, 66.7, 83.9, 73.9, 61

[R] how can I compare two vector by a factor

Hi, I used Wilcox.test to carry out mann whiteney test when paired=false. However, I want to see the comparison of two variables, e.g. pre and post, grouped by treatment. Anyone has this experience? Thanks! Xin [[alternative HTML version deleted

Re: [R] any function for package searching?

Thanks Wacek. Sorry that I didn't clarify my question. I want to search uninstalled packaged in repository. I just found RSitesearch could do that and open the results in system browser. It's not flexible, although it's better than none. Xin -Original Message- From: Wa

[R] any function for package searching?

Hi there, I wanted to search some packages related with 'genetics'. I know I can do it on CRAN webpage. I'm just wondering if there's some function in R could do that. Debian apt offers one way. It would be perfect if R has some builtin function like Debian apt. Appre

[R] how to read data from the different worksheets in the same Excel file

Dear All: I am trying to read the data from the different sheet but in the same Excel file. Does R provide this kind of command to do it? Anyone has experience on this? Thanks! Xin [[alternative HTML version deleted

[R] how to plot histogram plot and fitted distributions on the same graph

60, 65, 70, 75, 80, 85, 90, 95, 100, 105, 110, 115, 120, 125, 130, 135, 140, 145, 150, 155, 160, 165, 170, 175, 180, 185, 190, 200) Many Thanks! Xin [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https:

[R] keep function in stepAIC

Dear: Does anyone use keep function in stepAIC command? If so, could you give an example? I try to use this function to choose some variables in all of possible models. Many Thanks! Xin [[alternative HTML version deleted]] __ R

[R] multinomial losgitic regression--vglm()

working for my case or have alternative approach? Many thanks! Xin [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

[R] histogram plot default

Dear All: I am trying to plot a series of data using histogram plot. But I want to change the default of setting for histogram plot. For example, I want to set frequency plot starting with zero. Many Thanks Xin [[alternative HTML version deleted

[R] how to plots two pairwise data sets into a same graph

for "x-fits". But I want to plot x-fits against fits on the same graph with the obove one. The problem is the number of data points are different for these two data set. This is because one is barplot for frequency and the other is distribution. Anyone has expreience on this? Many thank

[R] how to make lien smoothly over barplot

Dear I try to make lines smoothly which are added on current barplot. I tried to use "smo" and it does not work. Has nayone have this experience? Many Thanks! Xin [[alternative HTML version deleted]] __ R-help@r-project.o

[R] how to add different type of lines (short dash, long dash) into current plot)

Dear I try to add three different lines (solid, short dash, long dash) in to current barplot. I saw there are types of lines ("p","b", ect). However, they are not what I request. Has nayone have this experience? Many Thanks! Xin [[alternative

[R] how to read data from table based a condition

ve a help? Thanks! Xin [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, mini

[R] multdrc error

it the reason to produce the error? Any clue about this error message? Thanks you in advance for your help, Xin Gong [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-h

[R] multdrc error---Error in mdrcOpt(opfct, startVec, optMethod, derFlag, constrained, warnVal

gly, is it the reason to produce the error? Thanks you in advance for your help, Xin Gong [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the post

[R] how to specify the location of tick mark on x axies

Dear: I want to plot barplot and let bar be in the middle of each x axis category. Do you have this experience? Many Thanks! Xin [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch

[R] anyone know how to plot histograma and line in the same graph

2.7168 2 2.3468 4 2.0283 2 1.7537 2 1.517 5 1.3127 1 1.1362 And I want to plot one of series by historgram and the other by line. such as the attachment. Thanks Xin __ R-help@r-project.org mailing list

[R] question about calling winbugs in R

your advice? Many Thanks! Xin [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

[R] help on R for loops

Hi, all. I need help on improving the efficiency of my R simulations. Below is a function that simulates a change point model. It first generates a sequence of three dimensional ARMA(1,1) observations, then calculates the one step ahead prediction errors, some statistic is calcualted and compare

[R] can optimize solve paired euqations?

s case? Thanks Xin [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal

[R] can R solve these paired equations

Dear: I have a paired equation below. Can I solve (x,y) using R. Thanks! Xin A=327.727 B=9517.336 p=0.114^10 (1-p)*y*(1-x)/x/(1-x^y)=A A(1+(1-x)*(1+y)/x-A))=B [[alternative HTML version deleted]] __ R-help@r-project.org

[R] can I solve this equation in R----29.040x+1=327.727x^(355.768x/(1-x))

-29.040x+1=327.727x^(355.768x/(1-x)) [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide c

[R] Syntax for Trim and Fill Method

Hi all, I am a graduate student doing some simulation study using R about publication bias in Meta Analysis. One method I need to implement is the Trim and Fill method. I know that there is a Trim and Fill "package" in S-Plus and Stata, so I am wondering wheter there is something in R that I c

[R] about Runtime Error

uot;. The memory count not be "read". Please let me know if there is anything I can do to solve the problem. Thank you very much. Xin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guid