[R] finding suitable distribution with default parameter
Dear R-users, Attach with is my data..what i want to do is finding a suitable distribution for my data..I want to run a few test like the poisson and the exponential distribution. Please help me on how to find the p-value for poisson as well as the exponential distribution without knowing the parameter (default setting). Is it possible?? Thanks in advance. love, Anisah length.breakp.start1 length.breakp.start2 length.breakp.start3 length.breakp.start4 length.breakp.start5 17358865 17966995 21306539 27880531 34166504 length.breakp.start6 length.breakp.start7 length.breakp.start8 length.breakp.start9 length.breakp.start10 36111044 36266288 36854306 43786190 44322336 length.breakp.start11 length.breakp.start12 length.breakp.start13 length.breakp.start14 length.breakp.start15 45529444 46302360 53479132 58567262 60564442 length.breakp.start16 length.breakp.start17 length.breakp.start18 length.breakp.start19 length.breakp.start20 72637088 79875476 93155112 94372260 96643396 length.breakp.start21 length.breakp.start22 length.breakp.start23 length.breakp.start24 length.breakp.start25 131781664 132968364 135945080 141788832 156539568 length.breakp.start26 length.breakp.start27 length.breakp.start28 length.breakp.start29 length.breakp.start30 157817896 162399496 168344072 173146584 176302744 length.breakp.start31 length.breakp.start32 length.breakp.start33 length.breakp.start34 breakp.start1 182878168 183946152 185068720 190791232 15463513 breakp.start2 breakp.start3 breakp.start4 breakp.start5 breakp.start6 23060906 23916705 24334228 25558572 26581623 breakp.start7 breakp.start8 breakp.start9 breakp.start10breakp.start11 29519924 32088916 32771112 34613932 35179096 breakp.start12breakp.start13breakp.start14 breakp.start15breakp.start16 37335390 53673024 84317660 93708872 104501808 breakp.start17breakp.start18breakp.start19 breakp.start20breakp.start21 106810172 110909604 122347404 12684 129217008 breakp.start22breakp.start23breakp.start24 breakp.start25 breakp.start1 148519056 150945112 155771432 181069984 12429060 breakp.start2 breakp.start3 breakp.start4 breakp.start5 breakp.start6 15619549 19992908 21354695 29557382 31963410 breakp.start7 breakp.start8 breakp.start9 breakp.start10breakp.start11 32954726 36395876 52052740 55235022 56624024 breakp.start12breakp.start13breakp.start14 breakp.start15breakp.start16 67608312 76061320 85912864 87104384 103784400 breakp.start17breakp.start18breakp.start19 breakp.start20 breakp.start1 104120600 106089480 122625964 142306944 35110312 breakp.start2 breakp.start3 breakp.start4 breakp.start5 breakp.start6 59132150 70759360 73599196 94416760 101484244 breakp.start7 breakp.start8 breakp.start9 breakp.start10breakp.start11 102740180 107703760 118153092 118787588 119000416 breakp.start12breakp.start13breakp.start14 breakp.start15breakp.start16 121491272 123326116 127732804 134564008 148007280 breakp.start1 breakp.start2 breakp.start3 breakp.start4 breakp.start5 23323764 23752341 30169538 31852312 32894227 breakp.start6 breakp.start7 breakp.start8 breakp.start9breakp.start10
[R] finding the suitable distribution
Dear R-users, Attach with is my data..what i want to do is finding a suitable distribution for my data..I want to run a few test like the poisson and the exponential distribution. Please help me on how to find the p-value for poisson as well as the exponential distribution without knowing the parameter. Is it possible?? Thanks in advance. love, Anisah length.breakp.start1 length.breakp.start2 length.breakp.start3 length.breakp.start4 length.breakp.start5 17358865 17966995 21306539 27880531 34166504 length.breakp.start6 length.breakp.start7 length.breakp.start8 length.breakp.start9 length.breakp.start10 36111044 36266288 36854306 43786190 44322336 length.breakp.start11 length.breakp.start12 length.breakp.start13 length.breakp.start14 length.breakp.start15 45529444 46302360 53479132 58567262 60564442 length.breakp.start16 length.breakp.start17 length.breakp.start18 length.breakp.start19 length.breakp.start20 72637088 79875476 93155112 94372260 96643396 length.breakp.start21 length.breakp.start22 length.breakp.start23 length.breakp.start24 length.breakp.start25 131781664 132968364 135945080 141788832 156539568 length.breakp.start26 length.breakp.start27 length.breakp.start28 length.breakp.start29 length.breakp.start30 157817896 162399496 168344072 173146584 176302744 length.breakp.start31 length.breakp.start32 length.breakp.start33 length.breakp.start34 breakp.start1 182878168 183946152 185068720 190791232 15463513 breakp.start2 breakp.start3 breakp.start4 breakp.start5 breakp.start6 23060906 23916705 24334228 25558572 26581623 breakp.start7 breakp.start8 breakp.start9 breakp.start10breakp.start11 29519924 32088916 32771112 34613932 35179096 breakp.start12breakp.start13breakp.start14 breakp.start15breakp.start16 37335390 53673024 84317660 93708872 104501808 breakp.start17breakp.start18breakp.start19 breakp.start20breakp.start21 106810172 110909604 122347404 12684 129217008 breakp.start22breakp.start23breakp.start24 breakp.start25 breakp.start1 148519056 150945112 155771432 181069984 12429060 breakp.start2 breakp.start3 breakp.start4 breakp.start5 breakp.start6 15619549 19992908 21354695 29557382 31963410 breakp.start7 breakp.start8 breakp.start9 breakp.start10breakp.start11 32954726 36395876 52052740 55235022 56624024 breakp.start12breakp.start13breakp.start14 breakp.start15breakp.start16 67608312 76061320 85912864 87104384 103784400 breakp.start17breakp.start18breakp.start19 breakp.start20 breakp.start1 104120600 106089480 122625964 142306944 35110312 breakp.start2 breakp.start3 breakp.start4 breakp.start5 breakp.start6 59132150 70759360 73599196 94416760 101484244 breakp.start7 breakp.start8 breakp.start9 breakp.start10breakp.start11 102740180 107703760 118153092 118787588 119000416 breakp.start12breakp.start13breakp.start14 breakp.start15breakp.start16 121491272 123326116 127732804 134564008 148007280 breakp.start1 breakp.start2 breakp.start3 breakp.start4 breakp.start5 23323764 23752341 30169538 31852312 32894227 breakp.start6 breakp.start7 breakp.start8 breakp.start9breakp.start10 33475734
[R] p-value
Dear R User, say I have this sample of data ( attach with). What i'm going to do is to test whether this data is uniformly distributed or not by finding the p-value. I've tried using the punif command but it gave me the value of 1 of all the data. Any suggestion on R command to find the p-value??Thanks in advance!! Cheers, Anisah 132968364 135945080 156539568 157817896 162399496 168344072 173146584 176302744 182878168 183946152 185068720 190791232 84317660 93708872 106810172 12684 148519056 150945112 155771432 181069984 87104384 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] how to perform a uniform distribution
Dear R-users, I have a data which is in an interval length of start-end. How am i going to perform a uniform distribution? and how to test if it is the uniform distribution?? Suggestion is needed. Many thanks!! start end 132955984 132980744 135891280 135998880 156158176 156920960 157797616 157838176 162055856 162743136 168221296 168466848 173125232 173167936 176267104 176338384 182826240 182930096 183742528 184149776 184401728 185735712 190671888 190910576 84233696 84401624 93636272 93781472 106763368 106856976 126700168 126899840 148511520 148526592 150626144 151264080 155756800 155786064 181004992 181134976 86872080 87336688 104097168 104144032 106053944 106125016 107679864 107727656 118753400 118821776 118963216 119037616 134486048 134641968 147994448 148020112 108509200 108553224 109017536 109795008 110344280 110376160 123771080 123773128 125765312 125911944 130148128 130193296 138017664 138183600 139008736 139229664 147641520 147741232 120330080 120428864 140628560 140696896 141190272 141289440 75657952 75710672 108716872 108717664 126404928 126540016 105848568 106011056 123527136 123638208 75392032 75398288 77032624 77278608 103147592 103181184 120196464 120263696 76775880 76842744 44009616 44088256 137564368 137927776 148377392 148480704 148838944 150553856 154876816 154971520 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] combine elements of list into a vector
Dear R, Can anybody help me on how to combine all the elements in a list into a single vector?? I've tried using unlist command but it gave me wrong output. Below is the example of the vector output that i want and attach with is my list of data. Any suggestion?? Many Thanks breakp.start 6978022 10249966 10955201 11045352 11814604 13847633 16022770 23625856 29676192 30164446 34558312 44009616 53582632 7428353 20225456 33046200 57970280 70011280 70677328 71362456 71824776 74014744 90487648 101058192 102176888 102409592 110591464 116407200 131185424 132395992 137564368 142716752 148377392 148838944 154876816 [[19]] breakp.start 1 6978022 2 10249966 3 10955201 4 11045352 5 11814604 6 13847633 7 16022770 8 23625856 9 29676192 10 30164446 11 34558312 12 44009616 13 53582632 [[20]] breakp.start 1 7428353 2 20225456 3 33046200 4 57970280 5 70011280 6 70677328 7 71362456 8 71824776 9 74014744 10 90487648 11101058192 12102176888 13102409592 14110591464 15116407200 16131185424 17132395992 18137564368 19142716752 20148377392 21148838944 22154876816 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] extracting columns from a list
Dear all, i have 2 lists of data with each of the list contain 14 columns. How am i going to extract column 12 and 13 from each of the list ?? and can i combine my extracted columns to form a single list. Attach with are my data. Your coorperation is highly appreciated. Many thanks Regards, Anisah [[1]] CS(O) id no.anchor ref loc.start loc.end CS(O).size CS(O)ref.density tested loc.start loc.end breakp.start breakp.end den of anchor 1 CS 2.0 3 mmu19 6465196 6978022 5128266 cfa18 55567952 55782336 6978022699273428 2 CS 3.057 mmu19 6992734 102499663257232 17 cfa18 55792632 57688808 10249966 1027743025 3 CS 4.021 mmu19 10277430 10955201 61 31 cfa18 57721864 58419812 10955201 1100852636 4 CS 5.0 2 mmu19 11008526 11045352 36826 54 cfa18 58462088 58518608 11045352 1131686532 5 CS 6.0 7 mmu19 11316865 11814604 497739 14 cfa21 53902028 53514536 11814604 1193722837 6 CS 7.014 mmu19 11937228 1384763319104057 cfa18 40111560 41204940 13847633 1451509117 7 CS 8.0 3 mmu19 14515091 1602277015076792 cfa1 82195232 83511824 16022770 16199850 5 8 CS 9.025 mmu19 16199850 2362585674260063 cfa1 83903856 90638880 23625856 2375455414 9 CS 10.033 mmu19 23754554 2967619259216386 cfa1 90773872 96913624 29676192 2976391816 10CS 11.0 6 mmu19 29763918 30164446 400528 15 cfa11 29919668 30510776 30164446 30611872 8 11CS 12.019 mmu19 30611872 3455831239464405 cfa26 38767664 41958808 34558312 3473440415 12 CSO 13.173 mmu19 34734404 4400961692752128 cfa28 7343952 15872122 44009616 4408825621 13 CSO 13.271 mmu19 44088256 5358263294943767 cfa28 15923283 24830712 53582632 5365479614 [[2]] CS(O) id no.anchor ref loc.start loc.end CS(O).size CS(O)ref.density tested loc.start loc.end breakp.start breakp.end den of anchor 1CSO 2.1 4 mmuX 7311438 7428353 116915 34 cfaX 41732964 41660008 7428353835693216 2CSO 2.241 mmuX 8356932 20225456 118685243 cfaX 32299338 41540632 20225456 2064158413 3CSO 2.3 8 mmuX 20641584 33046200 124046161 cfaX 91770240 94412912 33046200 3317518815 4CSO 3.173 mmuX 33175188 57970280 247950923 cfaX 94538728 114133200 57970280 64939220 3 5CSO 3.229 mmuX 64939220 7001128050720606 cfaX 119319152 124625688 70011280 7003978432 6 CS 4.021 mmuX 70039784 70677328 637544 33 cfaX 124652504 125280776 70677328 7124127229 7CSO 5.1 2 mmuX 71241272 71362456 121184 17 cfaX 125764816 125872960 71362456 7142557614 8CSO 6.1 4 mmuX 71425576 71824776 399200 10 cfaX 125917392 126261784 71824776 73894168 6 9CSO 6.2 2 mmuX 73894168 74014744 120576 17 cfaX 6363930 1898656 74014744 7470540811 10 CSO 6.317 mmuX 74705408 90487648 157822401 cfaX 32034930 19573208 90487648 9122816010 11 CSO 6.449 mmuX 91228160 10105819298300325 cfaX 51733740 61930260101058192 101967264 9 12CS 7.0 3 mmuX 101967264 102176888 209624 14 cfaX 62797768 63194052102176888 10223004011 13 CSO 8.1 2 mmuX 102230040 102409592 179552 11 cfaX 63307688 63460248102409592 103123352 9 14 CSO 8.216 mmuX 103123352 11059146474681122 cfaX 63896752 71000280110591464 116044144 2 15CS 9.0 2 mmuX 116044144 116407200 3630566 cfaX 71543016 71752808116407200 118512328 1 16CS 10.024 mmuX 118512328 131185424 126730962 cfaX 72291072 79468544131185424 13161668820 17CS 11.0 8 mmuX 131616688 132395992 779304 10 cfaX 79821128 80586408132395992 13246176017 18CS 12.021 mmuX 132461760 13756436851026084 cfaX 80674704 8593
[R] extracting elements by using logical values
dear lists,
My question is quite simple but i'm a new user in R and it's seem tough to
me. How am i going to recall back my values??ok..it could be easy if i
ilustrate my problem with the simple example. say that,
x
[1] 1 2
y
[1] 2 3
and my logical output from my function are (Below is my function):
overlapped(x,y)
[1] TRUE
[1] TRUE
suppose i want my output to be:
[1] 1 2
[1] 2 3
overlapped<-function(x,y){
z<-rbind(x,y)
for(i in 1:nrow(z)){
overlap<-(z[i,1]<=z[,2]&z[i,2]>=z[,1])[-i]
if(any(overlap)) print(TRUE)
else print(FALSE)
}
}
Thanks in advance!!
cheers,
Anisah
-
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
[R] help on loop function
Dear lists,
Anysuggestion on how to write a loop function which has nx2 matrix as an
input, where each of the row represents an interval. the function should return
a mx2 matrix containing the no matched interval of x?? i've tried my own
function but it's only applies for 2x2 matrix.
Here is my code:
overlap<-function(x,y){
m=length(x)
n=length(y)
for(i in 1:m){
for(j in 1:n){
if(x[i]>=y[j]) print(TRUE)
else print(FALSE)
}
}
}
Many thanks!!
Cheers,
Anisah
-
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
[R] writing a simple function
Dear lists,
any suggestion on how to write a function to return a TRUE if interval [a,b]
overlaps the interval [c,d]. I've tried it but an error occur saying that
'could not find function v ; in addition warning message occur'. Below is my
codes.Please help me sort this problem as i'm in process of learning of writing
a function. Many thanks
overlap<-function(m,n){
t=length(m)
v=length(n)
tt=logical(t)
tv=logical(v)
for(i in 1:m){
for(i in 1:n){
if(v(i,j)<=t(i,j)){
tt=T
tv=T
}
}
}
k=c(tt,tv)
return(k)
}
Regards,
Anisah
-
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
[R] reverse vector elements
Dear lists,
I want to write a function of a vector and reverse the order of its elements.
Here is my code:
revector<-function(n){
y=vector(length=n)
for(i in n:1){
y[i]=i
}
return(y)
}
i want my output to be like this:
y
[1] 10 9 8 7 6 5 4 3 2 1
Any suggestion?? Thanks!!
Cheers,
Anisah
-
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
[R] error in the function
Dear lists,
i want to find the non-overlapping interval values with this code:
mysetdiff=function(x,y){
m=length(x)
n=length(y)
bx = logical(m)
by = logical(n)
for(i in 1:m){
for(j in 1:n){
if(x[i]<=y[j+1]){
bx[i] = T
by[j] = T
NA= NA
}
}
}
sx = x[!bx]
sy = y[!by]
s=c(sx,sy)
return(s)
}
Below is my dataset. When i called back my function with the
code;>mysetdiff(f,e). An error had occur: Error in if (x[i] <= y[j + 1]) { :
missing value where TRUE/FALSE needed. How am i going to fix my function so
that i can get the values of my non-overlapping interval. Any suggestion??
Thanks a bunch!!
e
[1] 17130612 17712302 21225764 25012714 33852816 36012944 36252300
[8] 36737468 43693832 44148616 45318876 45852632 53258208 58530988
[15] 60437872 72516480 79673224 93128744 94269896 95868704 99651504
[22] 113688560 131101008 132955984 135891280 141318144 148257888 156158176
[29] 157797616 162055856 168221296 173125232 176267104 182826240 183742528
[36] 184401728 190671888 196639616 17587118 18221688 21387314 30748348
[43] 34480192 36209144 36280276 36971144 43878548 44496056 45740012
[50] 46752088 53700056 58603536 60691012 72757696 80077728 93181480
[57] 94474624 97418088 106596368 120128352 132462320 132980744 135998880
[64] 142259520 151591840 156920960 157838176 162743136 168466848 173167936
[71] 176338384 182930096 184149776 185735712 190910576
> f
[1] 17712302 21203780 25012714 33852816 34794536 36012944 37891284
[8] 43693832 44148616 45852632 53289188 61573112 63664928 72516480
[15] 79673224 94474624 95868704 99651504 113688560 125159688 127388568
[22] 131101008 154599216 176267104 181504912 182562720 182826240 183742528
[29] 196841904 18100404 21387314 30748348 34384588 35996440 36252300
[36] 37942556 43878548 44496056 46752088 53700056 62637560 63969952
[43] 72757696 80077728 94617360 97144032 106596368 120128352 127220456
[50] 127504536 132462320 154717312 176338384 181836032 182687824 182930096
[57] 184149776
-
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
[R] remove the missing value,NA
I have two sets of interval data.Below are my two dataset. In these dataset,
there is a missing values in each of the data. I want to find the
non-overlapping interval values. Here is my code:
mysetdiff=function(x,y){
m=length(x)
n=length(y)
bx = logical(m)
by = logical(n)
for(i in 1:m){
for(j in 1:n){
if(x[i]<=y[j]){
bx[i] = T
by[j] = T
}
}
}
sx = x[!bx]
sy = y[!by]
s=c(sx,sy)
return(s)
}
How am i going to remove the NA value in my function?? any sugesstion??
thanks a bunch
c
[1] 17130612 17712302 21225764 25012714 33852816 36012944 36252300
[8] 36737468 43693832 44148616 45318876 45852632 53258208 58530988
[15] 60437872 72516480 79673224 93128744 94269896 95868704 99651504
[22] 113688560 131101008 132955984 135891280 141318144 148257888 156158176
[29] 157797616 162055856 168221296 173125232 176267104 182826240 183742528
[36] 184401728 190671888 196639616 17587118 18221688 21387314 30748348
[43] 34480192 36209144 36280276 36971144 43878548 44496056 45740012
[50] 46752088 53700056 58603536 60691012 72757696 80077728 93181480
[57] 94474624 97418088 106596368 120128352 132462320 132980744 135998880
[64] 142259520 151591840 156920960 157838176 162743136 168466848 173167936
[71] 176338384 182930096 184149776 185735712 190910576NA
> d
[1] 17712302 21203780 25012714 33852816 34794536 36012944 37891284
[8] 43693832 44148616 45852632 53289188 61573112 63664928 72516480
[15] 79673224 94474624 95868704 99651504 113688560 125159688 127388568
[22] 131101008 154599216 176267104 181504912 182562720 182826240 183742528
[29] 196841904 18100404 21387314 30748348 34384588 35996440 36252300
[36] 37942556 43878548 44496056 46752088 53700056 62637560 63969952
[43] 72757696 80077728 94617360 97144032 106596368 120128352 127220456
[50] 127504536 132462320 154717312 176338384 181836032 182687824 182930096
[57] 184149776NA
-
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
[R] writing a function
Dear lists,
I'm in my process of learning of writing a function. I tried to write a
simple functions of a matrix and a vector. Here are the codes:
mm<-function(m,n){ #matrix function
w<-matrix(nrow=m, ncol=n)
for(i in 1:m){
for(j in 1:n){
w[i,j]=i+j
}
}
return(w[i,j])
}
v<-function(n){ #function of a vector
y=vector(length=n)
for(i in 1:n){
y[i]=i
}
return(y[i])
}
when i tried called on my matrix function; say mm(5,10). what i got is the
single value of 15..where's is my matrix?? same thing happened to my vector
function. Please help me figure out these problems. Thanks a bunch!!
Cheers,
Anisah
-
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
[R] matrix loop
Dear list,
I'm trying to make a loop of a (5x10) matrix and below are my codes. Could
anybody help me figure out why my loop is not working. Thanks in advance!!
m<-1:5
n<-1:10
for(i in 1:length(m))
{ for(j in 1:length(n))
{
y[i,j]=sum(i,j)
y<-as.matrix(y[i,j])
}
}
cheers,
Anisah
-
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
[R] Vector loop
hi, I'm in my learning process of doing a programming with "for" loop. How to make a loop of a vector of length 10 where elements are 1,2,3,4,5,6,7,8,9,10. Any suggestion needed!! Many thanks. Cheers, Anisah - [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] loops
hi!! Below I have 4 columns vector of c and d which are unequal in length.These c and d have 2 columns each where these 2 columns represent an interval values. I easily get these overlapping interval values by using these code: s<-apply(c,1,rev) aa<-apply(d,1,function(x)apply(s<=x,2,any)) My huge problem is i need to do this detection using the loops and i'm seriously not good in doing that. I've tried to read the R manual but still failed to do so. Can somebody help me as i'm a new user to R!! Thanks in advance.. c d 17130612 17587118 17712302 18100404 17712302 18221688 21203780 21387314 21225764 21387314 25012714 30748348 25012714 30748348 33852816 34384588 33852816 34480192 34794536 35996440 36012944 36209144 36252300 36280276 36737468 36971144 43693832 43878548 - [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] overlapping intervals
hi!! Below I have 4 columns vector of c and d which are unequal in length.These c and d have 2 columns each where these 2 columns represent an interval values. How am I going to get an overlapping over these interval values?? Please help me sort this problem!! Thanks in advance.. c d 17130612 17587118 17712302 18100404 17712302 18221688 21203780 21387314 21225764 21387314 25012714 30748348 25012714 30748348 33852816 34384588 33852816 34480192 34794536 35996440 36012944 36209144 36252300 36280276 36737468 36971144 43693832 43878548 - [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] remove similar values
hello!! say that i have the values of x1 and x2. my x1 has longer length than x2. how am i going to remove the similar values of x1 and x2. Any suggestion?? Thanks in advance!! Cheers, Anisah e f [1,] 17358865 17906353 [2,] 17966995 21295547 [3,] 21306539 27880531 [4,] 27880531 34118702 [5,] 34166504 35395488 [6,] 36111044 36132622 [7,] 36266288 37916920 [8,] 36854306 43786190 [9,] 43786190 44322336 [10,] 44322336 46302360 [11,] 45529444 53494622 [12,] 46302360 62105336 [13,] 53479132 63817440 [14,] 58567262 72637088 [15,] 60564442 79875476 [16,] 72637088 94545992 [17,] 79875476 96506368 [18,] 93155112 103123936 [19,] 94372260 116908456 [20,] 96643396 126190072 [21,] 103123936 127446552 [22,] 116908456 131781664 [23,] 131781664 154658264 [24,] 132968364 176302744 [25,] 135945080 181670472 [26,] 141788832 182625272 [27,] 149924864 182878168 [28,] 156539568 183946152 [29,] 157817896NA [30,] 162399496 [31,] 168344072 [32,] 173146584 [33,] 176302744 [34,] 182878168 [35,] 183946152 [36,] 185068720 [37,] 190791232 [38,]NA - [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] x-axis
Hi,
I want to plot a graph and here is my code:
ec<-rep(0,length(e))
fc<-rep(0,length(f))
plot(e,ec,type="p",col=1,pch=19)
points(f,fc,col=2,pch=20)
legend(1.0e+08,1.0,c("dog", "human"),text.col="green4",pch=c(19,20),col=c(1,2))
my major problem here is the x-axis is too large in scale and its very hard
for me to read the x-points from my graph. Can I change the x-label too??
Please help me figure out my problem!! Thanks in advance
e:
[1] 17358865 17966995 21306539 27880531 34166504 36111044 36266288
[8] 36854306 43786190 44322336 45529444 46302360 53479132 58567262
[15] 60564442 72637088 79875476 93155112 94372260 96643396 103123936
[22] 116908456 131781664 132968364 135945080 141788832 149924864 156539568
[29] 157817896 162399496 168344072 173146584 176302744 182878168 183946152
[36] 185068720 190791232NA
f<-[1] 17906353 21295547 27880531 34118702 35395488 36132622 37916920
[8] 43786190 44322336 46302360 53494622 62105336 63817440 72637088
[15] 79875476 94545992 96506368 103123936 116908456 126190072 127446552
[22] 131781664 154658264 176302744 181670472 182625272 182878168 183946152
[29]NA
cheers,
Anisah
-
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
[R] adjustment on the x-axis scale width
Hi,
I have sucessfully plotted the graph by using the following command:
ec<-rep(0,length(e))
fc<-rep(0,length(f))
plot(e,ec,type="p",col=1,pch=19)
points(f,fc,col=2,pch=20)
legend(1.0e+08,1.0,c("dog", "human"),text.col="green4",pch=c(19,20),col=c(1,2))
Below are the values for e and f...the questions is how am i going to adjust
the scale width of my x-axis and my y-axis will start at 0?? Any suggestion??
thanks in advance!!
e:
[1] 17358865 17966995 21306539 27880531 34166504 36111044 36266288
[8] 36854306 43786190 44322336 45529444 46302360 53479132 58567262
[15] 60564442 72637088 79875476 93155112 94372260 96643396 103123936
[22] 116908456 131781664 132968364 135945080 141788832 149924864 156539568
[29] 157817896 162399496 168344072 173146584 176302744 182878168 183946152
[36] 185068720 190791232NA
f<-[1] 17906353 21295547 27880531 34118702 35395488 36132622 37916920
[8] 43786190 44322336 46302360 53494622 62105336 63817440 72637088
[15] 79875476 94545992 96506368 103123936 116908456 126190072 127446552
[22] 131781664 154658264 176302744 181670472 182625272 182878168 183946152
[29]NA
cheers,
Anisah
-
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
[R] stripchart
hi, I want to plot a dotplot graph but unfortunately R does not have it any more. It suggest to use either the stripchart or the dotchart and i prefer to play around with the stripchart graph. Sadly, the result is not like i want where I'm actually wants the plots/points were scattered around the x-axis. Plus, from the results shows that and the x-axis gave a very large scale, i.e 5.0e+0.7. How am i going to change this scale and to change the colour or symbols where i can show the difference of the sample data. Any suggestion?? two sample data: e [1] 17358865 17966995 21306539 27880531 34166504 36111044 36266288 [8] 36854306 43786190 44322336 45529444 46302360 53479132 58567262 [15] 60564442 72637088 79875476 93155112 94372260 96643396 103123936 [22] 116908456 131781664 132968364 135945080 141788832 149924864 156539568 [29] 157817896 162399496 168344072 173146584 176302744 182878168 183946152 [36] 185068720 190791232NA > f [1] 17906353 21295547 27880531 34118702 35395488 36132622 37916920 [8] 43786190 44322336 46302360 53494622 62105336 63817440 72637088 [15] 79875476 94545992 96506368 103123936 116908456 126190072 127446552 [22] 131781664 154658264 176302744 181670472 182625272 182878168 183946152 [29]NA R-command: stripchart(e,method="stack",col="green",add=F,at=0.0) stripchart(f,method="stack",col="blue",add=T,at=0.0) Thanks in advance. Cheers, Anisah - [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] spacing problems
Dear all,
I try to do a matrix with this data but the warning massage had occur:
In matrix(x5, ncol = 6, byrow = T) :
data length [116] is not a sub-multiple or multiple of the number of rows [20]
I think the spacing problem cause this where double spacing and tabs are
everywhere!!
Secondly, I have a raw data (which i downloaded from the AutoGraph server)
where it has 2 data in one text file. How am i going to split it to two?? I
also need to make a command function so that i can do it for the other 38
similar types of data. Below show my command function and attach with is my
data. Kindly help me to solve these problems and thanks in advance.
fun<-function(filename)
{
x<-scan(file=filename,sep="\n",skip=12,what=character(0))
x1<-gsub("[][)(:\\,|\\-]","",x)
x2<-gsub("Telomere","NA",x1)
x3<-gsub("Decreasing order|Increasing order","",x2)
x4<-strsplit(x3,"\t")
x5<-unlist(x4)
y<-matrix(x5,ncol=6,byrow=T)
tc<-textConnection(apply(y,1,paste,collapse=" "))
w<-read.table(tc)
t<-as.data.frame(w)
attr(t,"names")=c("CS(O)","id","no.anchor","ref","loc.start","loc.end","CS(O).size","CS(O)ref.density","tested","loc.start","loc.end","breakp.start","breakp.end","den
of anchor")
return(t)
}
Cheers,
Anisah
-
AutoGRAPH analysis in a flat-file format ( CS(O) locations, size,
breakpoints...):
--
Reference Chromosome vs Dataset 3
-
CS(O) id (number of marker/anchor) Location(s) on reference CS(O)
size CS(O) density on reference chromosome Location(s) on tested
Breakpoints CS(O) locations (denstiy of marker/anchor)
CS 1 (27): cfa1: [ 3251712 - 12398289 ] 91465773 mmu18:
[ 24330828 - 90644456 ] ] 12398289, 13347136 [(2 )
CS 2 (19): cfa1: [ 13347136 - 18193820 ]48466844 mmu1:
[ 113688560 - 106596368 ]] 18193820, 19140840 [(3 )
CSO 3.1 (18):cfa1: [ 19140840 - 22178912 ] 3038072 6
mmu18: [ 66984412 - 63788168 ]- Decreasing order -] 22178912, 23188292 [
(2 )
CSO 3.2 (4): cfa1: [ 23188292 - 24126920 ] 938628 4
mmu18: [ 69469888 - 70578512 ]- Increasing order -] 24126920, 24190392 [(8
)
CS 4 (2):cfa1: [ 24190392 - 24265894 ]75502 26 mmu18:
[ 70634048 - 70693560 ] ] 24265894, 24823786 [(7 )
CSO 5.1 (6): cfa1: [ 24823786 - 27113036 ] 2289250 3
mmu18: [ 71384360 - 73984760 ]- Increasing order -] 27113036, 27418228 [
(13 )
CSO 5.2 (4): cfa1: [ 27418228 - 27578150 ] 159922 25
mmu18: [ 68532272 - 68058624 ]- Decreasing order -] 27578150, 28055666 [(9
)
CS 6 (77): cfa1: [ 28055666 - 47327576 ]19271910 4 mmu10:
[ 24284924 - 3134304 ] ] 47327576, 47940412 [(5 )
CSO 7.1 (15):cfa1: [ 47940412 - 51570228 ] 3629816 4
mmu17: [ 3283055 - 7577392 ]- Increasing order - ] 51570228, 51988448 [
(11 )
CSO 7.2 (16):cfa1: [ 51988448 - 57900888 ] 5912440 3
mmu17: [ 12850771 - 8003706 ]- Decreasing order - ] 57900888, 58482632 [
(11 )
CSO 7.3 (5): cfa1: [ 58482632 - 59714992 ] 1232360 4
mmu17: [ 13497139 - 14565909 ]- Increasing order -] 59714992, 59864308 [(15
)
CSO 8.1 (8): cfa1: [ 59864308 - 60129744 ] 265436 30
mmu10: [ 33840808 - 33608556 ]- Decreasing order -] 60129744, 60211840 [
(17 )
CSO 8.2 (21):cfa1: [ 60211840 - 65147456 ] 4935616 4
mmu10: [ 51273836 - 57483016 ]- Increasing order -] 65147456, 65273344 [
(7 )
CSO 8.3 (21):cfa1: [ 65273344 - 72417520 ] 7144176 3
mmu10: [ 33201750 - 24858506 ]- Decreasing order -] 72417520, 73256040 [(7
)
CS 9 (24): cfa1: [ 73256040 - 78879792 ]56237524 mmu13:
[ 64301608 - 58167304 ] ] 78879792, 79088752 [(10 )
CS 10 (3): cfa1: [ 79088752 - 80616144 ]15273922 mmu4:
[ 73484664 - 71603504 ] ] 80616144, 82195232 [(2 )
CS 11 (58): cfa1: [ 82195232 - 96913624 ]14718392 4 mmu19:
[ 14515091 - 29676192 ] ] 96913624, 97163832 [(8 )
CS 12 (10): cfa1: [ 97163832 - 99948816 ]27849844 mmu13:
[ 53424904 - 51664268 ] ] 99948816, 100087840 [(9 )
CSO 13.1 (4):cfa1: [ 100087840 - 100423328 ] 335488 12
mmu13: [ 51443540 - 51113392 ]- Decreasing order -] 100423328, 101013264
[ (11 )
CSO 13.2 (17): cfa1: [ 101013264 - 102120080 ] 1106816 15
mmu13: [ 48589808 - 49694080 ]- Increasing order -] 102271920, 102458192
[(25 )
CSO 14.1 (10): cfa1: [ 102458192 - 104863664 ] 2405472 4
mmu7: [ 11943329 - 5734576 ]- Decreasing order - ] 104863664, 105135376
[ (35 )
[R] splitting the column
Hi, I have a matrix data with 21 rows and 6 columns. Below and attach with is my matrix data. My problem is to split the column into 3 subcolumns (except for column 3 and 4) . I want my new matrix data to have 14 columns and 21 rows. Kindly help is highly appreciated. [,1] [,2] [,3] [,4] [,5] [,6] [1,] "CS 1 73 " " cfa1 3251712 24126920 ""20875208" "3" " hsa18 132170848 50139168 " " 24126920 24153560 8 " [2,] "CS 2 3 " " cfa1 24153560 24265894 " "112334" "27" " hsa18 50105060 49934572 " " 24265894 24823786 7 " [3,] "CSO 3.1 6 " " cfa1 24823786 27113036 " "2289250" "3" " hsa18 48121156 46579500 " " 27113036 27418228 13" [4,] "CSO 3.2 4 " " cfa1 27418228 27578150 " "159922" "25" " hsa18 13872043 13208795 " " 27578150 28055666 9 " [5,] "CS 4 4 " " cfa1 28055666 28835230 " "779564" "5" " hsa6 132311008 133132200 " " 28835230 29482792 7 " [6,] "CS 5 46 " " cfa1 29482792 40120672 " "10637880" "4" " hsa6 133604208 146227152 " " 40120672 40539680 8 " [7,] "CS 6 9 " " cfa1 40539680 43339444 " "2799764" "3" " hsa6 146390608 149867328 " " 43339444 43390788 13 " [8,] "CSO 7.1 74 " " cfa1 43390788 59714992 " "16324204" "5" " hsa6 149929104 169714432 " " 59714992 59864308 15" [9,] "CSO 7.2 52 " " cfa1 59864308 72417520 " "12553212" "4" " hsa6 116707976 131508152 " " 72417520 73256040 7 " [10,] "CSO 8.1 12 " " cfa1 73256040 75192808 " "1936768" "6" " hsa9 98441680 96360824 " " 75192808 75272528 6" [11,] "CSO 8.2 56 " " cfa1 75272528 91881664 " "16609136" "3" " hsa9 89530256 70341312 " " 91881664 92281272 5" [12,] "CSO 8.3 22 " " cfa1 92281272 96913624 " "4632352" "5" " hsa9 261625 5755076 " " 96913624 98067040 5" [13,] "CSO 8.4 15 " " cfa1 98067040 100692560 " "2625520" "6" " hsa9 93833248 89771184 " " 100692560 101013264 13" [14,] "CSO 8.5 18 " " cfa1 101013264 102120080 " "1106816" "16" " hsa9 95832896 94012312 " " 102271920 102458192 25 " [15,] "CS 9 55 " " cfa1 102458192 105936824 " "3478632" "16" " hsa19 63765096 59618416 " " 105936824 106097392 35 " [16,] "CSO 10.1 81 " " cfa1 106097392 110263696 " "4166304" "19" " hsa19 59386008 54256216 " " 110263696 110288752 60" [17,] "CSO 10.2 18 " " cfa1 110288752 110567608 " "278856" "65" " hsa19 54163196 53814360 " " 110567608 110575576 50" [18,] "CSO 10.3 60 " " cfa1 110575576 112727048 " "2151472" "28" " hsa19 53649284 50959884 " " 112727048 112775144 40 " [19,] "CS 11 173 " " cfa1 112775144 119848336 " "7073192" "24" " hsa19 50887772 40849556 " " 119848336 119880560 55 " [20,] "CS 12 33 "" cfa1 119880560 121690672 " "1810112" "18" " hsa19 40824500 38556448 " " 121690672 121820640 16 " [21,] "CS 13 22 "" cfa1 121820640 124798800 " "2978160" "7" " hsa19 38391408 34709332 " " 124798800 NA NA " Thanks in advance. Cheers, Anisah - [,1] [,2] [,3] [,4] [,5] [,6] [1,] "CS 1 73 " " cfa1 3251712 24126920 ""20875208" "3" " hsa18 132170848 50139168 " " 24126920 24153560 8 " [2,] "CS 2 3 " " cfa1 24153560 24265894 " "112334" "27" " hsa18 50105060 49934572 " " 24265894 24823786 7 " [3,] "CSO 3.1 6 " " cfa1 24823786 27113036 " "2289250" "3" " hsa18 48121156 46579500 " " 27113036 27418228 13" [4,] "CSO 3.2 4 " " cfa1 27418228 27578150 " "159922" "25" " hsa18 13872043 13208795 " " 27578150 28055666 9 " [5,] "CS 4 4 " " cfa1 28055666 28835230 " "779564" "5" " hsa6 132311008 133132200 " " 28835230 29482792 7 " [6,] "CS 5 46 " " cfa1 29482792 40120672 " "10637880" "4" " hsa6 133604208 146227152 " " 40120672 40539680 8 " [7,] "CS 6 9 " " cfa1 40539680 43339444 " "2799764" "3" " hsa6 146390608 149867328 " " 43339444 43390788 13 " [8,] "CSO 7.1 74 " " cfa1 43390788 59714992 " "16324204" "5" " hsa6 149929104 169714432 " " 59714992 59864308 15" [9,] "CSO 7.2 52 " " cfa1 59864308 72417520 " "12553212" "4" " hsa6 116707976 131508152 " " 72417520 73256040 7 " [10,] "CSO 8.1 12 " " cfa1 73256040 75192808 " "1936768" "6" " hsa9 98441680 96360824 " " 75192808 75272528 6" [11,] "CSO 8.2 56 " " cfa1 75272528 91881664 " "16609136" "3" " hsa9 89530256 70341312 " " 91881664 92281272 5" [12,] "CSO 8.3 22 " " cfa1 92281272 96913624 " "4632352" "5" " hsa9 261625 5755076 " " 96913624 98067040 5" [13,] "CSO 8.4 15 " " cfa1 98067040 100692560 " "2625520" "6" " hsa9 93833248 89771184 " " 100692560 101013264 13" [14,] "CSO 8.5 18 " " cfa1 101013264 102120080 " "1106816" "16" " hsa9 95832896 94012312 " " 102271920 102458192 25 " [15,] "CS 9 55 " " cfa1 102
[R] numerical data frame
Dear All,
I've successfully import my synteny data to R by using scan command. Below
show my results. My major problem with my data is how am i going to combine the
column names with the data( splt) where i have tried on cbind but a warning
message occur. I have realized that the splt data only have 5 column instead of
6. Please help me with this!!
I want my data to be a numerical data with a proper column and column names
and to replace CS with 1 and CSO with 0 and also to get remove all the
punctuations and the characters from the data.
Attach herewith is my original data. Your kindly help is highly appreciated
and thanks in advance.
Cheers,
Anisah
1)for col names
nms<-scan("C:/Users/user/Documents/cfa-1.txt",sep="\t",nlines=1,skip=10,what=character(0))
Read 6 items
> nms
[1] "CS(O) id (number of marker/anchor) "
[2] " Location(s) on reference "
[3] "CS(O) size"
[4] "CS(O) density on reference chromosome"
[5] "Location(s) on tested "
[6] "Breakpoints CS(O) locations (denstiy of marker/anchor)"
2) my data
x<-scan("C:/Users/user/Documents/cfa-1.txt",sep="\n",skip=12,what=character(0))
Read 21 items
> splt<-strsplit(x,"\t")
> splt
[[1]]
[1] "CS 1 (73): " " cfa1: [ 3251712 - 24126920 ] "
[3] " 20875208 " " 3 "
[5] " hsa18: [ 132170848 - 50139168 ] " "] 24126920, 24153560 [(8 ) "
[[2]]
[1] "CS 2 (3): " " cfa1: [ 24153560 - 24265894 ] "
[3] " 112334 "" 27 "
[5] " hsa18: [ 50105060 - 49934572 ] " "] 24265894, 24823786 [(7 ) "
[[3]]
[1] "CSO 3.1 (6): "
[2] " cfa1: [ 24823786 - 27113036 ] "
[3] " 2289250 "
[4] " 3 "
[5] " hsa18: [ 48121156 - 46579500 ]- Decreasing order - ] 27113036, 27418228 [
(13)"
[[4]]
[1] "CSO 3.2 (4): "
[2] " cfa1: [ 27418228 - 27578150 ] "
[3] " 159922 "
[4] " 25 "
[5] " hsa18: [ 13872043 - 13208795 ]- Decreasing order - ] 27578150, 28055666
[(9 ) "
[[5]]
[1] "CS 4 (4): "" cfa1: [ 28055666 - 28835230 ] "
[3] " 779564 " " 5 "
[5] " hsa6: [ 132311008 - 133132200 ] " "] 28835230, 29482792 [(7 ) "
[[6]]
[1] "CS 5 (46): " " cfa1: [ 29482792 - 40120672 ] "
[3] " 10637880 " " 4 "
[5] " hsa6: [ 133604208 - 146227152 ] " "] 40120672, 40539680 [(8 ) "
[[7]]
[1] "CS 6 (9): "" cfa1: [ 40539680 - 43339444 ] "
[3] " 2799764 "" 3 "
[5] " hsa6: [ 146390608 - 149867328 ] " "] 43339444, 43390788 [(13 ) "
[[8]]
[1] "CSO 7.1 (74): "
[2] " cfa1: [ 43390788 - 59714992 ] "
[3] " 16324204 "
[4] " 5 "
[5] " hsa6: [ 149929104 - 169714432 ]- Increasing order -] 59714992, 59864308 [
(15)"
[[9]]
[1] "CSO 7.2 (52): "
[2] " cfa1: [ 59864308 - 72417520 ] "
[3] " 12553212 "
[4] " 4 "
[5] " hsa6: [ 116707976 - 131508152 ]- Increasing order - "
[6] "] 72417520, 73256040 [(7 ) "
[[10]]
[1] "CSO 8.1 (12): "
[2] " cfa1: [ 73256040 - 75192808 ] "
[3] " 1936768 "
[4] " 6 "
[5] " hsa9: [ 98441680 - 96360824 ]- Decreasing order - "
[6] "] 75192808, 75272528 [ "
[7] " (6 )"
[[11]]
[1] "CSO 8.2 (56): "
[2] " cfa1: [ 75272528 - 91881664 ] "
[3] " 16609136 "
[4] " 3 "
[5] " hsa9: [ 89530256 - 70341312 ]- Decreasing order - "
[6] "] 91881664, 92281272 [ "
[R] How to import ENSEMBL text data using R
Dear all, I have a data which is in text file and i would like to import the data to R. From the manual, ive found the read.table command function is the most appropriate but when i wrote the command an error had occur. It say Error in read.table"C:/Users/user/Documents/cfa-1.txt", header = T, sep = "\t",skip=10) :more columns than column names. Please help me with this as im a first time user to R. Thanks in advance. Cheers, Anisah - Listing breakpoints, Conserved Segments and Conserved Segments Ordered locations : -- Reference Chromosome (1) vs Dataset 2 - CS(O) id (number of marker/anchor) Location(s) on reference CS(O) size CS(O) density on reference chromosome Location(s) on tested Breakpoints CS(O) locations (denstiy of marker/anchor) CS 1 (73): cfa1: [ 3251712 - 24126920 ] 20875208 3 hsa18: [ 132170848 - 50139168 ]] 24126920, 24153560 [(8 ) CS 2 (3):cfa1: [ 24153560 - 24265894 ]112334 27 hsa18: [ 50105060 - 49934572 ] ] 24265894, 24823786 [(7 ) CSO 3.1 (6): cfa1: [ 24823786 - 27113036 ]22892503 hsa18: [ 48121156 - 46579500 ]- Decreasing order - ] 27113036, 27418228 [ (13) CSO 3.2 (4): cfa1: [ 27418228 - 27578150 ]159922 25 hsa18: [ 13872043 - 13208795 ]- Decreasing order - ] 27578150, 28055666 [(9 ) CS 4 (4):cfa1: [ 28055666 - 28835230 ]779564 5 hsa6: [ 132311008 - 133132200 ]] 28835230, 29482792 [(7 ) CS 5 (46): cfa1: [ 29482792 - 40120672 ]10637880 4 hsa6: [ 133604208 - 146227152 ]] 40120672, 40539680 [(8 ) CS 6 (9):cfa1: [ 40539680 - 43339444 ]27997643 hsa6: [ 146390608 - 149867328 ]] 43339444, 43390788 [(13 ) CSO 7.1 (74):cfa1: [ 43390788 - 59714992 ]16324204 5 hsa6: [ 149929104 - 169714432 ]- Increasing order -] 59714992, 59864308 [ (15) CSO 7.2 (52):cfa1: [ 59864308 - 72417520 ]12553212 4 hsa6: [ 116707976 - 131508152 ]- Increasing order -] 72417520, 73256040 [(7 ) CSO 8.1 (12):cfa1: [ 73256040 - 75192808 ]19367686 hsa9: [ 98441680 - 96360824 ]- Decreasing order - ] 75192808, 75272528 [ (6 ) CSO 8.2 (56):cfa1: [ 75272528 - 91881664 ]16609136 3 hsa9: [ 89530256 - 70341312 ]- Decreasing order - ] 91881664, 92281272 [ (5 ) CSO 8.3 (22):cfa1: [ 92281272 - 96913624 ]46323525 hsa9: [ 261625 - 5755076 ]- Increasing order - ] 96913624, 98067040 [ (5 ) CSO 8.4 (15):cfa1: [ 98067040 - 100692560 ] 26255206 hsa9: [ 93833248 - 89771184 ]- Decreasing order - ] 100692560, 101013264 [ (13 ) CSO 8.5 (18):cfa1: [ 101013264 - 102120080 ] 110681616 hsa9: [ 95832896 - 94012312 ]- Decreasing order -] 102271920,102458192 [(25 ) CS 9 (55): cfa1: [ 102458192 - 105936824 ] 347863216 hsa19: [ 63765096 - 59618416 ] ] 105936824, 106097392 [(35 ) CSO 10.1 (81): cfa1: [ 106097392 - 110263696 ] 416630419 hsa19: [ 59386008 - 54256216 ]- Decreasing order - ] 110263696,110288752 [ (60 ) CSO 10.2 (18): cfa1: [ 110288752 - 110567608 ] 278856 65 hsa19: [ 54163196 - 53814360 ]- Decreasing order - ] 110567608,110575576 [ (50 ) CSO 10.3 (60): cfa1: [ 110575576 - 112727048 ] 215147228 hsa19: [ 53649284 - 50959884 ]- Decreasing order - ] 112727048,112775144 [(40 ) CS 11 (173): cfa1: [ 112775144 - 119848336 ] 707319224 hsa19: [ 50887772 - 40849556 ] ] 119848336, 119880560 [(55 ) CS 12 (33): cfa1: [ 119880560 - 121690672 ] 181011218 hsa19: [ 40824500 - 38556448 ] ] 121690672, 121820640 [(16 ) CS 13 (22): cfa1: [ 121820640 - 124798800 ] 29781607 hsa19: [ 38391408 - 34709332 ] ] 124798800, Telomere [(-NA-) __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
