Re: [R] Discrete choice model maximum likelihood estimation

2012-05-14 Thread Rui Barradas
Ok, I forgot to say that 't' is also an R function, the matrix transpose. Sorry, but after 'par' I thought (in my mind) I had said it when in fact I even talked about 't'! Use 'tt'. If 'tt' is a vector you must first define it, in your code it doesn't exist. That's why R searches for and finds an

Re: [R] Discrete choice model maximum likelihood estimation

2012-05-14 Thread infinitehorizon
Hello Rui, First of all, thanks a lot! 1. I changed par to param, 2. t is a variable too, a binary one, b is the parameter associated to it, 4. Yes, this is where I am stuck actually. I fixed the code for likelihood functions as follows, but still getting the same error: L3 -

[R] To summary data

2012-05-14 Thread mrzung
hi all, my data is here, -- data clss number 1A 1 2B 2 3C 3 4A 4 5B 5 6C 6 7A 7 8B 8 9C 9

[R] R i386 2.15.0 'gogarch' package issue

2012-05-14 Thread Cren
Hi all, I've just downloaded and installed the latest R 32-bit version plus RExcel and R Commander. I'm having several problems in loading gogarch package: The command *library(gogarch, pos=4)* returns *ERROR: package/namespace load failed for 'gogarch'* The command *require(gogarch)* returns

Re: [R] R i386 2.15.0 'gogarch' package issue

2012-05-14 Thread Cren
If I use gogarch_0.7-1 the command *require(gogarch)* returns the following error: *Error in get(.packageName, where) : cannot allocate memoby block of size 3.2 Gb* -- View this message in context: http://r.789695.n4.nabble.com/R-i386-2-15-0-gogarch-package-issue-tp4629888p4629889.html Sent from

Re: [R] Why can we combine design matrix and data-frame in R?

2012-05-14 Thread peter dalgaard
On May 14, 2012, at 02:24 , Luna wrote: Thanks! Do you think if the correctness of the such results could be generalized to other future cases? If correctly generalized, yes (Apologies for being slightly facetious; the point is that the properties you build on are part of the

Re: [R] R i386 2.15.0 'gogarch' package issue

2012-05-14 Thread Cren
Ok, solved. If anyone had the same problem, just install the last gogarch pacakge (vers. 0.7-1) and restart R + R Commander after the package installation. When you've restarted it, the command *require(gogarch)* should load fastICA package in the end and it will work ;) -- View this message in

Re: [R] To summary data

2012-05-14 Thread Berend Hasselman
On 14-05-2012, at 07:14, mrzung wrote: hi all, my data is here, -- data clss number 1A 1 2B 2 3C 3 4A 4 5B 5 6C 6 7A 7

[R] labels in hist and ranges of an axis

2012-05-14 Thread carol white
Hi, How is it possible to change the character size of labels in a histogram (labels parameter in hist)? How is it possible to display the ranges of an axis in 3 significant digits? Thanks  Carol [[alternative HTML version deleted]] __

[R] Aligning time series

2012-05-14 Thread Matthew Johnson
Sir, I have large data sets of economic indicators and would like to align them to a reference series - say the unemployment rate or industrial production. Is there a canned routine that returns the optimal lead / lag according to some (or a variety) of algorithims? not all series will be of the

Re: [R] file path

2012-05-14 Thread Wincent
Thanks for the suggestion. The file name in my case is Chinese, which makes the regular expression less useful. Anyway, I would like to pose a followup question. I have a character string of ABC\D, and want to strip away the \ and want a returned character of ABCD. How can I do it with gsub() ?

Re: [R] labels in hist and ranges of an axis

2012-05-14 Thread Gerrit Eichner
Hi, Carol! How is it possible to change the character size of labels in a histogram (labels parameter in hist)? See ?par - cex.axis How is it possible to display the ranges of an axis in 3 significant digits? See ?axis and its argument at together with ?signif Hth -- Gerrit

Re: [R] labels in hist and ranges of an axis

2012-05-14 Thread carol white
Thanks Gerrit. cex.axis changes the size of the axis label. I want to change the size of data labels that are displayed at the top of the bars in the histogram. What we obtain by histogram(...,labels=T) Any clue? Cheers, From: Gerrit Eichner

Re: [R] file path

2012-05-14 Thread Wincent
No, but make.names only makes syntactically valid names of R object. I think I have mentioned what I did in the first email, although no concrete example was provided. Let me explain a bit. I have a text file like this, Title:AA? Content: 1xx1 Title:BA Content: 2xx2 I want

Re: [R] file path

2012-05-14 Thread Baoqiang
This works on Mac: str - abc/d gsub(/, , str) Return: abcd Sent from my iPhone On May 14, 2012, at 4:28 AM, Wincent ronggui.hu...@gmail.com wrote: Thanks for the suggestion. The file name in my case is Chinese, which makes the regular expression less useful. Anyway, I would like to pose

Re: [R] file path

2012-05-14 Thread Wincent
Emm, my bad. I meant str - abc\d. Any ideas? On 14 May 2012 18:02, Baoqiang bqcaom...@gmail.com wrote: This works on Mac: str - abc/d gsub(/, , str) Return: abcd Sent from my iPhone On May 14, 2012, at 4:28 AM, Wincent ronggui.hu...@gmail.com wrote: Thanks for the suggestion. The

[R] How to interpret an ANOVA result?

2012-05-14 Thread Robert Latest
Hello all, here's a real-world example: I'm measuring a quantity (d) at five sites (site1 thru site5) on a silicon wafer. There is a clear site-dependence of the measured value. To find out if this is a measurement artifact I measured the wafer four times: twice in the normal position (posN), and

Re: [R] file path

2012-05-14 Thread Berend Hasselman
On 14-05-2012, at 12:07, Wincent wrote: Emm, my bad. I meant str - abc\d. Any ideas? gsub(, , str) Berend __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] update.formula simplify = FALSE

2012-05-14 Thread SUBIRANA CACHINERO, ISAAC
Dear R users, I am building a packaged where I want to use the update.formula function, but modifying the option simplify = FALSE to TRUE. Since this is not an argument, I modified the code and created a new function as: update.formula2 - function (old, new, ...) { tmp -

[R] Error using rect.hclust

2012-05-14 Thread Suhaila Haji Mohd Hussin
Hello. I'm trying to see clusters into a certain number of classes as seeing they are too many of them see entirely on dendrogram. So I did the following: d - dist(as.matrix(medData)) hc - hclust(d, method=average) plot(hc) plot(hc, hang = -1) But when it came to this: rect.hclust(hc, 3)Error

Re: [R] labels in hist and ranges of an axis

2012-05-14 Thread Gerrit Eichner
Carol, it is not clear to me which function histogram() you use (package lattice, package histogram ...?), but -- if it is not a lattice function -- a quick hack _might_ be to use par( cex = 2) or whatever magnification you wish followed by a reduction of cex.axis, cex.labels etc.

Re: [R] Aligning time series

2012-05-14 Thread Gabor Grothendieck
On Mon, May 14, 2012 at 5:09 AM, Matthew Johnson mcoog...@gmail.com wrote: Sir, I have large data sets of economic indicators and would like to align them to a reference series - say the unemployment rate or industrial production. Is there a canned routine that returns the optimal lead /

Re: [R] Plotmath bug or my misunderstanding?

2012-05-14 Thread Uwe Ligges
On 13.05.2012 14:37, Rui Barradas wrote: Hello, It seems like a feature. When trying the example below with more atops, the two bottom most lines, and only those two lines, feature character expansion relative to the default size and relative to the line before last. plot(1, type=n,

Re: [R] Overlapping area script

2012-05-14 Thread Uwe Ligges
On 11.05.2012 11:06, achughes wrote: Dear All I would really appreciate some help with a script which a colleague wrote for me (attached), but I am having problems running (and have not been able to contact my colleague). I think you have to contact your colleague or you have to provide

[R] Re : R i386 2.15.0 'gogarch' package issue

2012-05-14 Thread Pascal Oettli
Hello, Probably you should try: update.packages(checkBuilt=TRUE) install.packages('gogarch', dependencies=TRUE) Best Regards, Pascal - Mail original - De : Cren oscar.soppe...@bancaakros.it À : r-help@r-project.org Cc : Envoyé le : Lundi 14 mai 2012 16h06 Objet : [R] R i386 2.15.0

[R] Extract Variance Components

2012-05-14 Thread robgriffin247
Hi, I'm still having problems putting the variance components of my model in to a data frame, it is a continuation of this discussion, http://r.789695.n4.nabble.com/ANOVA-problem-td4609062.html, but now focussed on the problem of extracting variance components. I have got my mixed effects model

[R] Data read as labels

2012-05-14 Thread barb
Hey guys, i have a strange problem reading a .csv file. Seems not to be covered by the usual read.csv techniques. The relevant data i want to use, seems to be saved as the label of the data point. Therefore i can not really use it spec-EU2001

[R] rollmax.zoo : column names NULL

2012-05-14 Thread Giles
Hi I am comparing the output of rollmax in two versions of R. In the current version, the column names are 'lost' ie NULL in the output; in the earlier version they were retained. Function rollapply has the same behaviour as before. Is this change in rollmax behaviour deliberate? Thanks

[R] New to R

2012-05-14 Thread Ronald McDowell
I am new to R and starting to explore its functionality. I wondered if anyone could advise whether R supports non-linear canonical correlation and/or the specification of models using alternating least squares? Thanks Ron [[alternative HTML version deleted]]

Re: [R] Discrete choice model maximum likelihood estimation

2012-05-14 Thread infinitehorizon
Hello again, I changed the name to tt. and for a and tt actually I was getting them from data, I didn't put them here in the question. Now I restructured my code and below I copy the full code, I tried many things but still getting the same error, I don't understand where is the mistake. I also

Re: [R] Overlapping area script

2012-05-14 Thread achughes
Dear Uwe I did previously attach the script, but I include it again here in addition to current and future projections (which come from different directories).Help would very much be appreciated! Alice Date: Mon, 14 May 2012 04:52:04 -0700 From: ml-node+s789695n4629915...@n4.nabble.com To:

[R] cross validation in glmnet

2012-05-14 Thread yan
I am using cv.glmnet from glmnet package for logistic regression. my dataset is very imbalanced, 5% sample from one group, the rest from the other. I'm wondering when doing cv.glmnet for choosing lambda, is every fold having the same ratio for two groups(every fold has 5% sample from one group,

[R] Add column from other columns data.

2012-05-14 Thread Yellow
Hi everyone, I am having some problems with making a new colomn wit data in it. I have this one column named: Fulfilled Fulfilled 1 1 0 1 1 1 1 0 0 1 And now I would like to add another colum to my .csv file (Finished) In this Finished column I would like to have Yes or No. Where in colomn

Re: [R] Overlapping area script

2012-05-14 Thread Uwe Ligges
That script is simple broken, the object bat.asc is used but not defined anywhere. Since the code is also not really optimal, you may want to ask an expert who knows how to write such script and how to test if their code work in the end. Best, Uwe Ligges On 14.05.2012 14:12, achughes

[R] Confidence interval for Whittle method

2012-05-14 Thread Barun Saha
Hello, How could we get the confidence interval when using the whittleFit method from fArma package? -- Thanks, Barun Saha __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Re : R i386 2.15.0 'gogarch' package issue

2012-05-14 Thread Cren
Pascal Oettli-2 wrote Hello, Probably you should try: update.packages(checkBuilt=TRUE) install.packages('gogarch', dependencies=TRUE) Best Regards, Pascal Dear Pascal Oettli-2, thank you for your suggestment; I was not aware of that command and it will be very useful when new

Re: [R] Calculating all possible ratios

2012-05-14 Thread genome1976
Hi Rui, Thanks once again for all the help. I need to ask for one more help from you. I have two matrices, with probesets as rows and samples as columns. The samples in the two matrices are matched (from the same animal but two different tissues). I want to create a correlation matrix of

Re: [R] rollmax.zoo : column names NULL

2012-05-14 Thread Achim Zeileis
Giles, thanks for the bug report: I am comparing the output of rollmax in two versions of R. In the current version, the column names are 'lost' ie NULL in the output; in the earlier version they were retained. Yes, this was an error. I just fixed it in the devel version on R-Forge.

Re: [R] Add column from other columns data.

2012-05-14 Thread Sarah Goslee
Assuming you actually have a data frame or matrix, and not a csv file, ifelse() is the general solution to your problem. Sarah On May 14, 2012, at 8:45 AM, Yellow s1010...@student.hsleiden.nl wrote: Hi everyone, I am having some problems with making a new colomn wit data in it. I have

[R] Rjava on Ubuntu quantal

2012-05-14 Thread Hasan Diwan
I just upgraded to Ubuntu Quantal from Precise and RJava stopped working, log follows: 0}% /usr/bin/find $HOME/workspace/FinanceOCR/visualizations/ -name '*R' -print | /usr/bin/xargs -n 1 -i% /usr/bin/Rscript % $1 [~] Loading required package: RJDBC Loading required package: methods Loading

Re: [R] Why can we combine design matrix and data-frame in R?

2012-05-14 Thread Michael
Oh, so we can always combine model matrices and formulas in regression in R? Thanks! On Mon, May 14, 2012 at 2:41 AM, peter dalgaard pda...@gmail.com wrote: On May 14, 2012, at 02:24 , Luna wrote: Thanks! Do you think if the correctness of the such results could be generalized to

[R] How do I do group wise clustering in R?

2012-05-14 Thread Luna
How do I do group wise clustering in R? Hi all, I have N x K data matrix, where N is the number of observations, K is the number of variables. The N observations fall into M categories or groups. Now I want to cluster the groups, instead of the observations, how do I do that? i.e. the

Re: [R] Add column from other columns data.

2012-05-14 Thread John Kane
Something along the lines of dat2 - ifelse( dat1==1 , yes, no) should do it. John Kane Kingston ON Canada -Original Message- From: s1010...@student.hsleiden.nl Sent: Mon, 14 May 2012 05:45:38 -0700 (PDT) To: r-help@r-project.org Subject: [R] Add column from other columns data.

Re: [R] Random forests prediction

2012-05-14 Thread Liaw, Andy
I don't think this is so hard to explain. If you evaluate AUC using either OOB prediction or on a test set (or something like CV or bootstrap), that would be what I expect for most data. When you add more variables (that are, say, less informative) to a model, the model has to look harder to

Re: [R] Add column from other columns data.

2012-05-14 Thread Yellow
That worked. Thanks. :) -- View this message in context: http://r.789695.n4.nabble.com/Add-column-from-other-columns-data-tp4629921p4629937.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

Re: [R] To summary data

2012-05-14 Thread John Kane
One way is to use the reshape2 package === library(reshape2) dcast(xx, A ~ . , sum) John Kane Kingston ON Canada -Original Message- From: mrzun...@gmail.com Sent: Sun, 13 May 2012 22:14:02 -0700 (PDT) To:

[R] ggplot2: Dendrogram text position

2012-05-14 Thread Brian Smith
Hi, I was trying to create a dendrogram using ggplot2. Everything seems to be looking ok except that the text labels are too close to the dendrogram (in the example below, 'a','b', ..). Is there a way that I can put a little gap between where the dendrogram ends and the label begins? thanks!!

Re: [R] file path

2012-05-14 Thread David Winsemius
On May 14, 2012, at 6:35 AM, Berend Hasselman wrote: On 14-05-2012, at 12:07, Wincent wrote: Emm, my bad. I meant str - abc\d. Any ideas? gsub(, , str) #1: One cannot execute: str - abc\d , at least on my machine, since that throws an error because \d is an unrecognized escape.

[R] phyloclim could not be installed in linux - problems on tkrplot dependence

2012-05-14 Thread Mao Jianfeng
Dear R-helpers, Christoph (author of phyloclim) and Luke (author of tkrplot), I would like to get your helps on installing of phyloclim in Ubuntu linux. It seems a package named 'tkrplot' could not be installed at firstly, then packages depends on it could not be installed latter. As I have

Re: [R] range segment exclusion using range endpoints

2012-05-14 Thread Ben quant
Wow, I'll have to study this one for a bit. Thanks! Ben On Sat, May 12, 2012 at 3:09 PM, William Dunlap wdun...@tibco.com wrote: Here is some code that I've been fiddling with for years (since I wanted to provide evidence that our main office needed more modems and wanted to show how often

Re: [R] Data read as labels

2012-05-14 Thread David Winsemius
On May 14, 2012, at 5:33 AM, barb wrote: Hey guys, i have a strange problem reading a .csv file. Seems not to be covered by the usual read.csv techniques. The relevant data i want to use, seems to be saved as the label of the data point. Therefore i can not really use it spec-EU2001

Re: [R] range segment exclusion using range endpoints

2012-05-14 Thread Ben quant
Great solution! Thanks! Ben On Sat, May 12, 2012 at 12:50 PM, jim holtman jholt...@gmail.com wrote: Here is an example of how you might do it. It uses a technique of counting how many items are in a queue based on their arrival times; it can be used to also find areas of overlap. Note

Re: [R] No Data in randomForest predict

2012-05-14 Thread Liaw, Andy
It doesn't: You just get an error if there are NAs in the data; e.g., R rf1 = randomForest(iris[1:4], iris[[5]]) R predict(rf1, newdata=data.frame(Sepal.Length=1, Sepal.Width=2, Petal.Length=3, Petal.Width=NA)) Error in predict.randomForest(rf1, newdata = data.frame(Sepal.Length = 1, :

[R] glmnet speed

2012-05-14 Thread yan
I'm using glmnet for logistic regression, I got a fairly sparse dataset, 20,000 samples(very imbalanced too, 5% from one group), 1500 variables,. the code have beed running for 2 hours, still waiting for result, I am doing lasso here(alpha=1), my computer is core 2 due CPU @3Ghz, 4GB ram, why

[R] New version of cloudRmpi

2012-05-14 Thread Barnet Wagman
cloudRmpi v 1.2 is now available on CRAN. cloudRmpi is a means for doing parallel processing in R, using MPI on a cloud-based network. It currently supports the use of Amazon's EC2 cloud computer service. Changes in v 1.2: Support for RStudio. RStudio Server is available on new AMIs.

Re: [R] Help with writing data to csv

2012-05-14 Thread David L Carlson
If you have a data.frame and you want a table in Microsoft Word, the quickest path (without additional packages and assuming you are using Windows) is write.table(DataFrameName, file=clipboard, sep=\t, row.names=FALSE) # If the file is large, you may need clipboard-128 instead of clipboard. Now

Re: [R] how to write data using xlsReadWrite

2012-05-14 Thread David L Carlson
You didn't tell us what your problem is, but it probably relates to the fact that mydata is never defined. -- David L Carlson Associate Professor of Anthropology Texas AM University College Station, TX 77843-4352 -Original Message- From:

Re: [R] How do I do group wise clustering in R?

2012-05-14 Thread David L Carlson
Look at the aggregate function to create a new data.frame in which you have M rows that have the means of the K variables for each group. Then use cluster analysis to cluster the M groups. -- David L Carlson Associate Professor of Anthropology Texas AM

[R] select data

2012-05-14 Thread Andrea Sica
Dear all, I am sure it won't be difficult for you!! I need to calculate the average among variables for the single units of my dataset. But, while doing it, I need to do not consider some values. To better explain, think like there are two units and three variables: V1V2 V3 [1] 3

Re: [R] New to R

2012-05-14 Thread David L Carlson
You will find functions such as these in the thousands of packages that are available once you have installed R. You can use rseek.org to search for specific topics. Good overviews are found in the CRAN Task Views (from the main R webpage, click on CRAN, select a mirror host, and then select Task

Re: [R] Post stratification weights in survey package in R

2012-05-14 Thread Thomas Lumley
On Sun, May 13, 2012 at 7:10 PM, Ruijie breakaw...@gmail.com wrote: Hi all, I have data collected from a survey administered on a subset of the population. I also have the population proportions of variables such as gender, race and housing type. I would like to combine the weights from each

Re: [R] Discrete choice model maximum likelihood estimation

2012-05-14 Thread infinitehorizon
By the way, in my last post I forgot to return negative of llfn, hence the llfn will be as follows: llfn - function(param) { bx - param[1] b1 - param[2] b2 - param[3] b3 - param[4] b - param[5] lL1 - log(L1(bx,b1,b2,b,tt)) lL2 - log(L2(bx,b1,b2,b3,b,tt)) lL3 - log(L3(bx,b1,b2,b3,b,tt))

Re: [R] Random forests prediction

2012-05-14 Thread matt
But shouldn't it be resolved when I set mtry to the maximum number of variables? Then the model explores all the variables for the next step, so it will still be able to find the better ones? And then in the later steps it could use the (less important) variables. Matthijs -- View this message

Re: [R] Data read as labels

2012-05-14 Thread barb
Hey David, thanks for your fast reply, i really appreciate that you answer so many posts. Unfortunately it´s not that easy. Try to operate with the output: e.g file-read.csv2(tmp,sep=;,skip=5) a-(relevant-file[,2][1]) a*5 # or as.numeric(relevant-file[,2][1]) a is saved in the workspace as

Re: [R] Extract Variance Components

2012-05-14 Thread robgriffin247
Well, I'm going to reply to my own thread with a solution here, turns out one attempt we made last week nearly had it, slight adjustment made it work. For anyone that is interested / in the future wants to achieve the same thing *varcomp - matrix(nrow=0, ncol=3) for (i in

Re: [R] Discrete choice model maximum likelihood estimation

2012-05-14 Thread infinitehorizon
Hello again, You are absolutely right about probabilities.. Thanks for reminding me about that. I did exactly how you said but in the end I receive the error : objective function in optim evaluates to length 12 not 1. I checked how llfn give a vector instead of scalar, but couldn't figure it

Re: [R] Discrete choice model maximum likelihood estimation

2012-05-14 Thread Rui Barradas
Hello, again. Bug report: 1. Your densities can return negative values, 1 - exp(...) 0. Shouldn't those be 1 PLUS exp()? P3 - function(bx,b3,b,tt) { P - exp(bx*x+b3+b*(tt == 1))/(1+exp(bx*x+b3+b*(tt == 1))) return(P) } And the same for P2 and P1? 2. Include 'a' and 'tt' as

Re: [R] select data

2012-05-14 Thread R. Michael Weylandt
This was actually discussed about a week and a half ago with many good solutions offered, but I think the most idiomatic would be something like this: apply(dataset, 1, function(x) mean(x[x0])) The reasons I like it: i) It uses the apply function to do the same operation row-wise (that's what

[R] How to apply a function to a multidimensional array based on its indices

2012-05-14 Thread math_daddy
Hello. I have a 4 dimensional array and I want to fill in the slots with values which are a function of the inputs. Through searching the forums here I found that the function outer is helpful for 2x2 matrices but cannot be applied to general multidimensional arrays. Is there anything which can

Re: [R] Data read as labels

2012-05-14 Thread Krzysztof Mitko
On Mon, May 14, 2012, at 02:33, barb wrote: Hey guys, i have a strange problem reading a .csv file. Seems not to be covered by the usual read.csv techniques. The relevant data i want to use, seems to be saved as the label of the data point. Therefore i can not really use it

Re: [R] select data

2012-05-14 Thread David Winsemius
On May 14, 2012, at 12:32 PM, Andrea Sica wrote: Dear all, I am sure it won't be difficult for you!! I need to calculate the average among variables for the single units of my dataset. But, while doing it, I need to do not consider some values. To better explain, think like there are two

Re: [R] Discrete choice model maximum likelihood estimation

2012-05-14 Thread Rui Barradas
Once again, sorry. I had a different llfn in my R session and it messed with yours. llfn - function(param, a, tt) { llfn - sum((a==1)*lL1+(a==2)*lL2+(a==3)*lL3) # sum of logs, it's a log-likelihood. return(-llfn) } Rui Barradas infinitehorizon wrote Hello again, You are absolutely

Re: [R] select data

2012-05-14 Thread Andrea Sica
Thank you all. Really! I have used the following function: apply(dfrm, 1, function(x) mean(x[x=0]) ) Someone of you even gave me a few interesting explanations about why to use it. Still thank you all. Andrea On Mon, May 14, 2012 at 6:52 PM, David Winsemius dwinsem...@comcast.netwrote: On

Re: [R] How to apply a function to a multidimensional array based on its indices

2012-05-14 Thread Rui Barradas
Hello, Try K - array(0,dim=c(2,2,2,2)) #dimensions will be much larger for(x1 in 1:2){ for(y1 in 1:2) { for(x2 in 1:2){ for(y2 in 1:2) { K[x1,y1,x2,y2] - x1*y2 - sin(x2*y1) #this is just a dummy function. } } } } fun - function(x){ x1 - x[1]

Re: [R] Data read as labels

2012-05-14 Thread David Winsemius
On May 14, 2012, at 11:23 AM, barb wrote: Hey David, thanks for your fast reply, i really appreciate that you answer so many posts. Unfortunately it´s not that easy. Try to operate with the output: e.g file-read.csv2(tmp,sep=;,skip=5) a-(relevant-file[,2][1]) a*5 # or

Re: [R] Discrete choice model maximum likelihood estimation

2012-05-14 Thread infinitehorizon
Of course, that was the trick! It works now. Thank you very much Rui, I am very grateful. I hope this thread will help others as well. Best, Rui Barradas wrote Once again, sorry. I had a different llfn in my R session and it messed with yours. llfn - function(param, a, tt) { llfn

[R] using rcom to control Power Point - problem to set a property

2012-05-14 Thread Mark Heckmann
Hi all, I try to convert the VBA code below to run it from R using rcom and Power Point. The VBA code creates a shape and will move the handle of the the shape to another position. This fails using rcom and I do not understand what I am doing wrong... ### VBA ###

Re: [R] New to R

2012-05-14 Thread Bert Gunter
I usually try google searches first. While not always successful, I am frequently surprised by how well it does. I strongly second the use of CRAN task views. IMHO, some fine folks have volunteered their time and efforts to produce this very well written series of guides to what's in R. It and

Re: [R] Help with V function in igraph

2012-05-14 Thread Gábor Csárdi
Something weird must be going on in your s641_social object. Can you just simply check that the vertex names look OK with 'V(s641_social)$name'? If they look good, then can you send me the s641_social object in private? (Or part of it, assuming a part is enough to reproduce the problem.) Best,

Re: [R] select data

2012-05-14 Thread David L Carlson
This overwrites the data so you might want to create a copy first. example - data.frame(V1=c(3, -1), V2=c(-2, 4), V3=c(4, 1)) tf - ifelse(example0, TRUE, FALSE) example[tf] - NA apply(example, 1, mean, na.rm=TRUE) -- David L Carlson Associate Professor

Re: [R] How to apply a function to a multidimensional array based on its indices

2012-05-14 Thread David Winsemius
On May 14, 2012, at 10:09 AM, math_daddy wrote: Hello. I have a 4 dimensional array and I want to fill in the slots with values which are a function of the inputs. Through searching the forums here I found that the function outer is helpful for 2x2 matrices but cannot be applied to

Re: [R] Discrete choice model maximum likelihood estimation

2012-05-14 Thread Berend Hasselman
See below. On 14-05-2012, at 13:21, infinitehorizon wrote: Hello again, I changed the name to tt. and for a and tt actually I was getting them from data, I didn't put them here in the question. Now I restructured my code and below I copy the full code, I tried many things but still

Re: [R] New to R

2012-05-14 Thread Sarah Goslee
Hi, On Mon, May 14, 2012 at 12:44 PM, Bert Gunter gunter.ber...@gene.com wrote: I usually try google searches first. While not always successful, I am frequently surprised by how well it does. rseek.org is simply a custom Google search for R-related things. It does an even better job pulling

Re: [R] rgl: cylinder3d() with elliptical cross-section

2012-05-14 Thread Duncan Murdoch
On 10/03/2012 9:44 AM, Duncan Murdoch wrote: My first reply to this went privately, by accident. I've done a little editing to it, but mainly this is for the archives. On 12-03-09 2:36 PM, Michael Friendly wrote: For a paper dealing with generalized ellipsoids, I want to illustrate in 3D

Re: [R] as.function parameters

2012-05-14 Thread jackl
Hi, ~ well that seems to be a better solution. Question is how much an enviroment for each node costs in terms of save space.. The example code is hard to present, because it is really based on that problem. The frame of the problem is, that I have to write a program that gives each node in a

Re: [R] range segment exclusion using range endpoints

2012-05-14 Thread Ben quant
Turns out this solution doesn't work if the s range is outside the range of the x range. I didn't include that in my examples, but it is something I have to deal with quite often. For example s1_rng below causes an issue: x_rng = c(-100,100) s1_rng = c(-250.5,30) s2_rng = c(0.77,10) s3_rng =

Re: [R] as.function parameters

2012-05-14 Thread David Winsemius
On May 14, 2012, at 1:56 PM, jackl wrote: Hi, ~ well that seems to be a better solution. Question is how much an enviroment for each node costs in terms of save space.. Seems unlike that it would expand you space very much. Every function you create will have an environment. The

Re: [R] range segment exclusion using range endpoints

2012-05-14 Thread William Dunlap
To the list of function I sent, add another that converts a list of intervals into a Ranges object: as.Ranges.list - function (x, ...) { stopifnot(nargs() == 1, all(vapply(x, length, 0) == 2)) # use c() instead of unlist() because c() doesn't mangle POSIXct and Date objects x

Re: [R] range segment exclusion using range endpoints

2012-05-14 Thread Ben quant
Yes, it is. I'm looking into understanding this now... thanks! Ben On Mon, May 14, 2012 at 12:38 PM, William Dunlap wdun...@tibco.com wrote: To the list of function I sent, add another that converts a list of intervals into a Ranges object: as.Ranges.list - function (x, ...) {

Re: [R] Random forests prediction

2012-05-14 Thread Liaw, Andy
That's not how RF works at all. The setting of mtry is irrelevant to this. Andy -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of matt Sent: Monday, May 14, 2012 10:22 AM To: r-help@r-project.org Subject: Re: [R] Random forests

Re: [R] range segment exclusion using range endpoints

2012-05-14 Thread Steve Lianoglou
Hi all, Nice code samples presented all around. Just wanted to point out that I think the stuff found in the `intervals` package might also be helpful: http://cran.at.r-project.org/web/packages/intervals/index.html HTH, -steve On Mon, May 14, 2012 at 2:54 PM, Ben quant ccqu...@gmail.com

[R] Inf and lazy evaluation

2012-05-14 Thread J Toll
Hi, I have a question involving Inf, lazy evaluation, and maybe argument recycling. I have a directory where I am trying to check for the existence of files of a certain pattern, basically something like filename-#, where # is an integer. I can do something like this, which works.

Re: [R] Inf and lazy evaluation

2012-05-14 Thread Bert Gunter
?list.files -- Bert On Mon, May 14, 2012 at 12:34 PM, J Toll jct...@gmail.com wrote: Hi, I have a question involving Inf, lazy evaluation, and maybe argument recycling.  I have a directory where I am trying to check for the existence of files of a certain pattern, basically something like

Re: [R] Inf and lazy evaluation

2012-05-14 Thread R. Michael Weylandt
R is lazy, but not quite that lazy ;-) It's likely much easier to do this with regexps something like list.files()[grepl(paste0(filename, -[0123456789]+), list.files())] Michael On Mon, May 14, 2012 at 3:34 PM, J Toll jct...@gmail.com wrote: Hi, I have a question involving Inf, lazy

Re: [R] range segment exclusion using range endpoints

2012-05-14 Thread Ben quant
Thank you Steve! This does everything I need (at this point): (this excludes ranges y2 from range y1) library('intervals') y1 = Intervals(c(-100,100)) y2 = Intervals(rbind( c(-100.5,30), c(0.77,10), c(25,35), c(70,80.3), c(90,95) )) interval_difference(y1,y2) Object of class

[R] Plot

2012-05-14 Thread Kelly Cool
Hello, I am trying to make a plot of the rates of an enzyme against three different protein concentrations (there are 45 rates in total and split up into 3 groups of 15, each receiving one of the 3 protein concentrations). When I enter the following code I instead get 3 separate boxplots for

[R] Error in names(x) - value: 'names' attribute must be the same length as the vector

2012-05-14 Thread Priya Bhatt
Dear R-helpers, I am stuck on an error in R: When I run my code (below), I get this error back: Error in names(x) - value : 'names' attribute must be the same length as the vector Then when I use traceback(), R gives me back this in return: `colnames-`(`*tmp*`, value = c(Item, Color

Re: [R] Help with V function in igraph

2012-05-14 Thread bmccowan
Thank for your response. It is oddly working now. Thanks again, Brenda On 5/14/12 10:36 AM, Gábor Csárdi-2 [via R] ml-node+s789695n4629973...@n4.nabble.com wrote: Something weird must be going on in your s641_social object. Can you just simply check that the vertex names look OK with

Re: [R] range segment exclusion using range endpoints

2012-05-14 Thread William Dunlap
Yes, package:intervals uses the same idiom as my code. Mine allows Date and POSIXct objects as interval endpoints (which is why it represents the objects as data.frames instead of matrices). package:intervals has more than basic set opertions on the collections of intervals. Bill Dunlap

Re: [R] Inf and lazy evaluation

2012-05-14 Thread J Toll
Thank you all for the replies. On Mon, May 14, 2012 at 2:45 PM, R. Michael Weylandt michael.weyla...@gmail.com wrote: R is lazy, but not quite that lazy ;-) Oh, what is this world coming to when you can't count on laziness to be lazy. ;) I should probably stop reading about Haskell and their

  1   2   >