[R] fractal package

2012-07-29 Thread mary Joy Regidor
Greetings of Peace! I am new in R software. I want to use the correlation dimension corrDim in computing a map. I have installed the packages(fractal) in R. Now my problem is this, when I tried your example in the following site http://rgm2.lab.nig.ac.jp/RGM2/func.php?rd_id=fractal:corrDim  > l

[R] distance matrix and hclustering

2012-07-29 Thread eliza botto
Dear R Users,i am very new to R. I want your help on an issue regarding distance matrix and cluster analysis i had discharge data of 4 rivers(a,b,c,d) in 4 vectors each having 364 values > dput(qmu)structure(list(a = c(0.26, 0.25, 0.25, 0.25, 0.24, 0.23, 0.22, 0.21, > 0.21, 0.21, 0.2, 0.19, 0.19

Re: [R] Extracting results from the VAR output

2012-07-29 Thread Rmillan
Hi guys, Thank you very much :) It seems that both your suggestions works. I wasn't aware of the function [[argument]], it makes the programming procedure much easier! :) Best regards, Emil -- View this message in context: http://r.789695.n4.nabble.com/Extracting-results-from-the-VAR-outpu

Re: [R] extended R

2012-07-29 Thread Ben Bolker
hony gmail.com> writes: > > I need to create a software and use R engine in my software. but I don't know > how to use R code in my software same www.revolutionanalytics.com. > Rapid Miner has 2 section. one of them is using for application and the > other is used for researcher for extended it.

[R] rJava loading error on hp-ux

2012-07-29 Thread 안대영
Hi all, When I installed rJava from source, rJava installation was failed with below messges. $ R CMD INSTALL rJava_0.9-4.tar.gz ... ... ... gmake[1]: Leaving directory `/tmp/RtmpNwiCae/R.INSTALL4ebb76a3400d/rJava/jri' installing via 'install.libs.R' to /stat/R-2.15.1/lib/R/libr

[R] curve comparison

2012-07-29 Thread Luis Fernando García Hernández
Dear R users, I have seven regression lines I´d like to compare, in order to find out if these are significatively different. The main problem is that these are curves, non normal, non homogeneous data, I´ve tried to linearize them but it has not worked. So I´d like to know if you know any comman

[R] R- Help (looping)

2012-07-29 Thread Wellington G. Silva
Hi, I'm Wellington from Brazil and I have the following issue: I've been working on a project a for a while, and I'm having trouble in using the loop (for) I need to read a column (c1), and for each value of this column, I need to check if it's within the control limits So, I was t

[R] Zoo panel function

2012-07-29 Thread Gene Leynes
I would really like some help with understanding the panel function, in zoo. Thank you. R 15.1 and zoo 1.7-7. library(zoo) x = seq(0,3*pi,length.out=100) y = sin(x) zobj = zoo(y, x) ### ## EXAMPLE 1 - GLOBAL ARGUMENT ## This panel functio

Re: [R] Hints for the best package

2012-07-29 Thread MP
Thanks a million, Michael. I really appreciated it. Marco -Original Message- From: R. Michael Weylandt [mailto:michael.weyla...@gmail.com] Sent: Sunday, July 29, 2012 5:03 PM To: MP Cc: r-help@r-project.org Subject: Re: [R] Hints for the best package It sounds like you want multiple lo

[R] why querying Sybase IQ using RODBC returns error ?

2012-07-29 Thread Alok Jadhav
Hello, I am not able to query Sybase IQ database from R using RODBC package. Check the below code R> sybaseiq.query <- function (sql, ..., as.is = FALSE) { connstr <- "Driver={Sybase IQ};UID=ajadhav2;PWD=*;ServerName=PSGSESHR01A_IQ;CommLinks=SharedMemory,TCPIP{host=psgseshr01_iq.sg.cs

[R] problem plotting in a grid

2012-07-29 Thread Jonathan
Hi all, I'm trying to generate a grid of four plots. The first 2 appear just fine, but the final 2 will not appear in the grid, instead overwriting the first two.Any ideas on how to get them all in the same window would be greatly appreciated. Cheers, Jonathan library(fields) par(mfrow=

[R] Z score in gplots

2012-07-29 Thread Khush gg
Hi, Can anyone tell me how to set Z-score according to my own requirement as the below code is taking as per the file entries. Any help would be appreciable. library(gplots) x=read.table("final.txt", header=TRUE) mat=data.matrix(x) heatmap.2(mat, col=colorRampPalette(c("green","white","red"))(256

Re: [R] extended R

2012-07-29 Thread Steve Lianoglou
Hi, On Sun, Jul 29, 2012 at 5:58 AM, hony wrote: > I need to create a software and use R engine in my software. but I don't know > how to use R code in my software same www.revolutionanalytics.com. If your software is written in C++ perhaps the RInside package will be helpful: http://cran.r-pro

Re: [R] Understanding the intercept value in a multiple linear regression with categorical values

2012-07-29 Thread Joao Azevedo
Hi! You're right. I was misinterpreting the way the coefficients were calculated. Reading about the method of least squares helped me in clarifying some of my doubts. Thanks for your tips! -- Joao. On Fri, Jul 27, 2012 at 2:36 PM, Jean V Adams wrote: > Joao, > > Your intuition is correct, the

Re: [R] Reverse order of x-axis

2012-07-29 Thread David L Carlson
Show us what you actually did. Rui's example does not generate a problem with a log axis but it is an index plot. This also does not generate any problems: > set.seed(42) > x <- runif(10)*100 > y <- runif(10)*100 > plot(x, y, xlim=c(100, 1), log="x") You wouldn't have used xlim=c(100, 0) and the

Re: [R] Nearest Neighbors

2012-07-29 Thread David L Carlson
Actually if you are looking for neighbors of neighbors you only need the nearest neighbor for each point. The problem is that the nearest neighbor of the nearest neighbor of point 1 is often point 1. Did you want the nearest neighbor not counting any point twice? Sounds more like a traveling salesm

Re: [R] Reverse order of x-axis

2012-07-29 Thread Rolf Turner
(1) Please show context; r-help is *not* nabble. (2) You might like to investigate the revaxis() function from the "plotrix" package. cheers, Rolf Turner On 30/07/12 07:41, sappy wrote: Thanks, for the moment it works. But is there a possibility for a logarithmic x-Axis? When i

Re: [R] Simple x,y Plot

2012-07-29 Thread Rui Barradas
Hello, Inline Em 29-07-2012 17:41, arun escreveu: Hi, As there is not much information like str() or the data, it is difficult to comment. But, try this (considering your x and y are factors): x<-c("1","4","7","9","10","15") y<-c("3","5","8","6","5","14") dat1<-data.frame(x,y) str(dat

Re: [R] Hints for the best package

2012-07-29 Thread R. Michael Weylandt
It sounds like you want multiple logistic regression: I think most people find the mlogit package quite helpful there. See also: http://www.ats.ucla.edu/stat/r/dae/mlogit.htm Michael On Sun, Jul 29, 2012 at 9:02 AM, MP wrote: > I am new to this field and I am currently studying the most basic pr

Re: [R] fitting Markov Switching Model

2012-07-29 Thread Ingmar Visser
Searching markov switch on the R website gives a list of useful hints and packages: https://www.google.com/search?q=markov+switch&domains=r-project.org&sitesearch=r-project.org&btnG=Google+Search If you know the cut-off between states no complicated modeling is necessary, using something like thi

Re: [R] Reverse order of x-axis

2012-07-29 Thread Rui Barradas
Hello, Just reverse the xlim values, it defaults to c(min, max), use c(max, min). x <- 1:10 plot(x) plot(x, xlim=c(10, 1)) Hope this helps, Rui Barradas Em 29-07-2012 18:23, sappy escreveu: Hello, does anybody know, how to reverse the elements of the x-axis. At now my x-axis for example goe

Re: [R] Reverse order of x-axis

2012-07-29 Thread sappy
Thanks, for the moment it works. But is there a possibility for a logarithmic x-Axis? When i use log='x' i get only a vertical line. Greetings -- View this message in context: http://r.789695.n4.nabble.com/Reverse-order-of-x-axis-tp4638275p4638282.html Sent from the R help mailing list archiv

Re: [R] readRDS, In as.double.xts(fishReport$count) : NAs introduced by coercion

2012-07-29 Thread arun
HI Yolande, No problem.  I tried to get a solution to your old post yesterday after converting xts object to dataframe.  But, still I am getting the error message that the lengths differ.  I don't know how to fix that.  A.K. From: Yolande Tra To: arun

Re: [R] Simple x,y Plot

2012-07-29 Thread sappy
Thanks a lot. It Works. -- View this message in context: http://r.789695.n4.nabble.com/Simple-x-y-Plot-tp4638262p4638274.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/l

Re: [R] Simple x,y Plot

2012-07-29 Thread arun
Hi, As there is not much information like str() or the data, it is difficult to comment. But, try this (considering your x and y are factors): x<-c("1","4","7","9","10","15")  y<-c("3","5","8","6","5","14")  dat1<-data.frame(x,y)  str(dat1) 'data.frame':    6 obs. of  2 variables:  $ x: Factor

[R] Reverse order of x-axis

2012-07-29 Thread sappy
Hello, does anybody know, how to reverse the elements of the x-axis. At now my x-axis for example goes from 0, 0.5, 1, 1.5, ... to 99.5, 100.0 Is there a way to change the x-axis, like this: 100.0, 99.5, 99.0, ..., 0.5, 0. The 100.0 should start at left x-axis and the 0 shold be on the right. Gre

Re: [R] Nearest Neighbors

2012-07-29 Thread Petr Savicky
On Tue, Jul 24, 2012 at 09:26:49AM -0700, olemissrebs1123 wrote: > I was wondering if there is a way in R to find k nearest neighbors of various > orders, say order 2, 3, or 4. In otherwords neighbors of neighbors of > neighbors. You get the idea. I know that I can use knearneigh(matrix.data, > k)

Re: [R] taylor expansions with real vectors

2012-07-29 Thread Petr Savicky
On Sun, Jul 22, 2012 at 05:34:08PM -0700, bilelsan wrote: > Hi, > Thanks a lot for answer. It is what I mean. > But the code does not seem to work ( Hi. I am sorry for a late reply. I was on vacations one week. Can you specify the problem? I get e1 <- 2 e2 <- 3 k <- 3 f <- 0 for (r

Re: [R] readRDS, In as.double.xts(fishReport$count) : NAs introduced by coercion

2012-07-29 Thread Yolande Tra
Thanks. It was very helpful. Y On Sun, Jul 29, 2012 at 2:45 AM, arun wrote: > Hello, > > I tried to do the same thing that you have done. Without the timestamps > as rownames (xts object), there were no errors. > > dat1<-read.table(text=" > dive_id datetime speciescount si

Re: [R] problems with R startup in mac

2012-07-29 Thread Prof Brian Ripley
Please post on R-sig-mac (see the posting guide, footer of this email). They will need to know how you are running R (command-line or R.app) and your language settings. On 29/07/2012 09:32, 肖文迪 wrote: Dear sir/madam, I encountered some problems when running R on my computer. To be honest, I

Re: [R] Simple x,y Plot

2012-07-29 Thread Peter Ehlers
On 2012-07-29 08:18, sappy wrote: Dear Community, i have read in a table with two column from a *txt-File. What i get is a data-frame wit two column. After that i allocate this 2 column to x and y. Plotting this, i obtain a picture that looks like a bar plot (something like a grid). Is that beca

Re: [R] Simple x,y Plot

2012-07-29 Thread mlell08
On 29.07.2012 17:18, sappy wrote: > Dear Community, > > i have read in a table with two column from a *txt-File. What i get is a > data-frame wit two column. > After that i allocate this 2 column to x and y. Plotting this, i obtain a > picture that looks like a bar plot (something like a grid). >

Re: [R] Error in for-loop

2012-07-29 Thread jim holtman
Try this: You had some 'indexing' problems. You were accessing the columns of the dataframe incorrectly. You are also accessing 'list' so you might want to review the Intro to R on indexing: > datos.mx1 <- > read.table('http://r.789695.n4.nabble.com/file/n4638257/PETmx1.csv',head=TRUE,sep=';')

Re: [R] problems with R startup in mac

2012-07-29 Thread peter dalgaard
On Jul 29, 2012, at 10:32 , 肖文迪 wrote: > Dear sir/madam, > > I encountered some problems when running R on my computer. To be honest, I > am only a beginner of computer programming and I have no idea how to deal > with it. I have tried to change the font in "system preference" but it did > not w

[R] fitting particular distributions

2012-07-29 Thread Jonathan Hughes
Hi there, I got a question that is both about stats and R. Imagine two alternative stochastic processes: a. Each subject has an exponential distribution of lifetime, with a parameter lambda; If I simulate 100 such subjects and rank their lifetimes, they'd look like this: >plot(sort(rex

[R] MetaCost (and Random Forest)

2012-07-29 Thread Annelen Brunner
Hello! Is there an implementation of the MetaCost algorithm in R? I found a project here: https://r-forge.r-project.org/scm/viewvc.php/src/mlrCost/metacost.r?root=mlr&view=log But I'm not sure about its status - it seems unfinished to me. My goal is to use MetaCost for cost-sensitive learning of

[R] Simple x,y Plot

2012-07-29 Thread sappy
Dear Community, i have read in a table with two column from a *txt-File. What i get is a data-frame wit two column. After that i allocate this 2 column to x and y. Plotting this, i obtain a picture that looks like a bar plot (something like a grid). Is that because x and y is in factor-formatted?

[R] extended R

2012-07-29 Thread hony
I need to create a software and use R engine in my software. but I don't know how to use R code in my software same www.revolutionanalytics.com. Rapid Miner has 2 section. one of them is using for application and the other is used for researcher for extended it. -- View this message in context:

[R] problems with R startup in mac

2012-07-29 Thread 肖文迪
Dear sir/madam, I encountered some problems when running R on my computer. To be honest, I am only a beginner of computer programming and I have no idea how to deal with it. I have tried to change the font in "system preference" but it did not work. Can you tell me how to get it right? The syste

[R] Hints for the best package

2012-07-29 Thread MP
I am new to this field and I am currently studying the most basic problem of statistics, I think. I have an indicator that is 0 or 1 and many categorical and continuous variables. I have to create a predictive model for the indicator. Could anyone provide me a direction in terms of packages that

[R] Choice of model and specification of error terms.

2012-07-29 Thread David Labonte
Hi there! My question is not necessarily specific to R, bit I still hope that I can find help here. I appreciate any suggestions; how to do/improve my analysis, what stuff to read more on etc... I am interested in a dependent variable A. I measured A of 10 individuals at 7 different levels of

[R] Error in for-loop

2012-07-29 Thread Domi
Hello erverybody, I have a problem with my second for-loop. 1. First i read the tables. datos.mx1 <- read.table('PETmx1.csv',head=TRUE,sep=';') datos.min <- read.table('PETmin.csv',head=TRUE,sep=';') http://r.789695.n4.nabble.com/file/n4638257/PETmx1.csv PETmx1.csv http://r.789695.n4.nabble.

Re: [R] quantreg Wald-Test

2012-07-29 Thread Roger Koenker
On Jul 28, 2012, at 8:55 AM, stefan23 wrote: > Dear all, > I know that my question is somewhat special but I tried several times to > solve the problems on my own but I am unfortunately not able to compute the > following test statistic using the quantreg package. Well, here we go, I > appreciate

Re: [R] using save() to work with objects that exceed memory capacity

2012-07-29 Thread R. Michael Weylandt
On Sat, Jul 28, 2012 at 10:48 AM, David Romano wrote: > Context: I'm relatively new to R and am working with very large datasets. > > General problem: If working on a dataset requires that I produce more than > two objects of roughly the size of the dataset, R quickly uses up its > available mem

Re: [R] Mixed-model with paired design & covariates

2012-07-29 Thread Bert Gunter
Wrong list! Post on r-sig-mixed-models instead of here. However, this is really not an R question. It appears that you are looking for remote statistical consulting, which I consider rather hazardous. Sort of like having a virtual thesis advisor. I believe you would be better off talking with your

Re: [R] Appending the Column names

2012-07-29 Thread John Kane
You cannot have two sets of column names in a data.frame You can do this rbind(X, Y) Summary G Y R 1 Acc 12 12 13 2 Bcc 11 14 15 3 Ccc 13 15 16 4 Acc 10 11 12 5 Bcc 13 12 11 6 Ccc 11 16 20 John Kane Kingston ON Canada > -Original Message- > From: saileshch

Re: [R] Coloring Counties in a State Map

2012-07-29 Thread John Kane
I think that first you need to tell us :what kind of map you are using, what mappig package(s) you are using Why you want the colours and what data set you are drawing on to get these colours (Note it possible provide some sample data (See ?dput as a way to do this And provide any code that you

Re: [R] Return od functions

2012-07-29 Thread Rui Barradas
Hello, return( list(y, x, z, dy, dx) ) # or return( list(y=y, x=x, z=z, dy=dy, dx=dx) ) That's ONE variable returned, a list. In the first case you access the list members by number, ret <- dgp.par(...) ret[[ 1 ]] # 'y' value In the second case, by number or name. ret[[ 4 ]] ret[[ "dy" ]]

Re: [R] Problem with a regression - Dataset Workinghours

2012-07-29 Thread peter dalgaard
On Jul 28, 2012, at 17:37 , Giorgio Monti wrote: > I'm a student. I'm working on a research using the statistical program "R > 2.15.1". > Here's my problem: how i can do a regression considering only values over a > certain limit? > For example, considering the dataset "Workinghour" of the "Ecdat

[R] Return od functions

2012-07-29 Thread wangchi3
Hi! I have some questions about R function. I try to write a function for multi-returns. The function code is as attachment. dgp.par<-function(ai, bi, t, n) { t0<-t+20 y0<-matrix(0, nr=t0, nc=n) y0[1,]<-ai/(1-bi) for(tt in 2:t0) { y0[tt,]<-ai+bi*y0[tt-1,]+rnorm(n, 0, 1) } y<-y0[21:t0,] x<-y0[2

Re: [R] lm without intercept

2012-07-29 Thread peter dalgaard
On Feb 18, 2011, at 14:20 , Jan wrote: > One of the references you googled suggests that intercepts should never > be omitted. Is this true even if I know that the physical reality behind > the numbers suggests an intercept of zero? No. That'll be a piece of pragmatic advice caused by the experie

Re: [R] lm without intercept

2012-07-29 Thread peter dalgaard
On Feb 18, 2011, at 14:20 , Jan wrote: > Hello Achim, > >> Not quite. Consult your statistics textbook for the correct interpretation >> of p-values. Under the null hypothesis of a true intercept of zero, it is >> very likely to observe an intercept as large as 13.52 or larger. > thank you for

Re: [R] readRDS, In as.double.xts(fishReport$count) : NAs introduced by coercion

2012-07-29 Thread arun
Hello, I tried to do the same thing that you have done.  Without the timestamps as rownames (xts object), there were no errors. dat1<-read.table(text=" dive_id  date    time  species    count size  1 08/06/2008 8:49:00 S.OYT    15 6  1 08/06/2008 8:49:00 S.atrovirens 1 23 ",head

Re: [R] Beta-Binomial Regression in R

2012-07-29 Thread arun
Hello, Isn't the same question posted on stackoverflow (http://stackoverflow.com/questions/11705062/generating-beta-binomial-responses-with-covariates)?  I guess you got the answer there. In general, cross-posting is not that preferred in R. A.K. - Original Message - From: Anamika C

Re: [R] Beta-Binomial Regression in R

2012-07-29 Thread David Winsemius
On Jul 28, 2012, at 3:28 PM, Anamika Chaudhuri wrote: Hi All: I am trying to generate Beta-Binomial data with regressors using R. I have used the following code to generate Beta-Binomial data. Now I want to add a covariate to the equation. I would then like to use the simulated data to