Re: [R] solving simultaneous Equations in R

2013-12-12 Thread Berend Hasselman
On 11-12-2013, at 23:56, eliza botto eliza_bo...@hotmail.com wrote: Dear Berend, Thankyou very much indeed for you reply. By taking help from your previous reply @ http://r.789695.n4.nabble.com/Simultaneous-equations-td2524645.html i was able to generate the following loop for the

Re: [R] refline in forest() {metafor}

2013-12-12 Thread Viechtbauer Wolfgang (STAT)
Regarding your question: Do you want *another* line or do you just want to move the reference line to the value of the summary estimate? The latter can be done by passing the value of the summary estimate to the 'refline' argument. If you want another line, you could just use the abline

Re: [R] refline in forest() {metafor}

2013-12-12 Thread johnwilliams
I would like another line, and the solution below will work great. Thank you. - John Quoting Viechtbauer Wolfgang (STAT) wolfgang.viechtba...@maastrichtuniversity.nl: Regarding your question: Do you want *another* line or do you just want to move the reference line to the value of the summary

[R] Logistic Regression with 200K features in R?

2013-12-12 Thread Romeo Kienzler
Dear List, I'm quite new to R and want to do logistic regression with a 200K feature data set (around 150 training examples). I'm aware that I should use Naive Bayes but I have a more general question about the capability of R handling very high dimensional data. Please consider the

Re: [R] Logistic Regression with 200K features in R?

2013-12-12 Thread Eik Vettorazzi
it is simply because you can't do a regression with more predictors than observations. Cheers. Am 12.12.2013 09:00, schrieb Romeo Kienzler: Dear List, I'm quite new to R and want to do logistic regression with a 200K feature data set (around 150 training examples). I'm aware that I

Re: [R] solving simultaneous Equations in R

2013-12-12 Thread eliza botto
Dear Berend,Got it! Thankyou very much indeed for your kind support and elaborating the issue.I'm Extremely grateful!!! Eliza Subject: Re: [R] solving simultaneous Equations in R From: b...@xs4all.nl Date: Thu, 12 Dec 2013 09:01:02 +0100 CC: r-help@r-project.org To: eliza_bo...@hotmail.com

Re: [R] Logistic Regression with 200K features in R?

2013-12-12 Thread Eik Vettorazzi
I thought so (with all the limitations due to collinearity and so on), but actually there is a limit for the maximum size of an array which is independent of your memory size and is due to the way arrays are indexed. You can't create an object with more than 2^31-1 = 2147483647 elements.

Re: [R] Logistic Regression with 200K features in R?

2013-12-12 Thread Duncan Murdoch
On 13-12-12 6:51 AM, Eik Vettorazzi wrote: I thought so (with all the limitations due to collinearity and so on), but actually there is a limit for the maximum size of an array which is independent of your memory size and is due to the way arrays are indexed. You can't create an object with more

Re: [R] Logistic Regression with 200K features in R?

2013-12-12 Thread Eik Vettorazzi
thanks Duncan for this clarification. A double precision matrix with 2e11 elements (as the op wanted) would need about 1.5 TB memory, that's more than a standard (windows 64bit) computer can handle. Cheers. Am 12.12.2013 13:00, schrieb Duncan Murdoch: On 13-12-12 6:51 AM, Eik Vettorazzi wrote:

[R] Solving a normal distribution pnorm for q

2013-12-12 Thread Johannes Radinger
Hi, I found follwowing example of pnorm here: http://www.r-tutor.com/elementary-statistics/probability-distributions/normal-distribution Problem Assume that the test scores of a college entrance exam fits a normal distribution. Furthermore, the mean test score is 72, and the standard

Re: [R] Solving a normal distribution pnorm for q

2013-12-12 Thread Berend Hasselman
On 12-12-2013, at 14:56, Johannes Radinger johannesradin...@gmail.com wrote: Hi, I found follwowing example of pnorm here: http://www.r-tutor.com/elementary-statistics/probability-distributions/normal-distribution Problem Assume that the test scores of a college entrance

Re: [R] Solving a normal distribution pnorm for q

2013-12-12 Thread Dániel Kehl
Looks like homework. Try ?qnorm Feladó: r-help-boun...@r-project.org [r-help-boun...@r-project.org] ; meghatalmaz#243;: Johannes Radinger [johannesradin...@gmail.com] Küldve: 2013. december 12. 14:56 To: R help Tárgy: [R] Solving a normal distribution

Re: [R] Solving a normal distribution pnorm for q

2013-12-12 Thread Jeff Newmiller
?pnorm ... carefully... --- Jeff NewmillerThe . . Go Live... DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live Go... Live:

[R] labels on right y-axis

2013-12-12 Thread Fisher Dennis
OS X R 3.0.2 Colleagues, I am displaying two sets of values changing over time; as a result, I have two y-axes. I add a label for the right-side Y axis with mtext(side=3, line=1.2, TEXT). The text is parallel to the axis -- so far, so good. However, the text is rotated counterclockwise

Re: [R] Solving a normal distribution pnorm for q

2013-12-12 Thread Johannes Radinger
Thanks for the fast response. Of course I totally overlooked qnorm as I had a more complex task in my head. I wanted to reverse following equation: r=0.67 q=-150 sd1=100 sd2=1000 X - pnorm(q=q,sd=sd1)*r+pnorm(q=q,sd=sd2)*(1-r) Maybe its mathematically really easy, but somehow I don't get it how

Re: [R] radial.plot shaded region

2013-12-12 Thread Hodge, Steven
From: Jim Lemon [j...@bitwrit.com.au] Sent: Tuesday, December 10, 2013 21:57 To: Hodge, Steven Cc: r-help@r-project.org Subject: Re: [R] radial.plot shaded region On 12/11/2013 11:24 AM, Jim Lemon wrote: ... There may be a solution using the radial.pie function, and if I find it, I'll post it.

Re: [R] Solving a normal distribution pnorm for q

2013-12-12 Thread Eik Vettorazzi
Assuming r, sd1 and sd2 as known(fixed) components, this works: X - pnorm(q=q,sd=sd1)*r+pnorm(q=q,sd=sd2)*(1-r) uniroot(function(x)pnorm(q=x,sd=sd1)*r+pnorm(q=x,sd=sd2)*(1-r)-X,c(-1e6,1e6)) Cheers Am 12.12.2013 15:58, schrieb Johannes Radinger: Thanks for the fast response. Of course I

[R] censored counts and glmer/glmmADMB

2013-12-12 Thread Eric Elguero
dear R-users, I have to model counts where all counts above some threshold have been censored. In the same dataset I have too many zeroes for a Poisson or even a negative binomial distribution to make sense, so I would need a zero-inflated-censored negative binomial family for use in glmer (or

Re: [R] Solving a normal distribution pnorm for q

2013-12-12 Thread Johannes Radinger
Thanks Eik... ..the uniroot() was the function I was looking for. Best regards, /J On Thu, Dec 12, 2013 at 4:27 PM, Eik Vettorazzi e.vettora...@uke.de wrote: Assuming r, sd1 and sd2 as known(fixed) components, this works: X - pnorm(q=q,sd=sd1)*r+pnorm(q=q,sd=sd2)*(1-r)

Re: [R] censored counts and glmer/glmmADMB

2013-12-12 Thread Bert Gunter
This post should go the r-sig-mixed-models list, where you are much more likely to get useful help than here, which is a general R programming help list. Cheers, Bert On Thu, Dec 12, 2013 at 7:41 AM, Eric Elguero eric.elgu...@ird.fr wrote: dear R-users, I have to model counts where all counts

Re: [R] labels on right y-axis

2013-12-12 Thread David Carlson
I don't see an alternative to text(), but the positioning is not that difficult: oldpar - par(mar=c(5.1, 4.1, 4.1, 4.1)) plot(0) axis(4) text(par(usr)[2], mean(par(usr)[3:4]), TEXT, srt=-90, adj=c(.5,-4), xpd=TRUE) par(oldpar) - David L Carlson

Re: [R] Logistic Regression with 200K features in R?

2013-12-12 Thread Duncan Murdoch
On 12/12/2013 7:08 AM, Eik Vettorazzi wrote: thanks Duncan for this clarification. A double precision matrix with 2e11 elements (as the op wanted) would need about 1.5 TB memory, that's more than a standard (windows 64bit) computer can handle. According to Microsoft's Memory Limits web page

[R] function inside a function

2013-12-12 Thread eliza botto
Dear users of R, I am trying to inculcate a function inside a function. For that to be done, i copied following function from internet. library(nleqslv) fun - function(x) { f - numeric(length(x)) f[1] - A[,1]+x[2] - 1/x[1] f[2] - A[,2]+x[2] - sin(x[1]) f } x.start -

Re: [R] labels on right y-axis

2013-12-12 Thread Bert Gunter
Might axis() be useful here? ?axis Cheers, Bert On Thu, Dec 12, 2013 at 8:02 AM, David Carlson dcarl...@tamu.edu wrote: I don't see an alternative to text(), but the positioning is not that difficult: oldpar - par(mar=c(5.1, 4.1, 4.1, 4.1)) plot(0) axis(4) text(par(usr)[2],

[R] Heatmap, and heatmap.2 gave different figures for the same dataset

2013-12-12 Thread capricy gao
I have a huge dataset(15k X 18) and tried to use the heatmap in R to examine the patterns. However, I found that heatmap and heatmap.2 gave me completely different outputs. Here are the codes: dim(as.matrix(data.dcpm)) [1] 15462    18 heatmap(as.matrix(data.dcpm),

Re: [R] Logistic Regression with 200K features in R?

2013-12-12 Thread Romeo Kienzler
ok, so 200K predictors an 10M observations would work? On 12/12/2013 12:12 PM, Eik Vettorazzi wrote: it is simply because you can't do a regression with more predictors than observations. Cheers. Am 12.12.2013 09:00, schrieb Romeo Kienzler: Dear List, I'm quite new to R and want to do

[R] boxcox transformations

2013-12-12 Thread Yogendra
Hi, I am new to R. I need help with regards to box cox transformation. I have phenotypic data for e.g. plant height. data is non-normal. Skewness is 0.34. Could you please help me? Regards, Yogi -- View this message in context:

Re: [R] Logistic Regression with 200K features in R?

2013-12-12 Thread Romeo Kienzler
Dear Eik, thank you so much for your help! best Regards, Romeo On 12/12/2013 12:51 PM, Eik Vettorazzi wrote: I thought so (with all the limitations due to collinearity and so on), but actually there is a limit for the maximum size of an array which is independent of your memory size and is

Re: [R] Wilcoxon test output as a table

2013-12-12 Thread dc347
Hello there Was looking up a similar problem yesterday and came across the gMWT package (Generalized Mann-Whitney Type Tests). It addresses your problem, although I am still trying to figure out how it works. Kind regards, David -- View this message in context:

Re: [R] function inside a function

2013-12-12 Thread Berend Hasselman
On 12-12-2013, at 17:10, eliza botto eliza_bo...@hotmail.com wrote: Dear users of R, I am trying to inculcate a function inside a function. For that to be done, i copied following function from internet. library(nleqslv) fun - function(x) { f - numeric(length(x)) f[1] -

Re: [R] function inside a function

2013-12-12 Thread eliza botto
Dear Berend,Thankyou very much for your reply. I actually wanted to inserts each row value of A in f[1] anf f[2], with column 1 value in f[1] and column 2 values in f[2]. Once that been done, i should have in the end, 124 vaues of x[1] for column 1 and 124 for column 2.I hope i am clear this

Re: [R] boxcox transformations

2013-12-12 Thread Dániel Kehl
Dear Yogi, try the powerTransform function in the car package and boxcox function in the MASS package. there are nice examples for both functions. HTH, daniel Feladó: r-help-boun...@r-project.org [r-help-boun...@r-project.org] ; meghatalmaz#243;:

Re: [R] Heatmap, and heatmap.2 gave different figures for the same dataset

2013-12-12 Thread Kevin Wright
Read the help page for both and pay particular attention to the scale argument. Kevin Wright On Thu, Dec 12, 2013 at 9:45 AM, capricy gao capri...@yahoo.com wrote: I have a huge dataset(15k X 18) and tried to use the heatmap in R to examine the patterns. However, I found that heatmap and

Re: [R] labels on right y-axis

2013-12-12 Thread Jim Lemon
On 12/13/2013 01:32 AM, Fisher Dennis wrote: OS X R 3.0.2 Colleagues, I am displaying two sets of values changing over time; as a result, I have two y-axes. I add a label for the right-side Y axis with mtext(side=3, line=1.2, TEXT). The text is parallel to the axis -- so far, so good.

Re: [R] function inside a function

2013-12-12 Thread eliza botto
Dear Arun, Thankyou very much for your reply. You always come up to save the day. So kind of you...!!! Eliza Date: Thu, 12 Dec 2013 12:36:18 -0800 From: smartpink...@yahoo.com Subject: Re: [R] function inside a function To: eliza_bo...@hotmail.com Hi Eliza, Actually, I used the

[R] method default for hclust function

2013-12-12 Thread capricy gao
I could not figure out what was the default when I ran hclust() without specifying the method. For example: I just have a code like: hclust(dist(data)) Any input would be appreciated:) [[alternative HTML version deleted]] __

Re: [R] method default for hclust function

2013-12-12 Thread eliza botto
Absolute distance is the default distance in hclust. v-c(1,2,3,4,5,6) dist(v) 2 1 3 2 1 4 3 2 1 5 4 3 2 1 6 5 4 3 2 1 Eliza Date: Thu, 12 Dec 2013 15:09:19 -0800 From: capri...@yahoo.com To: r-help@r-project.org Subject: [R] method default for hclust function I could

Re: [R] method default for hclust function

2013-12-12 Thread Peter Langfelder
On Thu, Dec 12, 2013 at 3:09 PM, capricy gao capri...@yahoo.com wrote: I could not figure out what was the default when I ran hclust() without specifying the method. According to help(hclust), the default method is complete linkage. HTH, Peter __

Re: [R] Heatmap, and heatmap.2 gave different figures for the same dataset

2013-12-12 Thread capricy gao
Thank you very much for the hints. I fixed the problem. On Thursday, December 12, 2013 2:08 PM, Kevin Wright kw.s...@gmail.com wrote: Read the help page for both and pay particular attention to the scale argument. Kevin Wright On Thu, Dec 12, 2013 at 9:45 AM, capricy gao

[R] Coxme time-dependent covariate

2013-12-12 Thread Jared Duquette
Hi all, I am using the coxme function to fit random effects survival models. The base survival models show a typical J shaped curve where survival drops sharply then plateaus over time interval. I want to include a time covariate which could account for this non-linear decrease in survival

Re: [R] function inside a function

2013-12-12 Thread arun
Hi, May be this helps: t(sapply(seq_len(nrow(A)),function(i) nleqslv(x.start,fun,krow=i)$x)) A.K. On Thursday, December 12, 2013 1:54 PM, eliza botto eliza_bo...@hotmail.com wrote: Dear Berend,Thankyou very much for your reply. I actually wanted to inserts each row value of A in f[1] anf

[R] charToRaw(Œ) is not 8C in R console

2013-12-12 Thread 水静流深
in http://www.ascii-code.com/, you can see the the hex value of Œ is 8C, why in my R console ? charToRaw(Œ) [1] c5 92 is not 8C ? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list