[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 https://stat.ethz.ch/mailman/l

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 wrote: 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 i

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

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 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 wrote: I have a huge datas

Re: [R] method default for hclust function

2013-12-12 Thread Peter Langfelder
On Thu, Dec 12, 2013 at 3:09 PM, capricy gao 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 __ R-help@r-projec

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 c

[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]] __ R-help@r-project

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

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

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 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 heatmap.2 gave me

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ó: Yogendra [y

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 tim

Re: [R] function inside a function

2013-12-12 Thread Berend Hasselman
On 12-12-2013, at 17:10, eliza botto 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] <- A[,1]+x[2] - 1/x[1

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: http://r.789695.n4.nabble.com/

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 du

[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: http://r.789695.n4.nabble.com/boxcox-transformations-tp4682077.h

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 logis

[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), col=topo.

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 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], mean(par("usr")[3:4])

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

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 Dep

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 wrote: > dear R-users, > > I have to model counts where all counts above some thres

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 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) > > uniroot(function(x)pnorm(q=x,sd=s

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

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 total

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 pos

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

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

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

2013-12-12 Thread Jeff Newmiller
?pnorm ... carefully... --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead: OO#.. Playing

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ó: Johannes Radinger [johannesradin...@gmail.com] Küldve: 2013. december 12. 14:56 To: R help Tárgy: [R] Solving a normal distribution pnorm fo

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

2013-12-12 Thread Berend Hasselman
On 12-12-2013, at 14:56, Johannes Radinger 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 exam fits a norma

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

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:

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
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. https://s

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

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

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

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)" : > 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 pa

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 function

Re: [R] solving simultaneous Equations in R

2013-12-12 Thread Berend Hasselman
On 11-12-2013, at 23:56, eliza botto 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 calculation of x=x[1] and