Re: [R] why summary() does not work here???

2009-08-15 Thread Ronggui Huang
I cannot reproduce what you mentioned. One possible is that you use summary() without loading VGAM package. 2009/8/16 Hongwei Dong : > Hi, R users, >  I'm using the function "vglm" to estimate a multinomial logit model. Every > time I use "summary()" to ask for the coefficients and std error, I g

[R] why summary() does not work here???

2009-08-15 Thread Hongwei Dong
Hi, R users, I'm using the function "vglm" to estimate a multinomial logit model. Every time I use "summary()" to ask for the coefficients and std error, I got this: Length Class Mode 1 vglm S4 Anyone know what's wrong here? Thanks. Harry [[alternative HTML version delet

Re: [R] extracting text from wordnet using synonym()

2009-08-15 Thread Kelvin Lam
Thank you both! It works now. (Need to check out the parsing ) Gabor Grothendieck wrote: > > Or a variation: > >> unlist(lapply(parse(text = Str), eval)) > [1] "aid""assist" "assistance" "help" "aid" > [6] "assistance" "help" > > On Sat, Aug 15, 2009 at 5:55 PM, Henri

Re: [R] Estimate Discrete Choice Models with R

2009-08-15 Thread Ronggui Huang
glm can do probit regression, and glmer (lme4 package) and do mixed effect logit model. Ronggui 2009/8/16 Hongwei Dong : > Thanks for the replies. > What I'm really interest in are the functions that can do GEV, probit, and > mixed logit. I searched R, and I did find the function that can do the

Re: [R] Building 'R' from source for Windows.

2009-08-15 Thread Duncan Murdoch
rkevinbur...@charter.net wrote: Thank you for looking into this. It turns out the problem was "You are misinterpreting R_HOME. . . " I thought R_HOME was were I installed R not the directory where I was trying to compile the source. Once I moved the "extra" stuff that RTools.exe installed in w

[R] coefficient p-value in ridge regression

2009-08-15 Thread Maria
Hello. I'have a problem with RIDGE REGRESSION. I've used lm.ridge function to estimate coefficients of my model. Why in the summary of models not appears t value, Pr(>|t|) and significance stars? How I can calculate coefficient's p-value in ridge regression? Thanks! [[alternative HTML

Re: [R] How to plot 3-D surface graph from lmer mixed models?

2009-08-15 Thread willow1980
Hi, First of all, I thank Professor Murdoch and Dr Sarkar for providing important information. Now, this problem has been basically solved. The package "akima" helps to fultil the objective originally resorting to "expand.grid". Now, I can use contour or perspective 3-D plot. Hopefully, this exper

[R] question

2009-08-15 Thread amal saki
Hi, How to compute survival tree based on logrank statistics from R package? [[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-pr

[R] System is computationally singular and scale of covariates

2009-08-15 Thread Stephan Lindner
Dear all, I'm running a self-written numerical optimization routine (hazard model) which includes computing the inverse of the outer product of the score. I have been getting the above error message ("System is computationally singular"), and after some tweaking, I realized that these variables h

Re: [R] Building 'R' from source for Windows.

2009-08-15 Thread rkevinburton
Thank you for looking into this. It turns out the problem was "You are misinterpreting R_HOME. . . " I thought R_HOME was were I installed R not the directory where I was trying to compile the source. Once I moved the "extra" stuff that RTools.exe installed in what I thought was the R installat

Re: [R] extracting text from wordnet using synonym()

2009-08-15 Thread Gabor Grothendieck
Or a variation: > unlist(lapply(parse(text = Str), eval)) [1] "aid""assist" "assistance" "help" "aid" [6] "assistance" "help" On Sat, Aug 15, 2009 at 5:55 PM, Henrique Dallazuanna wrote: > Try this: > > Str <- c("c(\"aid\", \"assist\", \"assistance\", \"help\")", >         "c(\"

Re: [R] extracting text from wordnet using synonym()

2009-08-15 Thread Henrique Dallazuanna
Try this: Str <- c("c(\"aid\", \"assist\", \"assistance\", \"help\")", "c(\"aid\", \"assistance\", \"help\")") unlist(sapply(Str, function(x)dget(textConnection(x)), USE.NAMES = FALSE)) On Sat, Aug 15, 2009 at 6:33 PM, Kelvin Lam wrote: > > Dear group, > > There are postings on extract

[R] extracting text from wordnet using synonym()

2009-08-15 Thread Kelvin Lam
Dear group, There are postings on extracting partial string using sub()/gsub()/gsubfn() but I still couldn't figure this one out. Supposed I am pulling the synonyms of "help" using synonyms("help") from wordnet and get the followings: > synonyms("help") [1] "c(\"aid\", \"assist\", \"assistance\

Re: [R] Filling in a uniform dist plot

2009-08-15 Thread David Kaplan
Yep. That's it. Thanks. David Peter Dalgaard wrote: David Kaplan wrote: Hi all, Apologies in advance if this is really too simple, but I have given a good shot at this. I'm generating a standard uniform distibution y <- dunif(x) plot(x,y,type="l") This will draw a horizontal line

Re: [R] Error in running RWeka Clusteres

2009-08-15 Thread Farid Seifi
Hi all, The problem is solved uninstalling the R and restarting the system and installing it again. I don't know what caused the error but reinstalling the R was not enough and the system must be restarted. Thanks Farid -Original Message- From: "Farid Seifi" To: R-help@r-project.or

Re: [R] Filling in a uniform dist plot

2009-08-15 Thread Peter Dalgaard
David Kaplan wrote: Hi all, Apologies in advance if this is really too simple, but I have given a good shot at this. I'm generating a standard uniform distibution y <- dunif(x) plot(x,y,type="l") This will draw a horizontal line at y = 1. I want to fill the area from x=(0,1) with a color.

[R] Filling in a uniform dist plot

2009-08-15 Thread David Kaplan
Hi all, Apologies in advance if this is really too simple, but I have given a good shot at this. I'm generating a standard uniform distibution y <- dunif(x) plot(x,y,type="l") This will draw a horizontal line at y = 1. I want to fill the area from x=(0,1) with a color. I've tried the poly

Re: [R] Estimate Discrete Choice Models with R

2009-08-15 Thread Hongwei Dong
Thanks for the replies. What I'm really interest in are the functions that can do GEV, probit, and mixed logit. I searched R, and I did find the function that can do the regular multinomial logit model, but I did not see the function that can do GEV, probit, and mixed logit. Any further advices? Th

Re: [R] Solutions of equation systems

2009-08-15 Thread Ravi Varadhan
Since you have an under-determined system of linear equations, you have infinitely many solutions. One reasonable solution is the "minimum-norm" solution. You can obtain this from singular-value decomposition as follows: # Minimum norm solution `x': A %*% x = b, such that ||x|| is minumum d

[R] Error in running RWeka Clusteres

2009-08-15 Thread Farid Seifi
Hi, I have a question about using RWeka Clusterers.If you could supply answer or insight, I would really appreciate it. When I run a simple code which uses a clusterer from RWeka I get an error. the sample codes and errors are mentioned below Code: library(RWeka) Cobweb(iris[,-5],control=NUL

Re: [R] How to use R to perform prediction based on history data

2009-08-15 Thread Daniel Malter
Please have a look at the posting guide. For your problem of loading the data, we do not know what you have done and, therefore, cannot even try to guess what the reason for the error message may be. So at least we need information what you did in R (the code). Second, the posting guide generally r

Re: [R] aov p values for levels of a factor.

2009-08-15 Thread Dieter Menne
John Sorkin wrote: > > I am running an balanced aov, in which Time has 4 levels, Group has 2 > levels: > > fit2<-aov(Post ~ as.factor(Time) > +as.factor(Group)+as.factor(Group*Time)+Error(SS), alldata) > > I would like to get a p value for each level of time, each level of > Group and each le

Re: [R] re ading transposed data in read.csv

2009-08-15 Thread Dieter Menne
> tdf <- read.table(textConnection("ID samp1 samp2 samp3 samp4 + G1 2332 12 87 + G2 8545 49 76 + G3 1246 39 28 + G4 7326 18 13 + "), header=T) > library(reshape) > recast(tdf, ...~ID) Using ID as id

Re: [R] reading in mdb and outputting to sql (GIS application)

2009-08-15 Thread Frank E Harrell Jr
stephen sefick wrote: tables=TRUE does return a list of tables with no error. The problem is that I can't get into the database because I am on a macintosh. I would send you the database, but it is for a military base and I don't believe that I am allowed. I'll see if I can't figure it out. T

Re: [R] Solutions of equation systems

2009-08-15 Thread Berend Hasselman
Moreno Mancosu wrote: > > Dear Greg, > > I tried to use function "solve", but I have some problems with it. > Let's try with a simple equation: > > x + 3y -2z = 5 > 3x + 5y + 6z =7 > > > a<-matrix(c(1,3,-2,3,5,6),nrow=2,ncol=3) > > a > [,1] [,2] [,3] > [1,]1 -25 > [2,]3

[R] VAR (pckg: vars) and memory problem

2009-08-15 Thread herrdittmann
Hi all, When I tried to estimate a VAR (package vars) of a rather large dataset with 5 lags: > dim(trial.var) [1] 20388 2 I ran into memory troubles: > summary(VAR(trial.var, type="none", p=5)) Error: cannot allocate vector of size 3.1 Gb In addition: Warning messages: 1: In diag(r

Re: [R] reading in mdb and outputting to sql (GIS application)

2009-08-15 Thread stephen sefick
tables=TRUE does return a list of tables with no error. The problem is that I can't get into the database because I am on a macintosh. I would send you the database, but it is for a military base and I don't believe that I am allowed. I'll see if I can't figure it out. Thanks for all of the hel

Re: [R] reading in mdb and outputting to sql (GIS application)

2009-08-15 Thread Frank E Harrell Jr
stephen sefick wrote: library(Hmisc) # Read all tables in the Microsoft Access database blabla.mdb # perhaps better change to Locale in which the .mdb file was # generated to avoid that accented table column names are # getting lost.file first Sys.setenv(PATH=paste(Sys.getenv('PATH'),'/sw/b

Re: [R] Solutions of equation systems

2009-08-15 Thread Moreno Mancosu
Greg Snow wrote: ?solve Dear Greg, I tried to use function "solve", but I have some problems with it. Let's try with a simple equation: x + 3y -2z = 5 3x + 5y + 6z =7 > a<-matrix(c(1,3,-2,3,5,6),nrow=2,ncol=3) > a [,1] [,2] [,3] [1,]1 -25 [2,]336 > b<-matrix(c(5,

Re: [R] re ading transposed data in read.csv

2009-08-15 Thread Vitalie S.
On Sat, 15 Aug 2009 09:20:04 +0200, Dieter Menne wrote: sandsky wrote: I have a data set (csv); e.g., ID samp1 samp2 samp3 samp4 G1 2332 12 87 G2 8545 49 76 G3 1246 39 28 G4 7326 18 13 and

Re: [R] re ading transposed data in read.csv

2009-08-15 Thread Dieter Menne
sandsky wrote: > > I have a data set (csv); e.g., > > ID samp1 samp2 samp3 samp4 > G1 2332 12 87 > G2 8545 49 76 > G3 1246 39 28 > G4 7326 18 13 > > and read it: > > data1<-(read.csv("Datafolder

Re: [R] Building 'R' from source for Windows.

2009-08-15 Thread Peter Dalgaard
rkevinbur...@charter.net wrote: I know I am going to catch alot of comments for this question but I am really stuck. If there is some written documentation that I have missed please redirect me. I want to build 'R' from source on a Windows Platform. The main reasons are that I want to check o

[R] How to use R to perform prediction based on history data

2009-08-15 Thread Jack Su
Say I have a csv file, each row contains several fields, one of them are whether the row is success. In history data, I have all the fields including the result of whether it is success. In future data, I only have fields without the result. For example: history data: Field1 Field2 Field3 F

[R] re ading transposed data in read.csv

2009-08-15 Thread sandsky
I have a data set (csv); e.g., ID samp1 samp2 samp3 samp4 G1 2332 12 87 G2 8545 49 76 G3 1246 39 28 G4 7326 18 13 and read it: data1<-(read.csv("Datafolder/rawdata.csv",header=T)) It is fine with