similar to: Probable Numeric Bugs in svd function (PR#8781)

I am new to R and have written a function that clusters on subsets of a big data data set with 60,000 points. I am not sure why, but I keep getting a run-time error. Any suggestions would be greatly appreciated. Here is the code: library(cba) d<-read.csv("data.csv", header=TRUE) v<-c(53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,21,72 73,74,75,76,77,78)

Surprisingly, I played around with some test code and below actually creates equations that look correct. tempmat<-matrix(10,nrow=6,ncol=6) restrictmat<-diag(6) colnames(tempmat)<-c("AUD.l1","CHF.l1","CAD.l1","GBP.l1","EUR.l1","JPY.l 1")

Dear R users, I'm glad to announce the release of a new package on CRAN for exploring correlations in R: corrr <https://cran.rstudio.com/web/packages/corrr/index.html> . corrr provides a handful of methods to create and explore a correlation tbl (rather than a matrix). E.g., rearrange() the correlations based on their strength, or focus() on selected variables against others. All

Dear R users, I'm glad to announce the release of a new package on CRAN for exploring correlations in R: corrr <https://cran.rstudio.com/web/packages/corrr/index.html> . corrr provides a handful of methods to create and explore a correlation tbl (rather than a matrix). E.g., rearrange() the correlations based on their strength, or focus() on selected variables against others. All

I was wondering if there is a way of editting strings in R. I have a set of strings and each set is a row of numbers and paranthesis. For example the first row is: (0 2)(3 4)(7 9)(5 9)(1 5) and I have a thousand or so such rows. I was wondering how I could get the corresponding string obtained by adding 1 to all the numbers in the string above. Dursun [[alternative HTML version deleted]]

R-help - I have this set of aggregated tables (sample data below via dput()). And I would like to have delayValue as the column variables with the "temp" (temp1, temp2, temp3) values as the row variables. However I would like to have the temp variables *aggregated into single rows* so that I have the frequency ("Freq" | counts) of each time each "delayValue" occurs

Hi everyone, I'm a newbie to R and to linear mixed effects modeling, so please have mercy :-) Just wanted to check, whether what I'm doing is alright. I've collected data concerning tonotopic organization of the auditory cortex in humans, and I have approximately 1400-1800 data/time points per person (13 in total). Observations were made how the focus of neuronal activity changed

this seems l ike it shouldn't be that hard but i give up. if i have a string say, temp<-"01", I want to increase it by 1 so that it becomes "02". but the following code obviously won't work when the input string is say "10" because then it gives "011" when I just want "11". uuugh. does someone know an easy way ( i guess hard is

Please keep the conversation on the list: Others may be able to help you better than I can. On 2017-07-06 10:38, SEB140004 Student wrote: > Ya. I had successfully got the result. Thank you very much. :) > > Isn't possible for me to obtain the network from the correlation matrix? I know nothing about correlation networks, but by googling I found two R packages, ggraph and corrr,

Hi Ralph, I applied the patch to a new build of Sernet 4.4.6 SRPM. I could still open the file from both hosts, and in fact this time I could actually open it on the same host without any problems: host 1: smbclient -U ajc //172.31.0.120/ifa_v ... smb: \jmi\Corporate Accounting\ISL Ltd Accounts\2016\1609\> open "201609 213_ ISL Balance sheet variance analysis.xls" open file

Does anyone know how -log(x) can equal 743 but -log(x+0)=Inf? That's what the following stream of calculations suggest: Browse[2]> -log ( 1e-323+yMat2 - yMat1 * logitShape(matrix(parsList$Xs, nrow = numXs, ncol=numOfCurves), matrix(means, nrow = numXs, ncol=numOfCurves, byrow=TRUE), matrix(sigmas, nrow = numXs, ncol=numOfCurves, byrow=TRUE)) )[5,9] [1] Inf Yet: Browse[2]>

Hi, While playing with quantile-quantile plots, I wrote up some code which plots something strangely different. Here's the pseudocode: testhist <- hist(sample_data) refhist <- hist(rnorm(n, mean=0,sd=1)) # for some large-ish n cumtest <- cumsum(testhist) cumref <- cumsum(refhist) plot(cumref,cumtest) This produces a straight line of slope 1 for a sample with the same

Alex, On Sat, Oct 22, 2016 at 07:30:01PM +0100, Alex Crow wrote: > I applied the patch to a new build of Sernet 4.4.6 SRPM. I could still > open the file from both hosts, and in fact this time I could actually > open it on the same host without any problems: I don't see where you're opening the file twice from one smbclient below. Am I missing something? And then there's no

Hi I am trying to estimate bottom temperatures over a particular depth range, based on one dataset containing synthetic temperature profiles (for set depths) and another that contains information on bathymetry. I need to do this for multiple regions and thus would ideally like an automated solution that is quicker and more accurate than the manual option I currently have! I have not had much luck

Dear R experts: I just tried some simple test that told me that hand computing the OLS coefficients is about 3-10 times as fast as using the built-in lm() function. (code included below.) Most of the time, I do not care, because I like the convenience, and I presume some of the time goes into saving a lot of stuff that I may or may not need. But when I do want to learn the properties of an

Dear List, I would like to plot a simple legend with two math expressions, e.g. plot(0) legend(1, 0.5, expression(sigma[i], sigma[j])) The difficulty is that i and j should be variables rather than strings "i" and "j". In other words I'd like to do something like: i = "A" j = "B" legend(1, 0.5, expression(sigma[i], sigma[j])) and have "A"

I'm using R on Windows and wondering if there is an R package that allows you to do folder structure analysis.   Specifically, I would like to grab the name, file size, type, and creation data of files within a folder and then to means and sigmas on the file size and creationdation date.   Thank you again for your insights.     [[alternative HTML version deleted]]

R listers: I have been foiled by plotmath! (in R 2.01,Windows 2000) The task: Plot a normal density and label the ticks as mu - 3 sigma, mu - 2 sigma, ...., mu + 3 sigma, where the mu's and sigmas appear as Greek symbols, of course. The following code does this: x<-seq(-3,to=3,by=.01) y<-dnorm(x) plot(x,y,type='h',col='lightblue',axes=FALSE)

I tried to make 5 separate patches, but as they are constructive, they had information from the last one. So I will post one by one as it gets on the tree. 1. We had a function isUsedInTerminator that tested if a comparator was used in the terminator of its parent BasicBlock. This is wrong because a comparator can be created in a BasicBlock and used in the terminator of other BasicBlock, and

Dear R-users, I have longitudinal data and would like to fit a model where both the variance-covariance matrix of the random effects and the residual variance are conditional on a (binary) grouping variable. I guess the model would have the following form (in hierarchical notation) Yi|bi,k ~ N(XiB+Zibi, sigmak*Ident) bi|k ~ N(0, Dk) K~Bernoulli(p) I can obtain different sigmas (sigma0 and