Re: [R] Appending the Column names

2012-07-29 Thread Michael Weylandt
You can't have two sets of colnames on a single data.frame What are you trying to do (big picture wise)? Michael On Jul 28, 2012, at 12:56 PM, namit saileshchowd...@gmail.com wrote: Hi Freinds, I have two data frames X,Y. I want to append both the data frames into one, along with the

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

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

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

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 that

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

[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,]

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

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] 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

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:

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

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 roma...@grinnell.edu 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

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

[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

[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

[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 system

[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

[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] 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=mlrview=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] 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:

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

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] 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). Is

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

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] 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 smartpink...@yahoo.com 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 species

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 in

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

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

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 w/ 6 levels 1,10,15,..: 1

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

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

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

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

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+switchdomains=r-project.orgsitesearch=r-project.orgbtnG=Google+Search If you know the cut-off between states no complicated modeling is necessary, using something like this

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 ny2292...@yahoo.com wrote: I am new to this field and I am currently

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(dat1) 'data.frame':6 obs.

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] 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

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 then

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 jvad...@usgs.gov wrote: Joao, Your intuition

Re: [R] extended R

2012-07-29 Thread Steve Lianoglou
Hi, On Sun, Jul 29, 2012 at 5:58 AM, hony mn.nas...@gmail.com 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:

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

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