[R] Run-time error "521" in SciViews

2006-09-14 Thread Wuming Gong
Dear list, I use SciViews R Console 0.8.9 with R-2.2.1 under Windows XP SP2, and it works very well for most of time. However, sometimes, when commands were executed by clicking F5, the error "Run-time error 521, Can't open clipboard" pop out. After choosing "Yes", the both R console and SciView

[R] Error when loading odesolve

2006-08-04 Thread Wuming Gong
Dear list, I installed odesolve package (0.5-15) in R 2.3.1 in a Solaris server (Generic_118558-11 sun4u sparc SUNW,Sun-Blade-1000). The installing progress completed without errors, though several warnings like "Warning: Option -fPIC passed to ld, if ld is invoked, ignored otherwise" were output

[R] Polynomial kernel in SVM in e1071 package

2006-07-07 Thread Wuming Gong
Dear list, In some places (for example, http://en.wikipedia.org/wiki/Support_vector_machine) , the polynomail kernel in SVM is written as (u'*v + 1)^d, while in the document of svm() in e1071 package, the polynomial kernel is written as (gamma*u'*v + coef0)^d. I am a little confused here: When d

Re: [R] R-2.3.0 make error

2006-05-02 Thread Wuming Gong
2006 08:51:21 +0200, Peter Dalgaard <[EMAIL PROTECTED]> wrote: > "Wuming Gong" <[EMAIL PROTECTED]> writes: > > > Dear Peter, > > > > $ uname -a > > Linux bl3 2.6.15-1.1833_FC4smp #1 SMP Wed Mar 1 23:55:52 EST 2006 > > x86_64 x86_64 x86_64 GNU/L

Re: [R] R-2.3.0 make error

2006-05-01 Thread Wuming Gong
Dear Peter, $ uname -a Linux bl3 2.6.15-1.1833_FC4smp #1 SMP Wed Mar 1 23:55:52 EST 2006 x86_64 x86_64 x86_64 GNU/Linux When ./configure, I did not claim any options... Wuming On 01 May 2006 19:30:48 +0200, Peter Dalgaard <[EMAIL PROTECTED]> wrote: > "Wuming Gong" <[EM

[R] R-2.3.0 make error

2006-05-01 Thread Wuming Gong
Dear list, When compiling the R-2.3.0 on FC4 x86_64, I got the following errors: make[3]: Entering directory `/project/scratch3/ligroup/wuming/src/R-2.3.0/src/main' gcc -Wl,--export-dynamic -L/usr/local/lib64 -o R.bin Rmain.o CConverters.o CommandLineArgs.o Rdynload.o Renviron.o RNG.o apply.o ari

[R] Sub-vector

2006-03-29 Thread Wuming Gong
Dear list, Given a vector of logical values, say >a <- c(TRUE,TRUE,TRUE,FALSE,FALSE,FALSE,FALSE,TRUE,TRUE,TRUE) Are there any R functions that can tell whether there are two or more "TRUE" in a row in this vector? Thanks, Wuming __ R-help@stat.math.

[R] R package for computing state path using Viterbi algorithm

2006-03-22 Thread Wuming Gong
Dear list, This question is about Hidden Markov Model. Given a transition matrix, an emission matrix and a sequence of observed symbols (actually, nucleotide sequences, A, T, C and G), I hope to predict the sequence of state by Viterbi algorithm. I searched R repository for related packages. ms

Re: [R] Can anyone please tell me how to strip the white spaces from a character vector?

2005-10-25 Thread Wuming Gong
Try trim() in gdata package. Wuming On 10/25/05, roger bos <[EMAIL PROTECTED]> wrote: > for example: > > a$tic[1:10] > [1] "AIR " "ABCB " "ABXA " "ACMR " "ADCT " "ADEX " > [7] "ABM " "AFCE " "AG " "ATG " > Can anyone please tell me how to strip the white spaces from a$tic? > Thanks, > Roger >

Re: [R] R-help Digest, Vol 31, Issue 9

2005-09-10 Thread Wuming Gong
lt;=.05||coe[4]<=.05||coe[5]<=.05) { > cat("phenotype is = ",x[j] , "\n") > cat("snp marker is = ",x[i] , "\n") > sign<-call("fstatistic",zz) > call("p-value",zz) > >#print(coe) >print(zz

Re: [R] R-help Digest, Vol 31, Issue 9

2005-09-10 Thread Wuming Gong
?summary.lm and check the Value section. Wuming On 9/10/05, Ping Yao <[EMAIL PROTECTED]> wrote: > Hi: > I use lm (linear model) to analyze 47 variables , 8 responses > So I use loop to finish it . > I want the program to show the results that P-value is less than 0.05. > How can I cite the P-valu

Re: [R] "Too long to display" problem

2005-09-08 Thread Wuming Gong
Dear list, Please ignore this thread - the "TOO LONG TO DISPLAY" is brought by another tool when parsing data sets. Sorry for this ... Wuming On 9/8/05, Wuming Gong <[EMAIL PROTECTED]> wrote: > Dear list, > > I used read.xls in gdata package to read a worksheet

[R] "Too long to display" problem

2005-09-08 Thread Wuming Gong
Dear list, I used read.xls in gdata package to read a worksheet in which certain field contains very long character strings (nucleotides sequence, nchar > 10,000). Then, the values in these fields are automatically converted to "TOO LONG TO DISPLAY". How can I get those original characters instea

[R] Calculating Goodman-Kurskal's gamma using delta method

2005-09-01 Thread Wuming Gong
Dear list, I have a problem on calculating the standard error of Goodman-Kurskal's gamma using delta method. I exactly follow the method and forumla described in Problem 3.27 of Alan Agresti's Categorical Data Analysis (2nd edition). The data I used is also from the job satisfaction vs. income ex

Re: [R] Gamma for ordinal trends

2005-08-31 Thread Wuming Gong
> > > On 08/31/05 17:46, Wuming Gong wrote: > > > Dear list, > > > > > > Are there any functions for calculating gamma (and its standard > > > error), which measures the association of ordinal factors in I x J > > > contingency table. I did a R

[R] Gamma for ordinal trends

2005-08-31 Thread Wuming Gong
Dear list, Are there any functions for calculating gamma (and its standard error), which measures the association of ordinal factors in I x J contingency table. I did a RSiteSearch but did not find any clues... Thanks, Wuming __ R-help@stat.math.ethz

Re: [R] NA when using read.csv

2005-08-01 Thread Wuming Gong
Hi Haibo, Are there any character "#" in the cells? By default, the things at the right side of the character "#" will be masked as comments. Wuming On 8/2/05, Gabor Grothendieck <[EMAIL PROTECTED]> wrote: > On 8/1/05, Haibo Huang <[EMAIL PROTECTED]> wrote: > > Hi, when I used: > > > > Dist=read

[R] How to find the significant differences among interactions in logit model?

2005-06-24 Thread Wuming Gong
Hi, I have a question about interpret the results from logistic regression model. I used a dataset from the book Categorical Data Analysis (2nd Edition) by Alan Agresti. > summary(crabs) color spine width satell weightpsat 2:12 1: 37 Min. :21.0

Re: [R] how 'stepAIC' selects?

2005-06-19 Thread Wuming Gong
Hi Veronika, StepAIC does not do any test. But you may do a likelihood ratio test on the difference between deviance of two models to judge which model to use, simply like this, > 1 - pchisq(deviance.i - deviance.j, df.i - df.j) And then check whether the returned p value is small enough... Wu

Re: [R] Logistic regression with more than two choices

2005-06-14 Thread Wuming Gong
Hi Koskinen For response variables with multiple categories, you may try polr() in MASS package, which implement a proportional odds model. And you may search the R archives, several threads discussed this problem before... Wuming On 6/15/05, Ville Koskinen <[EMAIL PROTECTED]> wrote: > Dear all

Re: [R] logistic regressioin - course ornotes

2005-06-10 Thread Wuming Gong
Hi Stephen, I think the tutorial written by Brett Presnell (http://web.stat.ufl.edu/~presnell/Teaching/sta4504-2000sp/R/) is a good start point for categorical data analysis using R. Wuming On 6/10/05, Stephen Choularton <[EMAIL PROTECTED]> wrote: > Hi > > I am using R for logistic regression a

Re: [R] how to plot a fitted smooth line on histograms?

2005-05-22 Thread Wuming Gong
You may check the chapter 8 Probability distributions in An Introduction to R. Wuming On 5/22/05, Hu Chen <[EMAIL PROTECTED]> wrote: > for example, > >x <- read.table(); here x is a vector containing my data to be > >analyzed. > >hist(x,plot=TRUE,breaks=200) > > then I got a histogram gra

Re: [R] covariance analysis by using R

2005-05-20 Thread Wuming Gong
You may fit the model using lm() directly - R will set up a coding for qualitative predictor automatically (taking experiments as qualitative predictor). HTH Wuming On 5/18/05, 孟欣 <[EMAIL PROTECTED]> wrote: > Hello sir: > Here's a question on covariance analysis which needs your help. > There'r

Re: [R] simple question, i hope

2005-05-17 Thread Wuming Gong
See ?write.table Wuming On 5/17/05, BJ <[EMAIL PROTECTED]> wrote: > How do you output a list to a text file without the extra line numbers > and stuff? > > I have a list b, and tried > > zz<-textConnection("captest.txt","w") > sink(zz) > b > sink() > close(zz) > > but that isnt what i want, b

[R] How can I estimate parameters of probability distributions?

2004-12-09 Thread Wuming Gong
Hi list, I have a group of data. It looks like they follow a exponential distribution. In R, how can I esimate lamda, that is the rate in pexp, of the distribution and can I use Kolmogorov-Smirnov for hypothesis testing in such a situation? I have read the "8.2 Examing the distribution of a set o

[R] How can I generate a random network with given cluster coefficient using R?

2004-11-27 Thread Wuming Gong
Hi list, I have read the mannual of package graph, but it looks that three functions randomEGraph, randomGraph and randomNodeGraph cannot generate a random network with given cluster coefficient. Is there any exisiting R package for such kind job and if no, how can I perform this job with existi

[R] How can I generate a random network with given cluster coefficient using R?

2004-11-27 Thread Wuming Gong
Hi list, I have read the mannual of package graph, but it looks that three functions randomEGraph, randomGraph and randomNodeGraph cannot generate a random network with given cluster coefficient. Is there any exisiting R package for such kind job and if no, how can I perform this job with existi