Re: [R] R development master class: NYC, Dec 12-13

2011-11-15 Thread Stavros Macrakis
Last time, I was told that I couldn't list my R package and associated papers as a research activity with substantial impact because it was outside my official scope of work. (Even though I wrote it so I could *do* my work.) That seems wrong. My impression is that method papers were frequent

[R] [R-pkgs] mvProbit -- Multivariate Probit Models

2011-11-15 Thread Arne Henningsen
Dear R users, I am happy to announce the initial release of the mvProbit package on CRAN (version 0.1-0). This package provides tools for econometric analysis with Multivariate Probit Models. While these models can be estimated also by several other statistical software packages (e.g.

Re: [R] R development master class: NYC, Dec 12-13

2011-11-15 Thread Thomas Levine
Hmm... you know I'd love to run a study on how software and other information displays affect the speed, accuracy and reliability with which people make insights about data. Tom On Tue, Nov 15, 2011 at 4:44 PM, Stavros Macrakis macra...@alum.mit.eduwrote: Last time, I was told that I couldn't

[R] Convert back to lower triangular matrix

2011-11-15 Thread Juliet Ndukum
Given a vector; ab = seq(0.5,1, by=0.1) ab[1] 0.5 0.6 0.7 0.8 0.9 1.0 The euclidean distance between the vector elements is given by the lower triangular matrix  dd1 = dist(ab,euclidean) dd1    1   2   3   4   52 0.1                3 0.2 0.1            4 0.3 0.2 0.1        5 0.4 0.3 0.2 0.1    6

Re: [R] Convert back to lower triangular matrix

2011-11-15 Thread Jeff Newmiller
This is R-help, not the linear algebra hotline. Please stay on topic. --- Jeff NewmillerThe . . Go Live... DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live Go...

[R] Convert full matrix back to lower triangular matrix

2011-11-15 Thread Juliet Ndukum
Given a vector; ab = seq(0.5,1, by=0.1) ab[1] 0.5 0.6 0.7 0.8 0.9 1.0 The euclidean distance between the vector elements is given by the lower triangular matrix  dd1 = dist(ab,euclidean) dd1    1   2   3   4   52 0.1                3 0.2 0.1            4 0.3 0.2 0.1        5 0.4 0.3 0.2 0.1    

Re: [R] Convert back to lower triangular matrix

2011-11-15 Thread David Winsemius
On Nov 15, 2011, at 5:38 PM, Juliet Ndukum wrote: Given a vector; ab = seq(0.5,1, by=0.1) ab[1] 0.5 0.6 0.7 0.8 0.9 1.0 The euclidean distance between the vector elements is given by the lower triangular matrix dd1 = dist(ab,euclidean) dd11 2 3 4 5 2 0.1 3 0.2 0.1

Re: [R] Convert back to lower triangular matrix

2011-11-15 Thread Peter Langfelder
On Tue, Nov 15, 2011 at 2:38 PM, Juliet Ndukum jpnts...@yahoo.com wrote: Given a vector; ab = seq(0.5,1, by=0.1) ab[1] 0.5 0.6 0.7 0.8 0.9 1.0 The euclidean distance between the vector elements is given by the lower triangular matrix  dd1 = dist(ab,euclidean) dd1    1   2   3   4   52 0.1    

Re: [R] Convert back to lower triangular matrix

2011-11-15 Thread Bert Gunter
1. Please post in plain text, not HTML (as the posting guide asks!) 2. This might actually be an R question -- is ?lower.tri what you want? -- Bert On Tue, Nov 15, 2011 at 2:50 PM, Jeff Newmiller jdnew...@dcn.davis.ca.us wrote: This is R-help, not the linear algebra hotline. Please stay on

[R] if/else scope

2011-11-15 Thread Kevin Burton
What is wrong with the following? x - 1:2 if(x[1] 0) { if(x[2] 0) { print(1 2 0) } else { print(1 0) } } else {

Re: [R] if/else scope

2011-11-15 Thread Peter Langfelder
On Tue, Nov 15, 2011 at 3:46 PM, Kevin Burton rkevinbur...@charter.net wrote: What is wrong with the following? AFAIK the else needs to follow the end brace of if {} on the same line, at least at the main level. Peter x - 1:2 if(x[1] 0) {                if(x[2] 0)              

Re: [R] if/else scope

2011-11-15 Thread Ben Tupper
Hi, Check out this R-help thread from 2007. http://tolstoy.newcastle.edu.au/R/e2/help/07/06/19513.html Cheers, Ben On Nov 15, 2011, at 6:46 PM, Kevin Burton wrote: What is wrong with the following? x - 1:2 if(x[1] 0) { if(x[2] 0) {

[R] Forescasting using predict() in an object of class arimax when there is an outlier IO in the model.

2011-11-15 Thread roman_elessar1
Forescasting using predict() in an object of class arimax when there is an outlier IO in the model. Hi R users I have a problem when a use the predict() method in an object of class arimax ( These objects are the results of the implementation of the function arimax() from the TSA library) . The

Re: [R] grid.arrange, grid.layout - legend, global y axis title

2011-11-15 Thread baptiste auguie
Hi, Please don't cross post. It seems that ggplotGrob has been replaced by new functions. You can define it as ggplotGrob - function(x) ggplot2:::gtable_gTree(ggplot2:::ggplot_gtable(x)) and it seems to work as before with grid.arrange(). HTH, baptiste On Wed, Nov 16, 2011 at 3:26 AM,

[R] Conversion of symmetry matrix into a vector

2011-11-15 Thread ram basnet
Dear R users, I am not good in R-language programming. So, i need your help. I want to convert my lower-triangle value of symmetry matrix into a vector with their row and column name. I found a function called sm2vec in corpcor package but it give only a vector of values but not row and column

[R] R Program

2011-11-15 Thread PRI
Can someone help me out with this problem? Write an R program to draw 10,000 random samples using the Polya Urn Scheme where the initial contents of the urn is one white ball and one black ball. Accumulate the proportion of white balls after each draw into one vector. Upload your R script file.

Re: [R] Reading data/variables

2011-11-15 Thread Steven Yen
A follow-up on the data/variable issue I posted earlier: Here was what I did, which might was obviously causing the problem: I inserted the following line in my file Rprofile.site: setwd(z:/R) Then, as soon as I run R (before I read any data) I issue summary(mydata) I get summary

[R] Create a function with multiple object as an output

2011-11-15 Thread anaraster
I've seen some questions regarding the output of multiple objects from a function, however the suggestions all end up suggesting the use of return(list(result1=result1, result2=result2 , result3=result3)). How can I return multiple objects that are 2 big to be added to a list? -- View this

Re: [R] Create a function with multiple object as an output

2011-11-15 Thread PRI
Write an R program to draw 10,000 random samples using the Polya Urn Scheme where the initial contents of the urn is one white ball and one black ball. Accumulate the proportion of white balls after each draw into one vector. Upload your R script file. HINT: You can simulate sampling one ball from

Re: [R] R Program

2011-11-15 Thread David Winsemius
On Nov 15, 2011, at 3:57 PM, PRI wrote: Can someone help me out with this problem? Write an R program to draw 10,000 random samples using the Polya Urn Scheme where the initial contents of the urn is one white ball and one black ball. Accumulate the proportion of white balls after each

Re: [R] Difference between two time series

2011-11-15 Thread Sarwarul Chy
Hello Michael, Thanks for your reply. What I want to do is something like this? For example, I have a continuous time series y=x(t), and another discrete time series z=w(t). Xdiff(i)=Max. difference between x(t) and w(t) in interval i Ndiff(i)=Min. difference between x(t) and w(t) in interval

Re: [R] if/else scope

2011-11-15 Thread P.B. Lecavalier
On 11/15/2011 06:46 PM, Kevin Burton wrote: What is wrong with the following? x- 1:2 if(x[1] 0) { if(x[2] 0) { print(1 2 0) } else {

[R] Problem with substr

2011-11-15 Thread sophy1987
Hi, everyone When I ran this cript, There is Error in substring(tmp.subject, tmp.end[ex] + 1, tmp.start[ex + 1] - 1) : invalid substring argument(s) Could someone figure out what the problem is? for(i in 1:length(genebody[,1])){ tmp.id-as.vector(genebody[i,1])

[R] Reading data/variables

2011-11-15 Thread Steven Yen
Can someone help me with this variable/data reading issue? I read a csv file and transform/create an additional variable (called y). The first set of commands below produced different sample statistics for hw11$y and y In the second set of command I renameuse the variable name yy, and sample

Re: [R] gsDesign

2011-11-15 Thread Dongli Zhou
Thank you so much for the help! It is really appreciated! Dongli On Nov 15, 2011, at 9:19 AM, Marc Schwartz marc_schwa...@me.com wrote: Hi Dongli, Sorry for the delay in following up. You might want to read the dsDesignManual.pdf document, which is available in the 'inst/doc' folder in

Re: [R] Create a function with multiple object as an output

2011-11-15 Thread Sarah Goslee
Hi, On Tue, Nov 15, 2011 at 4:10 PM, anaraster rrast...@gmail.com wrote: I've seen some questions regarding the output of multiple objects from a function, however the suggestions  all end up suggesting the use of return(list(result1=result1, result2=result2 , result3=result3)). How can I

Re: [R] Reading data/variables

2011-11-15 Thread Sarah Goslee
Hi, The obvious answer is don't use attach() and you'll never have that problem. And see further comments inline. On Tue, Nov 15, 2011 at 6:05 PM, Steven Yen s...@utk.edu wrote: Can someone help me with this variable/data reading issue? I read a csv file and transform/create an additional

Re: [R] Problem with substr

2011-11-15 Thread Sarah Goslee
Hi, On Tue, Nov 15, 2011 at 4:31 PM, sophy1987 yatingcheng.mar...@gmail.com wrote: Hi, everyone When I ran this cript,  There is Error in substring(tmp.subject, tmp.end[ex] + 1, tmp.start[ex + 1] - 1) :  invalid substring argument(s) Could someone figure out what the problem is? Not

Re: [R] Create a function with multiple object as an output

2011-11-15 Thread Rolf Turner
On 16/11/11 10:10, anaraster wrote: I've seen some questions regarding the output of multiple objects from a function, however the suggestions all end up suggesting the use of return(list(result1=result1, result2=result2 , result3=result3)). How can I return multiple objects that are 2 big to

Re: [R] pairs(), no axis labels/values for upper panel?

2011-11-15 Thread B77S
Steffen, Did you ever have luck getting rid of the tick marks?... I like your idea and have modified it, but yes, the tick marks need to go. Steffen Fleischer wrote: Dear all, I want to draw a graph that contains the scatterplot matrix in the lower panel and coefficients in the upper

Re: [R] Question on exercise

2011-11-15 Thread Ben Bolker
Jan Kraner jan.kraner at gmx.ch writes: Dear Sir/Madam: Could you please help me with this exercise and how I could solve it in R?: Attachment (PastedGraphic-1.pdf): application/pdf, 84 KiB I especially have problems with the selection of the 40 datasets, once taken the 110 (because

Re: [R] Question on exercise

2011-11-15 Thread R. Michael Weylandt michael.weyla...@gmail.com
Copying David W from another thread: We don't do homework here. In fact, your instructor is probably watching. That said, try googling or, even better, using the rseek search. (Or even better-er RSiteSearch() in your terminal) Michael On Nov 15, 2011, at 4:16 PM, Jan Kraner jan.kra...@gmx.ch

Re: [R] Create a function with multiple object as an output

2011-11-15 Thread Ben Bolker
PRI priyanka.246 at gmail.com writes: Write an R program to draw 10,000 random samples using the Polya Urn Scheme where the initial contents of the urn is one white ball and one black ball. Accumulate the proportion of white balls after each draw into one vector. Upload your R script file.

[R] knn density estimation in R

2011-11-15 Thread lyh123985458
How to do the nearest neighbor density estimation in R? Is there any function/package to perform k-Nearest Neighbor based density estimation in R Thank you very much! Thanks to a dozen of volunteer moderators! [[alternative HTML version deleted]]

Re: [R] Conversion of symmetry matrix into a vector

2011-11-15 Thread R. Michael Weylandt
Try something like this: X = matrix(rnorm(25), 5); colnames(X) - rownames(X) - letters[1:5] Y - dist(X) library(reshape) melt(as.matrix(Y)) Michael On Tue, Nov 15, 2011 at 7:03 PM, ram basnet basnet...@yahoo.com wrote: Dear R users, I am not good in R-language programming. So, i need your

Re: [R] Difference between two time series

2011-11-15 Thread R. Michael Weylandt
Can you post working examples of your data using the dput() function? There are so many types of time series in R and so many different things you could mean that it's just easier to work with real data. Michael On Tue, Nov 15, 2011 at 4:28 PM, Sarwarul Chy sarwar.sha...@gmail.com wrote: Hello

Re: [R] Help with error: no acceptable C compiler found in $PATH

2011-11-15 Thread R. Michael Weylandt
Yes, you probably need some sort of C compiler, but why can't you just download the appropriate binary directly? I just did on OS X 10.5.8 (admittedly for R 2.13.2, not 2.14) with no problems. The output of sessionInfo() install.packages(RCurl) if you don't mind please. Thanks, Michael On

[R] plotting a double y axis when x and y lengths differ

2011-11-15 Thread Vinny Moriarty
Hello All, Many thanks to the help I have received so far. Here is an example data set I hope to plot Data1 Year Data SE 1 20052 0.01 2 20064 0.01 3 20075 0.01 4 20082 0.01 5 20093 0.01 6 20106 0.01 Data2 Year Data SE 1 2006 32 1 2 2007 100 2 3 2008 60 4

Re: [R] equal spacing of the polygons in levelplot key (lattice)

2011-11-15 Thread Dennis Murphy
Hi: Does this work? # library('lattice') levs - as.vector(quantile(volcano,c(0,0.1,0.5,0.9,0.99,1))) levelplot(volcano, at = levs, colorkey = list(labels = list(at = levs, labels = levs) )) HTH, Dennis On Tue, Nov 15, 2011 at 1:12

Re: [R] plotting a double y axis when x and y lengths differ

2011-11-15 Thread R. Michael Weylandt michael.weyla...@gmail.com
Does plotCI take an xlim argument? Something like plotCI(..., xlim=range(Data1$Year, Data2$Year)) Michael On Nov 15, 2011, at 11:51 PM, Vinny Moriarty vwmoria...@gmail.com wrote: Hello All, Many thanks to the help I have received so far. Here is an example data set I hope to plot

Re: [R] Convert full matrix back to lower triangular matrix

2011-11-15 Thread Dennis Murphy
Hi: (1) Here is why your e-mails look mangled on this list: [[alternative HTML version deleted]] R-help is a text-only list, so please change your mailer's settings to send ASCII text rather than HTML. (2) The print method you see displayed in dd1 is equivalent to the following: ddm *

[R] Theil decomposition

2011-11-15 Thread Kitty Lee
I came across the package 'ineq' that computes a variety of inequality measures (e.g. gini, theil etc). I want to compute the Theil index (racial segregation) and decompose the total into sub-components (by geog levels). I think the package doesn't report the decomposition (correct me if I'm

[R] Determining second coordinate from distance

2011-11-15 Thread Hadassa Brunschwig
Hi, I am looking for a function which takes as input the spatial coordinate in longitute/latitude and the distance in metric and returns the second coordinate (of course depending on an angle as well). So far I have only found functions which calculate the distance between two coordinates.

[R] fitted.values in lm function

2011-11-15 Thread arunkumar1111
Hi i have few clarification regarding the output from the fitted(object). I'm creating a object using the following formula dataset - read.csv(~/data.csv) obj = lm (formula=log(y+1)~log(x1+1)+(log(x2+1)+log(x3+1), data=dataset) fitted(obj) here the output of the fitted(obj) considers the

Re: [R] R forum for only Statistics

2011-11-15 Thread geeknick
Hi there I need some [URL=http://calculate-conditional-and-distributio.blogspot.com/2011/11/blog-post.html]Statistics Help[/URL]!! What statistical tests should i use/consider (pref in SPSS) in regards to an invertebrate survey i have carried out comparing urban street trees to urban park

Re: [R] fitted.values in lm function

2011-11-15 Thread Joshua Wiley
log(x1+1) just like your formula best, Josh On Tue, Nov 15, 2011 at 9:54 PM, arunkumar akpbond...@gmail.com wrote: Hi  i have few clarification regarding the output from the fitted(object). I'm creating a object using the following formula dataset  - read.csv(~/data.csv) obj = lm

Re: [R] global optimisation with inequality constraints

2011-11-15 Thread geeknick
Please advice on the package I should use to run a linear regression model (weighted least squared) with linear equality constraint. I initially tried constrOptim but it turned out that it only supported http://solve-graph-linear-inequalities.blogspot.com/2011/11/blog-post.html linear inequality

Re: [R] plotting a double y axis when x and y lengths differ

2011-11-15 Thread Rolf Turner
Any time you use ``par(new=TRUE)'', you're probably doing something wrong. Avoid it like the plague. Following up on Michael Weylandt's post: plotCI() indeed has an xlim argument (i.e. the ... arguments which get passed to plot.default() can include xlim). Note that plotCI() also has an

Re: [R] global optimisation with inequality constraints

2011-11-15 Thread Dennis Murphy
(1) Please do not hijack another thread to ask an off-topic question - start a new one instead. (2) Your question refers to a topic called 'restricted least squares'. Search the R-help archives on that subject and you should get a number of answers. Someone answered a similar question here within

<    1   2