[julia-users] How quickly subtract the two large arrays ?

2014-08-15 Thread paul analyst
How quickly subtract the two large arrays ? A = random (10 ^ 5.10 ^ 5); B = rand (10 ^ 5.10 ^ 5); @time C = A-B ; very long: / Paul

[julia-users] How to save sparse in HDF5? ERROR: no method write(HDF5File, ASCIIString, SparseMatrixCSC{Int64,Int64})

2014-08-14 Thread paul analyst
How to save sparse in HDF5? julia> D=sparse(I,J,V); julia> h5write("D.h5","D", D); ERROR: no method write(HDF5File, ASCIIString, SparseMatrixCSC{Int64,Int64}) in h5write at C:\Users\SAMSUNG2\.julia\v0.3\HDF5\src\plain.jl:527 Paul

[julia-users] What's wrong? A simple condition is not working. [L. == 1] ERROR: BoundsError() in getindex_bool_1d at array.jl:285

2014-08-13 Thread paul analyst
F is dense julia> size(F) (6237437,284) julia> size(L) (6237437,1) julia> F[:,1][L.==0] ERROR: BoundsError() in getindex_bool_1d at array.jl:285 julia> L 6237437x1 sparse matrix with 20869 Int16 entries: [66 , 1] = 1 [104, 1] = 0 Paul

Re: [julia-users] How to divide DArray vertically and not horizontally?

2014-08-10 Thread Paul Analyst
.. and beter is distributed set by columns.. Is any idea ? Paul W dniu 2014-08-10 17:20, Kevin Squire pisze: Actually, Julia uses column major ordering of arrays. Cheers, Kevin On Sunday, August 10, 2014, Paul Analyst <mailto:paul.anal...@mail.com>> wrote: Thx Gael Is it on

Re: [julia-users] How to divide DArray vertically and not horizontally?

2014-08-10 Thread Paul Analyst
Thx Gael Is it one way ? I have next proceses optimized for cols ... Paul W dniu 2014-08-10 15:32, gael.mc...@gmail.com pisze: From my understanding, Julia being row major, it makes little sense to split your arrays vertically. If your algorithm makes use of a complete row instead of a complet

Re: [julia-users] simply parallel , please help

2014-08-10 Thread paul analyst
; > > On Wed, Aug 6, 2014 at 5:39 AM, paul analyst > wrote: > >> How to count consecutive parallel loops at work on a large matrix? >> D = big sparse >> for i = 0: ... >> c$i = cor (D [:, and * 1000 + 1: and * 1000 + 1000]) >> end >> Paul >> > >

[julia-users] DArray. What is wrong ? ERROR: no method +(DArray{Float64,2,Array{Float64,2}}, Array{Int64,1}) in sum at tuple.jl:118

2014-08-10 Thread paul analyst
What is wrong? Why "dd" array is of type "Any" ? ulia> addprocs(4) 4-element Array{Any,1}: 2 3 4 5 julia> D=readcsv("v4dane.txt"); julia> d=distribute(D) 2888x4 DArray{Float64,2,Array{Float64,2}}: julia> d 2888x4 DArray{Float64,2,Array{Float64,2}}: 0.0 0.0 0.00.0 0.0 50.0

[julia-users] How to divide DArray vertically and not horizontally?

2014-08-10 Thread paul analyst
How to divide DArray vertically and not horizontally? I need 2 columns with 512 rows ... Paul julia> addprocs(2) 2-element Array{Any,1}: 2 3 julia> julia> remotecall_fetch(2, whos) From worker 2: Base Module From worker 2: Core

[julia-users] simply parallel , please help

2014-08-06 Thread paul analyst
How to count consecutive parallel loops at work on a large matrix? D = big sparse for i = 0: ... c$i = cor (D [:, and * 1000 + 1: and * 1000 + 1000]) end Paul

[julia-users] Re: What wrong ? Df[:,i]=D[:,i].-E

2014-07-31 Thread paul analyst
No help ? W dniu czwartek, 31 lipca 2014 18:12:24 UTC+2 użytkownik paul analyst napisał: > > i=1 > Df=zeros(D); # sparse array 10^6 x 10^3 > ee=mean(D[:,i]); > Y=(D[:,i].-ee); > temp=mean(Y) > E=ee-temp; > Df[:,i]=D[:,i].-E > ERROR: InexactError() > in copy!

[julia-users] What wrong ? Df[:,i]=D[:,i].-E

2014-07-31 Thread paul analyst
i=1 Df=zeros(D); # sparse array 10^6 x 10^3 ee=mean(D[:,i]); Y=(D[:,i].-ee); temp=mean(Y) E=ee-temp; Df[:,i]=D[:,i].-E ERROR: InexactError() in copy! at abstractarray.jl:149 in convert at sparse/sparsematrix.jl:160 in vcat at abstractarray.jl:503 Paul

[julia-users] Re: Correlation sparse array is very slow.

2014-07-28 Thread paul analyst
day, July 28, 2014 1:41:55 PM UTC+5:30, paul analyst wrote: >> >> Correlation sparse array is very slow. Out of memory on a dense array when >> we have 30,000 columns. How quickly it calculated? >> >> julia> I=int32((rand(10^7)*999).+1); >> >> julia

[julia-users] Correlation sparse array is very slow.

2014-07-28 Thread paul analyst
Correlation sparse array is very slow. Out of memory on a dense array when we have 30,000 columns. How quickly it calculated? julia> I=int32((rand(10^7)*999).+1); julia> J=int32((rand(10^7)*2).+1); julia> V=int8((rand(10^7)*9).+1); julia> D=sparse(I,J,V); julia> @time cor(D[:,1:30]);

[julia-users] How to save sparse array to csv

2014-07-27 Thread paul analyst
How to save sparse array to csv [1, 4] = 1 [4, 7] = 2 [5, 9] = 3 [3, 18] = -5 ? Paul

[julia-users] Re: How to read data from a file to a new sparse matrix?

2014-07-27 Thread paul analyst
rse(I,J,V)`, where I is the vector of row indices, J of > column indices, and V of values. > > -viral > > On Friday, July 25, 2014 1:29:07 AM UTC+5:30, paul analyst wrote: >> >> how to read data from a file to a new sparse matrix? date = readcsv (" >> data.txt") >> Paul >> >>

[julia-users] How to read data from a file to a new sparse matrix?

2014-07-24 Thread paul analyst
how to read data from a file to a new sparse matrix? date = readcsv (" data.txt") Paul

[julia-users] How cartesianmap take into variable ?

2014-07-13 Thread paul analyst
Julia Documentation: cartesianmap(f, dims) Given a dims tuple of integers (m, n, ...), call f on all combinations of integers in the ranges 1:m, 1:n, etc. Example: julia> cartesianmap(println, (2,2)) 11 21 12 22 julia>x= cartesianmap(2,2) ERROR: no method cartesianmap(Int64, Int64) How to take

Re: [julia-users] Sort error, How to find reserved word in big text ?

2014-07-13 Thread Paul Analyst
THX, it work! Paul W dniu 2014-07-13 09:31, Stefan Karpinski pisze: readcsv("data.txt", UTF8String) should help. On Sun, Jul 13, 2014 at 12:15 AM, Paul Analyst <mailto:paul.anal...@mail.com>> wrote: ok, now : julia> data=readcsv("data.txt")

Re: [julia-users] Sort error, How to find reserved word in big text ?

2014-07-13 Thread Paul Analyst
ot;. How you defend against such situations? Paul W dniu 2014-07-13 09:03, John Myles White pisze: parse isn’t a cure for anything. parse interprets inputs under the assumption that they are valid Julia code. Your data isn’t valid Julia code, so you should not use parse. — John On Jul 13,

Re: [julia-users] Sort error, How to find reserved word in big text ?

2014-07-13 Thread Paul Analyst
At first I did and it was a problem with the "Inf". I realized that parse a cure for "Inf" and the like. Paul W dniu 2014-07-13 08:57, Stefan Karpinski pisze: You should NOT use parse to parse your data. Is that how you're doing it? On Sat, Jul 12, 2014

Re: [julia-users] Sort error, How to find reserved word in big text ?

2014-07-12 Thread Paul Analyst
aries from the network. Paul W dniu sobota, 12 lipca 2014 18:14:39 UTC+2 użytkownik Stefan Karpinski napisał: The problem here seems to be that you're comparing types of values that are incomparable like a float and a string.

[julia-users] how use find ?

2014-07-12 Thread paul analyst
I need all indexes where a[:] >5 julia> a=rand(10)*10 10-element Array{Float64,1}: 4.84005 8.29994 8.8531 3.42319 2.60318 7.25313 0.816263 4.44463 6.71836 4.65337 julia> a.>5 10-element BitArray{1}: false true true false false true false false true false julia> find(a.>5

[julia-users] Sort error, How to find reserved word in big text ?

2014-07-12 Thread paul analyst
I have a very long vector (5*10^6 ) with the words in a foreign language. I can not sort julia> sort (y1u) ERROR: no method isless (Float64, substring {UTF8String}) in the sort! at sort.jl: 233 in the sort! at sort.jl: 277 in next at string.jl: 630 Probably there is a word reserved f

[julia-users] Please help, I can't to deploy a simple parallelization .

2014-06-20 Thread paul analyst
Despite many attempts, I can't to deploy a simple parallelization and ask for help on a specific the example. Array WS can be broken down into clones, and at the finish may then add them, but WS is great and it would be better without cloning it. D=readcsv("D.txt"); k,l=size(D) 10,1 WS

Re: [julia-users] Re: How to convert vector to a series of possible pairs of elements ?

2014-06-20 Thread Paul Analyst
Vector[[i,j] for i in x, j in x] |> vec | On Thursday, June 19, 2014 12:16:32 PM UTC-7, paul analyst wrote: How to convert vector to a series of possible pairs of elements ? [1,2,3] I expect [1,1] [1,2] [1,3] [2,1] [2,2] [2,3] [3,1] [3,2] [3,3] ​ ​

[julia-users] How to convert vector to a series of possible pairs of elements ?

2014-06-19 Thread paul analyst
How to convert vector to a series of possible pairs of elements ? [1,2,3] I expect [1,1] [1,2] [1,3] [2,1] [2,2] [2,3] [3,1] [3,2] [3,3]

[julia-users] how to get vec with index (position) of a particular value (eg 1) in vector

2014-06-19 Thread paul analyst
How to Do a vector indicating the position of a particular value (eg 1) in another vector? by this question: (1)? [0 0 1 01] expected vector [3,5] Paul

Re: [julia-users] Re: How Delete row or col of array? like deleteat! in vectors

2014-06-18 Thread Paul Analyst
Great lesson thanks, I'll try :) Paul W dniu 2014-06-18 16:51, Tim Holy pisze: On Wednesday, June 18, 2014 06:06:24 AM paul analyst wrote: >Thanks, but I'm an analyst market, not a programmer;/ Still, I really like >Julia. Historically, only SPSS. I can so deeply programmed

Re: [julia-users] Re: How Delete row or col of array? like deleteat! in vectors

2014-06-18 Thread paul analyst
to vectors (if that > fits > your usual usage case), it would be possible to have a special version of > setdiff operating on ranges that would be much faster. Would be great if > you > could write & submit that. > > --Tim > > On Wednesday, June 18, 2014 02

Re: [julia-users] Re: How Delete row or col of array? like deleteat! in vectors

2014-06-18 Thread Paul Analyst
-18 13:38, Tim Holy pisze: setdiff --Tim On Wednesday, June 18, 2014 01:29:05 PM Paul Analyst wrote: number removal line is dynamic. This time 4 then maybe 2 maybe 3 etc Question now : How to construct fromvector [1:10] new vector without [2,4] like this [1,3,5:10] Paul W dniu 2014-06-18 1

Re: [julia-users] Re: How Delete row or col of array? like deleteat! in vectors

2014-06-18 Thread Paul Analyst
AM Paul Analyst wrote: Is possible or not ? Delete rows in array ? Paul W dniu 2014-06-17 07:49, Paul Analyst pisze: Thx, for this info But not about items I need delete some rows ... How fast delete rows in arrray, one [3,:]or more [[2,4],:] ? Paul W dniu 2014-06-16 21:22, Stefan Schwarz

Re: [julia-users] Re: How Delete row or col of array? like deleteat! in vectors

2014-06-18 Thread Paul Analyst
Is possible or not ? Delete rows in array ? Paul W dniu 2014-06-17 07:49, Paul Analyst pisze: Thx, for this info But not about items I need delete some rows ... How fast delete rows in arrray, one [3,:]or more [[2,4],:] ? Paul W dniu 2014-06-16 21:22, Stefan Schwarz pisze: deleting in place

Re: [julia-users] Re: How Delete row or col of array? like deleteat! in vectors

2014-06-16 Thread Paul Analyst
Thx, for this info But not about items I need delete some rows ... How fast delete rows in arrray, one [3,:]or more [[2,4],:] ? Paul W dniu 2014-06-16 21:22, Stefan Schwarz pisze: deleting in place a matrix is not supported and does not make sense. why is that? you've for instance a 5x5 matr

Re: [julia-users] Re: How to find the index of maximum but omitting to check any field eg not x [5,7]

2014-06-16 Thread Paul Analyst
Big thx, it work, Paul W dniu 2014-06-16 21:24, Dahua Lin pisze: |u =trues(length(x)) u[[5,7]]=false ir =0 vr =-Inf fori =1:length(x) ifu[i]&&x[i]>vr ir =i vr =x[i] end end |

[julia-users] Re: How to find the index of maximum but omitting to check any field eg not x [5,7]

2014-06-16 Thread paul analyst
No sugestion ? Paul W dniu poniedziałek, 16 czerwca 2014 14:06:23 UTC+2 użytkownik paul analyst napisał: > > I have a random vector x = rand (10) > how to find the index of maximum but omitting to check any field eg not x > [5,7] > some like: > indmax y = (x [but not read [5,7]]) > Paul >

[julia-users] How Delete row or col of array? like deleteat! in vectors

2014-06-16 Thread paul analyst
a=rand(5,5) how Delete row(s) or col(s) of array ? Paul

[julia-users] How to find the index of maximum but omitting to check any field eg not x [5,7]

2014-06-16 Thread paul analyst
I have a random vector x = rand (10) how to find the index of maximum but omitting to check any field eg not x [ 5,7] some like: indmax y = (x [but not read [5,7]]) Paul

[julia-users] Re: How use two contition to filtred array?

2014-06-10 Thread paul analyst
czerwca 2014 20:50:03 UTC+2 użytkownik Andrea Pagnani napisał: > > I think > > klient[ 500 .< klient[:,2] .<= 10,:] > > should work > > Andrea > > On Tuesday, June 10, 2014 8:41:17 PM UTC+2, paul analyst wrote: >> >> How use two conditions

[julia-users] How use two contition to filtred array?

2014-06-10 Thread paul analyst
How use two conditions to filtred array? Left and right side this code separatly are good. But how use booth conditions in filter ? (take rows with <1 and >500 in col 2) klient[klient[:,2].<=10;:] & klient[klient[:,2].>500;:] Paul

Re: [julia-users] Re: How to read and change the content of web pages to the vector ?

2014-06-09 Thread Paul Analyst
Thx James, i will try Paul. W dniu 2014-06-09 19:04, James Porter pisze: https://github.com/porterjamesj/GumboParser.jl

Re: [julia-users] Error Pkg.add("LightXML"), where is problem ?on my comp or on server?

2014-06-05 Thread Paul Analyst
the most usefull way:) Big Thx, working Paul W dniu 2014-06-05 17:33, Isaiah Norton pisze: Probably need to restart Julia. Otherwise, no idea. It works for me on Windows. On Thu, Jun 5, 2014 at 11:24 AM, Paul Analyst <mailto:paul.anal...@mail.com>> wrote: julia> Pkg.bui

Re: [julia-users] Error Pkg.add("LightXML"), where is problem ?on my comp or on server?

2014-06-05 Thread Paul Analyst
") ERROR: parse_file not defined file is in this dir... Paul W dniu 2014-06-05 17:21, Isaiah Norton pisze: try Pkg.build("LightXML") On Thu, Jun 5, 2014 at 11:03 AM, Paul Analyst <mailto:paul.anal...@mail.com>> wrote: new verison and new errors (

Re: [julia-users] Error Pkg.add("LightXML"), where is problem ?on my comp or on server?

2014-06-05 Thread Paul Analyst
b.jl, in expression starting on line 4 while loading C:\Users\SAMSUNG2\.julia\LightXML\src\LightXML.jl, in expression starting on line 24 julia> Paul W dniu 2014-06-05 15:45, Jameson Nash pisze: Version 0.3.0-prerelease+2599 (2014-04-11 23:52 UTC) Commit bf7096c (54 days old master) y

[julia-users] Error Pkg.add("LightXML"), where is problem ?on my comp or on server?

2014-06-05 Thread paul analyst
A fresh approach to technical computing Documentation: http://docs.julialang.org Type "help()" to list help topics Version 0.3.0-prerelease+2599 (2014-04-11 23:52 UTC) Commit bf7096c (54 days old master) x86_64-w64-mingw32 julia> Pkg.add("LightXML") INFO: Cloning cache of LightXML from git://gith

Re: [julia-users] How to read and change the content of web pages to the vector ?

2014-06-05 Thread Paul Analyst
table to STDOUT? -- John On Jun 4, 2014, at 9:33 AM, Kevin Squire <mailto:kevin.squ...@gmail.com>> wrote: I would think "web page" = "HTML document". On Wednesday, June 4, 2014, John Myles White mailto:johnmyleswh...@gmail.com>> wrote: I don't rea

Re: [julia-users] Re: Dynamic creation and naming of variables, what wrong?

2014-06-04 Thread paul analyst
Sorry, is ok , i have vectors .. Paul W dniu środa, 4 czerwca 2014 20:09:46 UTC+2 użytkownik paul analyst napisał: > > Big txh, it works, > but now i heve 100 arrays (in real more... ) with 99 empty columns. Is > posible read m1 , m2 to vector , not to array ? > Paul >

Re: [julia-users] Re: Dynamic creation and naming of variables, what wrong?

2014-06-04 Thread paul analyst
` is being >>> interpreted as the name of its argument, but this is not a valid identifier. >>> >>> Finally, although you almost certainly don't want to do what you're >>> trying to do, you can if you separate creation of the symbol from >>&g

[julia-users] Dynamic creation and naming of variables, what wrong?

2014-06-04 Thread paul analyst
dynamic creation and naming of variables, what wrong? I need 100 new variables : m1 to m 100 julia> for i=1:100 "m" * string(i)=readcsv(string("m",i,".txt")) end ERROR: syntax: "#" is not a valid function argument name

Re: [julia-users] How to read and change the content of web pages to the vector ?

2014-06-04 Thread Paul Analyst
nice will be in 2 varaints: 1 All inside html 2. Only text displayed by browser Paul W dniu 2014-06-04 18:23, John Myles White pisze: I don't really understand what you mean by web page. -- John On Jun 4, 2014, at 9:20 AM, Paul Analyst wrote: Any help ? How to read and chang

Re: [julia-users] How to read and change the content of web pages to the vector ?

2014-06-04 Thread Paul Analyst
Any help ? How to read and change the content of web pages to the vector ["word1", "word2", "word3", ",,,", "wordlast"]? W dniu 2014-06-04 12:47, paul analyst pisze: How to read and change the content of web pages to the vector ["word1", "word2", "word3", ",,,", "wordlast"]? Paul

Re: [julia-users] How to calculate all Euclidean distances between two sets of points in a 3-dimensional space

2014-06-04 Thread paul analyst
It is there. If D aand D1 ar two arrays with vectors in columns You have this in Distance package. julia> r = pairwise(Euclidean(), D,D1) 4x4 Array{Float64,2}: 4579.3 4528.87 2650.71 7326.62 5658.24 3646.92 2587.83 7348.04 5729.42 4581.69 1941.65 7266.23 5733.62 4656.58 2488.6

[julia-users] How to read and change the content of web pages to the vector ?

2014-06-04 Thread paul analyst
How to read and change the content of web pages to the vector ["word1", " word2", "word3", ",,,", "wordlast"]? Paul

Re: [julia-users] pmap(eig,M) . How to separate result of pmap ?

2014-06-01 Thread paul analyst
2014-06-01 16:56 GMT+02:00 paul analyst > > : > >> julia> M = {cor(rand(4,4)) for i=1:2} >> 2-element Array{Any,1}: >> 4x4 Array{Float64,2}: >> 1.0 -0.227554 0.507924 -0.800516 >> -0.227554 1.00.220113 0.75924 >> 0.507924

[julia-users] pmap(eig,M) . How to separate result of pmap ?

2014-06-01 Thread paul analyst
julia> M = {cor(rand(4,4)) for i=1:2} 2-element Array{Any,1}: 4x4 Array{Float64,2}: 1.0 -0.227554 0.507924 -0.800516 -0.227554 1.00.220113 0.75924 0.507924 0.220113 1.0 -0.272765 -0.800516 0.75924 -0.272765 1.0 4x4 Array{Float64,2}: 1.0 -0.65399

Re: [julia-users] Parallel, strange behavior of the loop after changing the value to a variable, 6 time slower

2014-06-01 Thread paul analyst
not defined in anonymous at no file:1379 ERROR: no method +(UndefVarError, UndefVarError) in mr_pairwise at reduce.jl:534 julia> toc() elapsed time: 2.493525081 seconds 2.493525081 julia> W dniu niedziela, 1 czerwca 2014 16:27:14 UTC+2 użytkownik paul analyst napisał: > > > > &

Re: [julia-users] Parallel, strange behavior of the loop after changing the value to a variable, 6 time slower

2014-06-01 Thread paul analyst
const k=2 > > the timing results should be similar. > > > 2014-06-01 16:09 GMT+02:00 paul analyst > > : > >> >> D:\install\Julia\Julia 0.3.0-prerelease-win64-ver3\Julia 0.3.0-prerelease >> ver 3>bin\julia.exe -p 8 >>_

[julia-users] Parallel, strange behavior of the loop after changing the value to a variable, 6 time slower

2014-06-01 Thread paul analyst
D:\install\Julia\Julia 0.3.0-prerelease-win64-ver3\Julia 0.3.0-prerelease ver 3>bin\julia.exe -p 8 _ _ _ _(_)_ | A fresh approach to technical computing (_) | (_) (_)| Documentation: http://docs.julialang.org _ _ _| |_ __ _ | Type "help()" to list

Re: [julia-users] Code running back? what happens?

2014-05-30 Thread Paul Analyst
not the same as a (or b, from the prior stage, for that matter). It doesn't matter whether a and b are numbers or matrices, large or small On Wednesday, May 28, 2014, paul analyst > wrote: Thx, I'm surprised. I suspect, and I did: a = ones (5,5)

[julia-users] Re: Code running back? what happens?

2014-05-28 Thread paul analyst
e identifier FSbis to the same memory that FS is bound to, so > the identifiers are aliases for one another. > > If you want FSbis to start out initialized to the same values as those in > FS, but be a separate container, use `FSbis = copy(FS)`. > > On Wednesday, May 28, 2014 2:07:

[julia-users] Code running back? what happens?

2014-05-28 Thread paul analyst
Code running back? what happens? I have a array FS julia> println(sum(FS)); 9.8267205e7 julia> l,m=size(FS); julia> FSbis=FS; julia> julia> julia> for i=1:l; #println(i) if w[i]==1 us=hcat(F[i,:]',[1:1:m]); us=sortrows(us, by=x->x[1],rev=true); for j=1

Re: [julia-users] Re: How quickly and automatically find the cut-off point of vectors with long tails ?

2014-05-27 Thread paul analyst
t; > will give you the portion of a for which the slope is steeper than 45 > degrees. (I might be off by one here, so you might want to experiment a > little and see if you need to include one more, or one less, data point to > get exactly what you want...) > > // T > > On Tue

Re: [julia-users] Re: How quickly and automatically find the cut-off point of vectors with long tails ?

2014-05-27 Thread paul analyst
an, but > that still depends on the source, distribution, and meaning of your data. > In any case, you'd be better off asking this somewhere like CrossValidated ( > http://stats.stackexchange.com/) -- this isn't a Julia question. > > > On Tue, May 27, 2014 at 8:22 AM, paul ana

[julia-users] Re: How quickly and automatically find the cut-off point of vectors with long tails ?

2014-05-27 Thread paul analyst
Does not work always, distributions are different. How to find the number of elements of the vector from which the chart is getting flat (where is the beginning of the tail?) Paul W dniu wtorek, 27 maja 2014 12:19:09 UTC+2 użytkownik Yuuki Soho napisał: > > The simplest way to do it is probably

[julia-users] Re: How quickly and automatically find the cut-off point of vectors with long tails ?

2014-05-27 Thread paul analyst
W dniu wtorek, 27 maja 2014 12:19:09 UTC+2 użytkownik Yuuki Soho napisał: > > The simplest way to do it is probably to use a quantile: > > > a=[5 3 2 1.5 1.1 1 0.8 0.25 0.2 0.16] > > q = quantile(vec(a),0.1) > > a = a[1:cut] > > > Does not work always, distributions are different. > How to find th

[julia-users] How quickly and automatically find the cut-off point of vectors with long tails ?

2014-05-27 Thread paul analyst
How quickly and automatically cut off long tails We have often sorted vectors with long tails. How quickly and automatically find the cut-off point where the graph flattens? Vectors are long and have a variety of unpleasant characteristics. Simple sample: a=[5 3 2 1.5 1.1 1 0.8 0.25 0.2 0.16]

Re: [julia-users] when the loop variable takes on values from a vector?

2014-05-26 Thread paul analyst
Thx, W dniu poniedziałek, 26 maja 2014 19:23:24 UTC+2 użytkownik Jameson napisał: > > yes, unless i did not understand the question correctly > > On Mon, May 26, 2014 at 1:15 PM, paul analyst > > > wrote: > > if vector a=[1, 4, 56, 67, 99] > > > >

Re: [julia-users] when the loop variable takes on values from a vector?

2014-05-26 Thread paul analyst
if vector a=[1, 4, 56, 67, 99] cod is for i in a F(i) end ? Paul ? W dniu poniedziałek, 26 maja 2014 19:09:03 UTC+2 użytkownik Jameson napisał: > > for i in [1, 4, 56, 67, 99] > F(i) > end > > On Mon, May 26, 2014 at 12:52 PM, paul analyst > > > wrote: &

[julia-users] when the loop variable takes on values from a vector?

2014-05-26 Thread paul analyst
I have some vector : 1 4 56 67 99 and some loop: for i F(i) do somethink end How to assign a variable (i) only value from a vector? Paul

[julia-users] How quickly calculate distances arrays etc.

2014-05-25 Thread paul analyst
How quickly calculate distances arrays etc. This is an array of quotients, on the diagonal must be zero. Are there any features that quickly count this type of arrays? Data in F... k=size(F,2); O=zeros(k,k); for o=0:k-2; for i=1:k-1-o j=i+1+o; O[i,j]=sum(F[:,i])/sum(F[:,j]); O[j,i]=1/O[i,j]; en

Re: [julia-users] Re: how to apply vector on diagonal of matrix ?

2014-05-20 Thread Paul Analyst
16.48.30 UTC+2 skrev paul analyst: b- not is zeros sometimes... c= diagm(a) b=b+c. Thx Paul W dniu wtorek, 20 maja 2014 16:17:11 UTC+2 użytkownik Oliver Lylloff napisał: Hi Paul, if b is just zeros, then b = diagm(a). Best, Oliver

[julia-users] Re: how to apply vector on diagonal of matrix ?

2014-05-20 Thread paul analyst
b- not is zeros sometimes... c= diagm(a) b=b+c. Thx Paul W dniu wtorek, 20 maja 2014 16:17:11 UTC+2 użytkownik Oliver Lylloff napisał: > > Hi Paul, > > if b is just zeros, then b = diagm(a). > > Best, > Oliver > > > > Den tirsdag den 20. maj 2014 15.04.33 UTC+2

[julia-users] how to apply vector on diagonal of matrix ?

2014-05-20 Thread paul analyst
a=rand(5) b=zeros(5,5) how to apply a on diagonal b ? Paul

[julia-users] Algebra, Eigen, Varimax. Precision of calculations.

2014-05-17 Thread paul analyst
Varimax rotation loadings matrix sizes of about 3000x3000, unless as a result of errors of machine completely loses ortogonality. Please help with the rotation of large matrices. Paul

[julia-users] how to Replicate array, translate from matlab

2014-05-16 Thread paul analyst
in matlab : x=[1,4] B=repmat(A,x) /(B==) how to replicate A in Julia Paul

[julia-users] sort of two arrays without concatenate them

2014-05-10 Thread paul analyst
At the same time as a sort of two arrays? k = 1. Matrix A=10 × k matrixB =10 x k I do not want to connect matrix. Im sorting the rows of the matrix B and I would like to sort matrix A by thesame key. Is it possible to link the matrix without concatenate (hcat) them? Paul

[julia-users] Re: How to choose from the array to the new array only selected columns?

2014-04-29 Thread paul analyst
evenly spaced indexes like that, you can > use a range: > > b = a[:, 1:2:5] > > On Tuesday, April 29, 2014 1:23:36 PM UTC-5, paul analyst wrote: >> >> a =rand(5,5) >> b=slicedim(a,2,1,3,5) >> ERROR: no method slicedim(Array{Float64,2}, Int64, Int64, Int64, Int64) >> >> how to choose from matrix a column 1,3,5 wihout while etc.? >> Paul >> >

[julia-users] How to choose from the array to the new array only selected columns?

2014-04-29 Thread paul analyst
a =rand(5,5) b=slicedim(a,2,1,3,5) ERROR: no method slicedim(Array{Float64,2}, Int64, Int64, Int64, Int64) how to choose from matrix a column 1,3,5 wihout while etc.? Paul

Re: [julia-users] how to convert csv file to bin and how to read bin file in Julia

2014-04-29 Thread Paul Analyst
W dniu 2014-04-28 22:34, Tim Holy pisze: Read as CSV, and write as HDF5? --Tim On Monday, April 28, 2014 11:24:06 AM paul analyst wrote: How to convert csv big file (utf8) to bin format ? csv file: 12356.5,12455.5,12545.4,124784.2 12564.2,25488.01,25544.6,125847.2 ... How to read new file

Re: [julia-users] Octave diag and Julia diag

2014-04-28 Thread paul analyst
to make new bug? Paul W dniu poniedziałek, 28 kwietnia 2014 07:26:51 UTC+2 użytkownik Viral Shah napisał: > > I filed https://github.com/JuliaLang/julia/issues/6676 > > -viral > > On Monday, April 28, 2014 4:04:33 AM UTC+5:30, Simon Kornblith wrote: >> >> If diag is passed a vector rather than a

[julia-users] how to convert csv file to bin and how to read bin file in Julia

2014-04-28 Thread paul analyst
How to convert csv big file (utf8) to bin format ? csv file: 12356.5,12455.5,12545.4,124784.2 12564.2,25488.01,25544.6,125847.2 ... How to read new file (bin ) in Julia ? Paul

Re: [julia-users] Re: different eigen value ordering from eigs

2014-04-26 Thread paul analyst
Hi Ethan Changing the subject: Do you have errors in your own serials vectors? With the array size is around 2500 and more last eigenvectors unfortunately often goes out not perpendicular. I am looking for a solution to this problem. Paul W dniu sobota, 26 kwietnia 2014 16:34:52 UTC+2 użytko

Re: [julia-users] how today to create horizontal string vector

2014-04-25 Thread paul analyst
n"].' > 1x3 Array{ASCIIString,2}: > "Andreas" "Noack" "Jensen" > > or use white space > > julia> ["Andreas" "Noack" "Jensen"] > 1x3 Array{ASCIIString,2}: > "Andreas" "Noack&qu

[julia-users] how today to create horizontal string vector

2014-04-25 Thread paul analyst
in old Julia version works this code: stringvector=["string1","string2","string3"]'; today : ERROR: no method conj(ASCIIString) in ctranspose at array.jl:1283 how today to create horizontal string vector ? Paul

[julia-users] How to read full separated column from HDF5 data (array)

2014-04-20 Thread paul analyst
I can read part fo array saved in HDF5 file. Is ok , nice! julia> data = h5read("X.h5", "grupa", (1:3,1)) 3x1 Array{Float64,2}: 0.269948 0.561731 0.894856 But I cant read all first column (:,1) : julia> data = h5read("X.h5", "grupa", (:,1)) ERROR: no method h5read(ASCIIString, ASCIIString, (C

[julia-users] Re: resultant vector, angle between the vectors

2014-04-20 Thread paul analyst
Thanks, I have a lot of these calculations. Check the command-and-loop performance. W dniu sobota, 19 kwietnia 2014 23:33:52 UTC+2 użytkownik Steven G. Johnson napisał: > > > > On Saturday, April 19, 2014 5:32:10 PM UTC-4, Steven G. Johnson wrote: > >> norm(a + b) >> acos(dot(a, b) / (norm(a) *

[julia-users] resultant vector, angle between the vectors

2014-04-19 Thread paul analyst
Is Julia command to do: 1 length of the resultant vector of vectors a, b 2 angle between the vectors a, b​​? Paul

Re: [julia-users] Conditional works like on the selected columns of the matrix

2014-04-19 Thread paul analyst
Big thx, what is "-2" before t ? Paul W dniu sobota, 19 kwietnia 2014 14:47:40 UTC+2 użytkownik Stefan Karpinski napisał: > > scale!(v, 1 .- 2t) > > > On Sat, Apr 19, 2014 at 5:11 AM, paul analyst > > wrote: > >> I have a logical vector t >>

[julia-users] Re: how transform data to matrix x,y

2014-04-19 Thread paul analyst
W dniu sobota, 19 kwietnia 2014 11:24:02 UTC+2 użytkownik paul analyst napisał: > > > > W dniu piątek, 18 kwietnia 2014 21:35:10 UTC+2 użytkownik Douglas Bates > napisał: >> >> On Friday, April 18, 2014 11:43:37 AM UTC-5, paul analyst wrote: >>> >>&g

[julia-users] Re: how transform data to matrix x,y

2014-04-19 Thread paul analyst
W dniu piątek, 18 kwietnia 2014 21:35:10 UTC+2 użytkownik Douglas Bates napisał: > > On Friday, April 18, 2014 11:43:37 AM UTC-5, paul analyst wrote: >> >> How quickly convert data from 3 columns to an array of x, y >> Data are the two columns of text, and the third

[julia-users] Conditional works like on the selected columns of the matrix

2014-04-19 Thread paul analyst
I have a logical vector t julia> t 1x4 BitArray{2}: true false false false and matrix v julia> v 4x4 Array {Float64, 2}: -0.0018187 -0.18174 0.982684 0.0360478 -0.018408 0.98314 0.1819 -0.00297265 -0.995394 -0.0167994 -0.00839779 0.0940118 0.0940681 0.0111097 -0.0342677 0.99491

[julia-users] two vectors , how indexes of second vector assign to first vector ?

2014-04-18 Thread paul analyst
How to assign to vector : e c c d a c d a c d d Indexes of vector: e a c b d I need : e 1 c 3 c 3d 5 a 2 c 3 d 5 a 2 c 3 d 5 d 5 Paul

[julia-users] Re: array with different column types

2014-04-18 Thread paul analyst
Each column is a separate vector? Maybe so? It has quickly afterwards. Paul W dniu czwartek, 17 kwietnia 2014 23:16:00 UTC+2 użytkownik Stéphane Laurent napisał: > > Hello, > > I need to deal with some objects represented as arrays whose some columns > are BigFloat, some columns are Int, some c

[julia-users] how transform data to matrix x,y

2014-04-18 Thread paul analyst
How quickly convert data from 3 columns to an array of x, y Data are the two columns of text, and the third value: x y z a r 0,31 a r 0,00 a s 0,86 a s 0,80 a t 0,28 b r 0,52 b s 0,79 b s 0,86 b t 0,25 b u 0,15 … … 0,29 q r 0,41 q t 0,61 q v 0,62 The data are: millions of rows

Re: [julia-users] BigFloat matrix

2014-04-07 Thread Paul Analyst
Noack Jensen pisze: That is not possible yet. Please copy and paste the commands and errors from the terminal. It makes it easier to figure out what is wrong. 2014-04-07 18:03 GMT+02:00 paul analyst <mailto:paul.anal...@mail.com>>: Friends , case back: R = convert (Matrix

Re: [julia-users] BigFloat matrix

2014-04-07 Thread paul analyst
; of precision > > julia> BigFloat("0.1") > 1.02e-01 > > with 256 bits of precision > > The better option is to initialize with a string. > > > kl. 10:04:54 UTC+1 lørdag 4. januar 2014 skrev Andreas Noack

[julia-users] Re: how to slice vector from 3D array?

2014-03-29 Thread paul analyst
whatever i'm no 1, but it is not important, How to slice vector b not SubArrary... > > a=zeros(10,10,10) > for i=1:10^3 > a[i]=i > end > > OK! I have 3D array. > > Now I need one vecotr from any dimension. > If > b=slice(a,1,1,:) > Julia return SubArray not vector ! > How to slice this (and othe

[julia-users] how to slice vector from 3D array?

2014-03-29 Thread paul analyst
(create 3d array) a=zeros(10,10,10) for i=1:10^3 a[i]=1 end OK! I have 3D array. Now I need one vecotr from any dimension. If b=slice(a,1,1,:) Julia return SubArray not vector ! How to slice this (and other) *vector* from 3D array? Paul

Re: [julia-users] How to use weight in the data.

2014-03-26 Thread Paul Analyst
Thx. It is grate ! This is great! Does the panel does not have all the documentation packages? Paul W dniu 2014-03-25 18:14, Stefan Karpinski pisze: StatsBase provides extensive support for weights: http://statsbasejl.readthedocs.org/en/latest/ On Tue, Mar 25, 2014 at 12:51 PM, paul

[julia-users] How to use weight in the data.

2014-03-25 Thread paul analyst
How to use weight in the data. Often in social research data are used in weight. It is one variable (column) indicates how much a given row in the data set is weighed (average weight = 1). Typical weights are numbers in the range 0.01 - 10.0 but can be any positive numbers. Sums, averages, and

[julia-users] Re: calendar , how to calculate the distance in time ?

2014-03-16 Thread paul analyst
I have data from a base in my calendar like this: 20130628,2 20130701,3 19970224,4 19990230,5 How to calculate the distance in time (in days) between dates and other similar treatments on the calendar? Paweł

<    1   2   3   4   >