Related
Question:
Below can find the data, codes, and output. How can I get the plot sorted according to the original order in the row names, instead of the alphabetic order?
Note: I understand there are many existing questions and answers saying we can use as. character() and specify the orders in factor levels. However, the codes for the heatmap are a result of StackOverflow help, and I don't know how to follow those guidances to modify the codes.
Data:
data <- structure(c(1.56415665339664, 1.88622232430646, 2.7670527911374,
1.92703163268293, 1.92815958820144, 1.53777698134781, 1.85785080937573,
2.73057928043023, 1.91265151742415, 1.91392443413493, 1.51188863842245,
1.82978143286307, 2.69398442644669, 1.89813137541848, 1.89955330864784,
1.48650974337831, 1.80204149133308, 2.65731564638771, 1.88348256348151,
1.8850572392635, 1.46165563934645, 1.77465595548961, 2.62061921610674,
1.86871647260771, 1.87044727645794, 1.4373390243787, 1.74764749626653,
2.58394008473444, 1.85384450032525, 1.85573446770023, 1.41357009067807,
1.72103652543546, 2.5473217079521, 1.83887802369739, 1.84092983214118,
1.39035667008188, 1.69484124901252, 2.51080590056645, 1.82382837308096,
1.82604433604935, 1.36770438358188, 1.66907773176055, 2.47443270873259,
1.80870680674341, 1.81108886908555, 1.34561679288903, 1.64375997112343,
2.4382403018842, 1.79352448643103, 1.79607422149822, 1.3240955522757,
1.61889997899129, 2.40226488416451, 1.77829245397141, 1.78101106231344,
1.30314055915218, 1.59450786977524, 2.36654062491116, 1.76302160898326,
1.7659099185843, 1.28275010205015, 1.57059195336362, 2.33109960753632,
1.74772268775752, 1.7507811557552, 1.26292100489292, 1.54715883163515,
2.29597179596148, 1.73240624336288, 1.73563495918782, 1.24364876662612,
1.52421349731482, 2.26118501761437, 1.71708262702006, 1.72048131688628,
1.22492769546226, 1.50175943407199, 2.22676496187464, 1.70176197077873,
1.70533000345043, 1.20675103715734, 1.47979871687523, 2.19273519276265,
1.68645417152199, 1.69019056527759, 1.18911109688617, 1.45833211173251,
2.15911717460239, 1.67116887631423, 1.67507230702478, 1.17199935441569,
1.43735917405612, 2.12593030935165, 1.65591546909929, 1.65998427933569,
1.15540657239252, 1.41687834499844, 2.09319198427861, 1.64070305874773,
1.6449352678292, 1.13932289766278, 1.39688704520578, 2.06091762867183,
1.6255404684443, 1.62993378333913, 1.12373795562949, 1.37738176553192,
2.02912077829616, 1.6104362263992, 1.61498805338851, 1.10864093772681,
1.35835815434096, 1.99781314634913, 1.5953985578604, 1.60010601487567,
1.09402068215196, 1.33981110110915, 1.96700469972745, 1.58043537839775,
1.58529530794391, 1.07986574804581, 1.32173481610856, 1.93670373947838,
1.56555428842454, 1.57056327100153, 1.06616448335306, 1.30412290602036,
1.90691698438448, 1.55076256891683, 1.55591693685494, 1.05290508662397,
1.28696844538424, 1.87764965670914, 1.53606717828693, 1.54136302991323,
1.04007566304237, 1.27026404384093, 1.84890556921326, 1.52147475036339,
1.52690796441975, 1.02766427498045, 1.25400190916953, 1.82068721263795,
1.50699159342693, 1.51255784366333, 1.01565898739092, 1.23817390615927,
1.79299584293262, 1.49262369024903, 1.49831846011925, 1.0040479083513,
1.22277161138752, 1.76583156759166, 1.47837669907795, 1.48419529646884,
0.992819225075914, 1.2077863640027, 1.73919343054343, 1.46425595551545,
1.47019352744474, 0.98196123570733, 1.19320931263264, 1.71307949511363,
1.45026647522637, 1.45631802244833, 0.9714623771931, 1.17903145855638,
1.68748692465827, 1.43641295742272, 1.4425733488853, 0.96131124954466,
1.16524369529089, 1.66241206053094, 1.42269978906386, 1.42896377616518,
0.951496636765001, 1.15183684475424, 1.63785049711299, 1.40913104971445,
1.41549328031103, 0.942007524719625, 1.13880169017355, 1.61379715369429,
1.39571051700261, 1.40216554912609, 0.932833116212497, 1.12612900591028,
1.59024634304597, 1.38244167262157, 1.38898398786489, 0.923962843515001,
1.11380958437724, 1.56719183657471, 1.36932770881932, 1.37595172535787,
0.915386378581994, 1.10183426022178, 1.54462692599157, 1.35637153532236,
1.3630716205394, 0.907093641174831, 1.09019393194762, 1.52254448146634,
1.34357578664111, 1.35034626933109, 0.899074805097146, 1.07887958114475,
1.50093700627199, 1.33094282970668, 1.33777801183386, 0.89132030273522,
1.06788228949262, 1.47979668795232, 1.31847477179054, 1.32536893978401,
0.883820828081197, 1.05719325369659, 1.45911544607085, 1.30617346866088,
1.31312090423067, 0.876567338404172, 1.04680379851168, 1.4388849766195,
1.29404053293171, 1.30103552339406, 0.869551054721558, 1.03670538800144,
1.41909679318253, 1.28207734256312, 1.28911419066604, 0.862763461210944,
1.02688963517278, 1.39974226496544, 1.2702850494734, 1.27735808271664,
0.856196303691171, 1.01734831012073, 1.38081265180883, 1.25866458822613,
1.26576816767268, 0.849841587290416, 1.00807334680985, 1.36229913631591,
1.24721668475813, 1.25434521333641, 0.843691573408777, 0.99905684861186,
1.34419285322773, 1.23594186511587, 1.24308979541463, 0.837738776073223,
0.990291092711936, 1.3264849161844, 1.22484046417117, 1.23200230573094,
0.831975957773722, 0.981768533488905, 1.3091664420122, 1.21391263428855,
1.22108296039557, 0.826396124860962, 0.973481804967925, 1.29222857267723,
1.20315835391962, 1.21033180790988, 0.820992522578248, 0.965423722437341,
1.27566249504498, 1.19257743610158, 1.1997487371841, 0.815758629792936,
0.957587283315014, 1.25945945858363, 1.18216953683933, 1.18933348544938,
0.810688153486053, 0.949965667343234, 1.24361079114564, 1.17193416335278,
1.17908564604669, 0.805775023052616, 0.942552236185361, 1.22810791295855,
1.1618706821728, 1.16900467607732, 0.801013384459458, 0.935340532491705,
1.21294234895183, 1.15197832707131, 1.15908990390117, 0.796397594302229,
0.928324278496812, 1.19810573954163, 1.14225620681276, 1.14934053647087,
0.791922213798429, 0.921497374205249, 1.18358984999023, 1.13270331271594,
1.13975566649139, 0.787582002749033, 0.914853895218226, 1.16938657845174,
1.12331852601676, 1.13033427939609, 0.783371913497309, 0.908388090248925,
1.15548796281005, 1.11410062502406, 1.12107526013177, 0.779287084909823,
0.902094378370215, 1.14188618640925, 1.10504829206212, 1.11197739974653,
0.775322836401382, 0.895967346034515, 1.12857358277151, 1.09616012019459,
1.10303940177539, 0.771474662022726, 0.89000174390193, 1.11554263939148,
1.08743461972608, 1.09425988841997, 0.767738224627101, 0.884192483509395,
1.10278600069095, 1.07887022447868, 1.08563740651938, 0.764109350129453,
0.878534633810428, 1.09029647021224, 1.07046529784161, 1.07717043331055,
0.760584021869845, 0.873023417612164, 1.07806701212338, 1.0622181385934,
1.06885738197716, 0.757158375090736, 0.867654207933682, 1.06609075210328,
1.05412698649671, 1.06069660698705, 0.753828691536061, 0.862422524307139,
1.05436097767014, 1.04619002766682, 1.05268640921891, 0.750591394178479,
0.857324029040955, 1.04287113801187, 1.03840539971537, 1.04482504087941,
0.747443042079812, 0.852354523462208, 1.03161484337265, 1.03077119667177,
1.03711071021304, 0.744380325388455, 0.847509944153447, 1.02058586404586,
1.02328547368523, 1.02954158600691, 0.741400060476472, 0.84278635919742,
1.00977812901932, 1.01594625151082, 1.02211580189372, 0.738499185218138,
0.838179964441562, 0.999185724315382, 1.00875152078335, 1.0148314604562,
0.735674754410845, 0.833687079792661, 0.988802891064574, 1.00169924608359,
1.00768663713703, 0.732923935338581, 0.829304145550759, 0.978624023348409,
0.994787369801221, 1.00067938395813, 0.730244003477527, 0.82502771879015,
0.968643665843796, 0.988013815799624, 0.993807733054007, 0.727632338342786,
0.820854469794238, 0.958856511298594, 0.981376492887521, 0.987069700023661,
0.725086419474772, 0.816781178550003, 0.949257397865192, 0.974873298102964,
0.980463287105968, 0.722603822563399, 0.812804731306956, 0.939841306316456,
0.968502119815104, 0.973986486183567, 0.72018221570783, 0.808922117204624,
0.930603357166104, 0.96226084064944, 0.967637281620371, 0.717819355809315,
0.805130424971885, 0.921538807713388, 0.956147340242261, 0.961413652937977,
0.715513085094331, 0.801426839700833, 0.912643049029987, 0.950159497830092,
0.955313577336244, 0.713261327765121, 0.797808639697236, 0.903911602905182,
0.944295194679964, 0.94933503206341, 0.711062086774516, 0.794273193409158,
0.895340118763694, 0.938552316366345, 0.943475996641087, 0.708913440721818,
0.79081795643482, 0.886924370569017, 0.932928754900547, 0.937734454949501,
0.706813540866439, 0.787440468610379, 0.878660253723646, 0.927422410718391,
0.932108397178287, 0.704760608255926, 0.784138351177934, 0.870543781976314,
0.922031194531849, 0.926595821648137, 0.702752930964936, 0.780909304033727,
0.862571084345134, 0.916753029050316, 0.921194736508511, 0.700788861441736,
0.777751103056251, 0.854738402064476, 0.911585850577083, 0.915903161316576,
0.698866813958788, 0.77466159751369, 0.847042085562406, 0.906527610486473,
0.910719128502428, 0.696985262163986, 0.771638707549953, 0.839478591474624,
0.901576276586997, 0.905640684725585, 0.695142736729139, 0.768680421748302,
0.832044479699982, 0.896729834375788, 0.90066589212762, 0.693337823092351,
0.765784794771495, 0.824736410501997, 0.891986288189406, 0.895792829485681,
0.691569159290967, 0.762949945077154, 0.817551141659994, 0.887343662256013,
0.891019593271567, 0.689835433881839, 0.760174052707016, 0.810485525673035,
0.882800001653769, 0.886344298620863, 0.6881353839457, 0.757455357148564,
0.803536507019115, 0.878353373180141, 0.88176508021654, 0.686467793172536,
0.754792155267527, 0.796701119471732, 0.874001866136701, 0.877280093091286,
0.684831490024889, 0.7521827993096, 0.789976483475408, 0.869743593033824,
0.872887513352707, 0.683225345976117, 0.749625694969754, 0.783359803581409,
0.865576690219552, 0.868585538835388, 0.681648273820716, 0.74711929952741,
0.776848365944517, 0.86149931843675, 0.864372389683703, 0.680099226053889,
0.744662120045772, 0.770439535881419, 0.857509663312525, 0.860246308869087,
0.678577193317622, 0.742252711633559, 0.764130755491001, 0.853605935783729,
0.856205562645381, 0.677081202910631, 0.739889675767395, 0.75791954133658,
0.849786372462227, 0.852248440945717, 0.67561031735961, 0.737571658673098,
0.751803482189895, 0.846049235943467, 0.84837325772426, 0.674163633049295,
0.735297349764132, 0.745780236836504, 0.842392815061733, 0.844578351246036,
0.67274027890897, 0.733065480135467, 0.739847531942042, 0.838815425095341,
0.840862084327889, 0.671339415153075, 0.730874821111148, 0.734003159978691,
0.835315407924883, 0.837222844533538, 0.66996023207371, 0.728724182843864,
0.728244977211044, 0.831891132147497, 0.83365904432554, 0.66860194888288,
0.726612412964853, 0.72257090174049, 0.828540993150013, 0.830169121176865,
0.667263812602402, 0.724538395282482, 0.716978911607115, 0.825263413143681,
0.826751537644663, 0.6659450969995, 0.722501048527904, 0.71146704294807,
0.822056841163083, 0.823404781408691, 0.664645101566158, 0.720499325146197,
0.706033388211285, 0.818919753031688, 0.820127365276744, 0.663363150540402,
0.718532210131447, 0.700676094423347, 0.815850651296407, 0.81691782715934,
0.662098591967731, 0.716598719904252, 0.695393361510352, 0.812848065133384,
0.813774730015795, 0.660850796801012, 0.714697901230176, 0.690183440670466,
0.809910550227147, 0.810696661773713, 0.659619158037194, 0.712828830177725,
0.685044632796959, 0.807036688625131, 0.807682235223833, 0.65840308988929,
0.710990611114427, 0.679975286950419, 0.804225088569508, 0.804730087892066,
0.657202026992117, 0.709182375739667, 0.674973798878868, 0.801474384308113,
0.801838881890468, 0.656015423640367, 0.707403282152963, 0.670038609584492,
0.798783235886211, 0.799007303748793, 0.654842753057613, 0.705652513956377,
0.665168203935704, 0.796150328920716, 0.796234064228217, 0.653683506694951,
0.703929279389845, 0.66036110932326, 0.793574374358407, 0.793517898118688,
0.652537193557989, 0.702232810498201, 0.655615894359165, 0.791054108219602,
0.79085756402133, 0.651403339560998, 0.700562362328734, 0.650931167617108,
0.788588291328647, 0.788251844117223, 0.650281486907035, 0.698917212158154,
0.646305576413202, 0.786175709032535, 0.785699543923806, 0.649171193492958,
0.697296658747865, 0.641737805625809, 0.78381517090886, 0.783199492040099,
0.64807203233824, 0.695700021626499, 0.63722657655324, 0.781505510464259,
0.780750539881858, 0.646983591036597, 0.694126640398684, 0.632770645808186,
0.779245584824428, 0.778351561407705, 0.645905471229425, 0.692575874079065,
0.628368804247696, 0.777034274416723, 0.776001452837246, 0.644837288100137,
0.691047100450625, 0.624019875937613, 0.7748704826463, 0.773699132362083,
0.643778669888508, 0.689539715446396, 0.619722717150336, 0.772753135566687,
0.771443539850621, 0.642729257424162, 0.688053132553661, 0.615476215394868,
0.77068118154563, 0.769233636547476, 0.641688703678407, 0.686586782239811,
0.611279288478086, 0.768653590926998, 0.767068404768254, 0.640656673333622,
0.685140111399023, 0.607130883596222, 0.766669355689472, 0.764946847590445,
0.63963284236946, 0.683712582818977, 0.603029976455569, 0.764727489102711,
0.76286798854108, 0.638616897665165, 0.682303674666839, 0.598975570421445,
0.762827025381632, 0.760830871281802, 0.637608536617296, 0.680912879993784,
0.594966695694488, 0.760967019339396, 0.758834559291943, 0.636607466772238,
0.67953970625735, 0.591002408513372, 0.759146546039657, 0.75687813555014,
0.635613405472859, 0.678183674860939, 0.587081790383065, 0.757364700448588,
0.754960702215028, 0.634626079518729, 0.67684432070982, 0.583203947327774,
0.755620597087154, 0.753081380305464, 0.633645224839319, 0.675521191782987,
0.579368009167761, 0.75391336968409, 0.751239309380749, 0.632670586179664,
0.674213848720291, 0.575573128819215, 0.752242170829969, 0.749433647221233,
0.631701916797938, 0.672921864424234, 0.571818481616423, 0.750606171632756,
0.747663569509711, 0.630738978174469, 0.671644823675892, 0.568103264655477,
0.749004561375194, 0.745928269513948, 0.629781539731715, 0.67038232276441,
0.564426696158795, 0.747436547174325, 0.74422695777065, 0.62882937856474,
0.66913396912956, 0.560788014859765, 0.745901353643456, 0.742558861771208,
0.627882279181762, 0.667899381016863, 0.557186479406819, 0.74439822255683,
0.740923225649459, 0.626940033254362, 0.666678187144791, 0.553621367786292,
0.742926412517243, 0.73931930987174, 0.626002439376954, 0.665470026383597,
0.550091976763438, 0.741485198626834, 0.737746390929453, 0.62506930283514,
0.664274547445332, 0.54659762134097, 0.740073872161253, 0.736203761034359
), dim = c(5L, 151L), dimnames = list(c("highest_pi_count", "second_pi_count",
"highest_si_count", "max_overhangz_minus", "max_overhangz_plus"
), NULL))
Codes:
library(latticeExtra)
library(ggplot2)
library(tidyverse)
rownames(data)
as.data.frame(data) |>
rownames_to_column() %>%
pivot_longer(-rowname) %>%
ggplot(aes(factor(name, unique(name)), rowname, fill = value)) +
scale_x_discrete(labels = ~., breaks = ~ c(6:29/10)) + ggtitle("Title") + theme(plot.title = element_text(hjust = 0.5)) +
geom_tile() +
scale_fill_gradientn(colours = c("blue4", "white", "red3")) +
scale_y_discrete(position = "right") +
theme(legend.position = "none",
axis.ticks = element_blank(),
axis.title = element_blank(),
text = element_text(face = 2))
Output:
> rownames(data)
[1] "highest_pi_count" "second_pi_count" "highest_si_count" "max_overhangz_minus" "max_overhangz_plus"
To fix your issue convert row name to factor with the order of the levels set according to the row names of your original dataset:
as.data.frame(data) |>
rownames_to_column() |>
pivot_longer(-rowname) |>
mutate(rowname = factor(rowname, rownames(data))) |>
ggplot(aes(factor(name, unique(name)), rowname, fill = value)) +
scale_x_discrete(labels = ~., breaks = ~ c(6:29 / 10)) +
ggtitle("Title") +
theme(plot.title = element_text(hjust = 0.5)) +
geom_tile() +
scale_fill_gradientn(colours = c("blue4", "white", "red3")) +
scale_y_discrete(position = "right") +
theme(
legend.position = "none",
axis.ticks = element_blank(),
axis.title = element_blank(),
text = element_text(face = 2)
)
I've run a PCA using prcomp in R and I am trying to produce a variable plot that has 1) a subset of the variables (arrows) in a different color (black) than the rest of the variables, 2) sort those variables prior to plotting so the black arrows aren't covered up by any of the other arrows, and 3) label the black arrows with their TUXXXX number.
Here is a truncated version of my data:
structure(list(sdev = c(21.7106794138444, 15.6885074594869, 11.9124316528111,
10.155277241318, 9.31528828036412, 7.56876266263865, 7.19938201515987,
5.81620977435434, 5.00785424840699, 4.57228787327195, 4.51525494575488,
3.20601607034873, 2.91477640215067, 2.48737967730048, 2.13230488376163,
1.74923754200417, 1.32745772948038, 1.27373216417502, 0.924437474777366,
0.749074623004602, 2.499709597053e-15), rotation = structure(c(-0.0710441966458092,
0.091894514828866, -0.0892433537986534, -0.269473709517009, -0.270455466278075,
0.217492458575054, 0.104541973199297, 0.198858094257877, 0.0222680112919805,
-0.220704163347643, -0.0144913885562279, 0.191255085890651, -0.0639203495167002,
0.156262929972648, 0.184067836594737, -0.221797618857792, 0.152932751853774,
0.218733634550932, 0.213989867048259, -0.219209176045661, 0.162232843055052,
-0.182649403408074, 0.0732706002149817, -0.0636694683826245,
0.197950156891487, -0.201332798849071, -0.0908133550739104, 0.204561640782066,
0.144543461720936, 0.128722984809536, 0.152014925749702, -0.219186163224382,
0.137427869376807, -0.0385728465064452, -0.210836096156869, 0.295086522896984,
0.0791956927556792, 0.154250304318938, 0.212537035577304, 0.165554637810987,
0.0559168852882808, 0.0389301896007549, 8.50109782205011e-06,
0.0386719091643949, -0.223181316780788, 0.0201295766763838, -0.0993945550528442,
-0.0797657364873552, 0.290392094157839, -0.210104310439075, 0.190507620860241,
0.288216604566037, -0.128091623622577, 0.0523741236254804, 0.165286635170839,
0.26567738477018, 0.067304658063138, -0.185395359345644, -0.0824392053628523,
0.0232331903282599, 0.206712718290849, -0.0599454190622519, 0.0694655529225049,
0.27934713843206, 0.120105146108152, 0.276987279291797, 0.137068573563571,
0.290227550805898, 0.0299855300238849, -0.048045392789825, 0.217209812823655,
-0.245207498414819, -0.0969050397538762, 0.0152958298727508,
0.0611874790392221, 0.244015616675892, -0.341439467152496, -0.0574074669754257,
-0.187473105527714, 0.0859725296785313, -0.230526618631008, 0.164238788534203,
-0.422116092865793, -0.0847447867791298, 0.0984005915083912,
0.0400031956936472, -0.0768787993932683, -0.127305543959199,
0.174417860300221, 0.00642366483942508, -0.0644270895895563,
-0.0445557418009, -0.206781215165792, -0.398796252833277, 0.198729650467593,
0.0351312161529938, -0.149474677035792, -0.0671760700655275,
-0.0776573448858521, -0.0276268199179851, -0.0622850132226139,
-0.0133848915206553, -0.0412751078511757, -0.196695235499305,
0.035533391699955, -0.0471461480863407, -0.471887698745297, 0.385807939659813,
-0.0521508210675749, -0.106054753553622, -0.0520098951237357,
-0.328301752147492, -0.0885218905714762, -0.392528395483506,
-0.0535243088811493, -0.0012232141956296, -0.143968515821013,
0.110020791115645, 0.105789408005642, -0.0173550741395983, -0.0754711706175302,
0.105789800532602, 0.035696880367874, -0.0235700328573188, -0.086664772131814,
0.140659513401983, -0.147738109860937, 0.158253547879328, 0.096508294839766,
-0.0405645836586841, -0.0308947063294566, 0.196563611460491,
-0.163272613036823, 0.102159062108236, -0.242801481572191, 0.100060947435785,
-0.122408620775568, 0.0879253395566385, 0.0964448760240853, -0.10574290240921,
0.473324985332234, 0.237791221156891, 0.121312056042028, -0.0156952213549262,
-0.0509545579371996, -0.188663582956394, -0.0296470996543493,
0.259217775431442, -0.151966417414955, -0.0246227447960926, 0.412107323913548,
0.135772975480753, -0.0902168853917747, -0.0478519874649824,
-0.0293395316078502, 0.036873372126565, -0.0600803315985544,
0.0951569802806819, -0.174539345831377, -0.0567399129758035,
-0.0522218071944251, -0.0208963041931566, 0.15376867866167, -0.0906491458699936,
-0.13257887196968, -0.0137277293203767, -0.329880055323168, -0.227898895128165,
-0.0402178711138279, -0.191724054954671, -0.0457418642155951,
-0.113122585201102, -0.0523451728177034, -0.0619052030873901,
-0.232037793152732, -0.286608825250971, -0.337575786471338, -0.106930203013502,
-0.0488020861881439, -0.0342694484533206, -0.316010730818033,
-0.138362452107825, 0.000204654351180751, 0.34062110315846, -0.0174472565351551,
0.321216705099307, 0.109402271713093, -0.195953316697484, 0.106147314965255,
0.0965509971711528, -0.0619150877056984, 0.0823900429544073,
0.107331603831664, -0.196505739977505, -0.106531330626324, 0.083316396319378,
-0.0150328553267955, -0.131557865292901, -0.197676924192316,
-0.302438247351333, 0.0850937235677142, -0.0588088540505814,
-0.0137660099545295, 0.0981614485249126, 0.064387690302654, 0.0870877447014138,
0.115369776587117, 0.0743160202050587, -0.302341999379974, 0.0990657432090211,
-0.109416851119016, 0.0566207941620112, -0.563376445425201, 0.022987878311214,
0.00914995477596945, -0.0698458973541682, -0.339551294424273,
0.328653723240687, -0.311621856766325, 0.0451475830297207, 0.0486778657760599,
-0.203890915518391, -0.00733314039012801, 0.0302380986964549,
-0.187740946530342, 0.107505967223115, 0.0112058761308914, 0.0124189091072754,
-0.0946013978939246, -0.0108516955206006, -0.0224544300789198,
-0.200447923165361, 0.00227210575436315, -0.0230219082173417,
-0.115276809069961, 0.113599893064251, 0.0547021334782484, 0.0750295713780988,
0.0318176533771422, 0.270756307602172, 0.0158888930662066, -0.0321100878211108,
0.0137635276916354, 0.285629078757821, -0.0707869220022327, -0.196194173921016,
0.0691894805409008, 0.196425060869134, 0.0609662656649608, 0.0112726648905893,
0.0763548314916228, -0.0203154117969338, -0.154183879477363,
-0.263435512251247, 0.0904584427059006, 0.443081105454767, 0.145760707944433,
-0.124076693641366, -0.00926033658765893, 0.322008337751576,
0.038976848601001, -0.00243771729397693, 0.0682300992336678,
0.101208051579165, -0.0044801862648156, -0.0288479608546419,
0.338293188340527, -0.148520145256712, -0.114001067101895, 0.0429068793268258,
-0.0808048267081677, 0.0147906601414549, 0.00024012812721683,
-0.0245592922986743, 0.26189600932665, 0.0089069926197661, 0.102640316156789,
0.0293398861736338, 0.461748240785653, 0.0268162505216955, -0.0262617915066989,
0.128477574136056, 0.304269602566599, -0.00833894425455883, 0.0183612622798514,
-0.0694598653855358, 0.143473468921974, 0.494741840190451, 0.11710470731812,
0.0938975501248067, -0.111819356177393, 0.074487393087207, -0.236744656433776,
-0.00888263838265717, -0.0411308668141714, -0.0205247425377322,
0.00682064532990936, -0.0657589971316789, -0.0854778766938379,
-0.0446898232128541, 0.0595158665026168, -0.164975323219787,
-0.21297106776461, -0.24328292167101, -0.0421999896796148, -0.0454152302659223,
0.536798625888851, -0.107184764029374, -0.134710620292912, -0.0486878301757869,
0.0112145455041351, -0.0888773206829193, -0.0231311387150055,
0.160224502964451, 0.116642580196556, 0.132842779999407, 0.0815222556677455,
0.0846715273158058, -0.168358316697627, -0.141086214309714, -0.0812000194350307,
0.108662664398591, -0.0391921007274459, -0.231315089883372, 0.25650378734853,
-0.0459946623780035, -0.366253502200094, 0.0485143397750607,
0.0417032700002287, -0.133038776710167, -0.230298287475923, 0.0369759420468897,
-0.00906710613565646, 0.0342935061381997, -0.273232725280753,
0.0602621389601039, 0.206808600842623, -0.32050527301245, 0.0622454407253566,
-0.0606525149750042, -0.0660628546307464, -0.038482743495315,
-0.198081930190034, -0.0732284456644516, 0.405606839140943, 0.0459681574877996,
-0.0809991829793473, -0.0121246261443008, -0.326485629157016,
0.0750251561011149, -0.0511696428776832, -0.0776271119484428,
0.0114894367280979, 0.014021492503364, 0.118849898460583, -0.0350044617639092,
0.0269961041673681, 0.239744113644073, -0.233879279575839, -0.115473220682409,
-0.139126309932733, 0.11002945246376, -0.0830577400425309, -0.0642289884705754,
-0.0817307298240943, 0.0291046627888711, -0.0466920006759583,
0.256672775663942, -0.0137954035992047, 0.17957481447807, -0.0303657675236596,
0.518263963411537, 0.344005787563606, -0.0655688570898937, -0.0265294401946153,
-0.0263107898868621, 0.289506920843548, 0.118475732210224, 0.053567958881615,
0.146877133183997, -0.0549414997788912, 0.143851895777107, 0.0428308701652828,
-0.348959761946987, -0.189096210288467, -0.00385035890156567,
0.150120708161381, 0.0607203081654835, 0.000166036942981724,
-0.0397011322296128, 0.0239760961564877, -0.00913095901818713,
-0.136862964262375, 0.102111106897839, 0.0687169707815706, 0.239013538410316,
-0.39892902681439, -0.197510750116376, 0.0499804963638917, 0.128640553281137,
0.110040614593457, -0.0490679379082568, -0.234025706976098, 0.157764576030734,
0.0571664700394519, -0.074307378291659, -0.231764546118876, 0.565447067892318,
-0.21480675567977, 0.0854469138131175, -0.0976097155414531, 0.160801495571792,
-0.0624041394929365, -0.0127140410117323, 0.196633195933416,
0.0340211274197625, 0.0195595932420065, 0.0964345475133984, -0.174817763577626,
0.118153342593374, 0.225094093283825, 0.0276798513722757, -0.39244132040106,
0.0738620568768489, 0.027031867335908, 0.0584738274201628, 0.00728266929277696,
-0.200928813049723, 0.100479151028423, 0.217133686900663, 0.152350478604359,
0.0798552081200005, 0.00566909653422267, 0.15071802028838, -0.0828967915619975,
0.0689991484295119, 0.118221935665031, 0.418906970577945, 0.161468821548893,
-0.0443917405359232, 0.127346300562372, 0.0854487046042784, 0.0086243432217815,
0.0186528057922736, 0.143224492276055, -0.113522100459258, 0.365748845830972,
-0.353026210924788, 0.0928434654809545, -0.217256921150109, 0.0926292235900054,
-0.0360081503134767, -0.0761247497474667, 0.0137710983238269,
0.135604749840715, -0.366785665744512, 0.174928851514231, 0.0364877873595186,
0.0808239023932249, 0.057715463934493, 0.103601691734105, -0.109666857605533,
0.38393023627746, -0.117248332196572, 0.134270206368149, 0.0580313640704774,
0.0649484539445183, 0.215353645432973, 0.0977108936702795, 0.16612143102929,
0.152770485581267, -0.0516344648811396, -0.0620940862015108,
0.114473317921704, -0.0383202829927206, 0.0820220075422196, 0.0556449561108403,
0.0053103511581772, 0.141569483443786, 0.264761186392854, 0.231379833347928,
-0.186458908266326, -0.280612868599116, -0.0134977078164976,
-0.190682413873154, 0.103622962625572, 0.048589128670959, 0.189228510250879,
-0.267141675131746, 0.24968486643315, 0.000515248155113446, 0.0518714761513495,
-0.0123619986636973, -0.0400774164248301, 0.14696585867848, 0.00893605318819822,
-0.115310054522583, 0.0730556770081751, -0.0415182843892857,
-0.100197819082248, 0.366565649542296, -0.280217314004845, 0.0742548639941826,
-0.0598573375691124, 0.0773367591528449, 0.00440008401952238,
-0.133823453930746, 0.0129306082883093, 0.249237874757721, 0.178566150191566,
0.0236542331383419, -0.121983928027758, -0.407639692892733, 0.105625218150781,
0.0606170741355181, -0.205658012407281, 0.198411895098939, 0.314706002474467,
-0.0699572649574542, -0.256054196724344, -0.0597932960838045,
0.12552805883941, -0.0738282459049305, -0.270044343022259, -0.241030766471284,
0.101198016298546, 0.207020228203116, -0.105134856921501, -0.0677078010226698,
0.0597983546937891, -0.181549492585234, -0.0707807644693492,
-0.242992185702685, -0.164512180143508, 0.0345701932987703, -0.0663154906257285,
0.302765631613708, 0.118224308922645, 0.0673423264385151, 0.417740310978627,
0.111554431260226, 0.16022442648127, -0.25635936207313, -0.35265959952747,
0.0429263261414316, 0.140480243375922, -0.130382613873709, -0.0446328161712632,
0.120679204177165, 0.100894280068327, 0.102396715497024, 0.292226072514004,
0.0170232777619504, -0.0854499461046589, 0.150649869168952, 0.054457847943396,
-0.0755734422376126, -0.0485946974863497, -0.214957242429833,
-0.19078567042026, -0.265677083691519, 0.11169177584997, -0.0546945656068388,
-0.133918556644411, -0.0815433043760734, 0.0419081230722499,
-0.119806690064691, -0.129001645490559, 0.164330559962565, -0.494999436348284,
-0.0104711534954494, -0.151933262861035, -0.0616048245349729,
0.0510091356535412, -0.15286661063535, 0.286795952380893, -0.00227878893487643,
-0.121875652782209, -0.0463279998907995, -0.08440989195395, 0.0161463685432679,
0.283239464786081, -0.0997334879966379, -0.107077349937415, 0.0797083319878824,
-0.023817698507559, 0.0308411589066249, -0.33639040110823, -0.22283646210506,
-0.0472382156514798, -0.0383954704026199, -0.279873414947505,
-0.110928067466998, 0.0224738031700519, -0.240448693328791, -0.0234165327843327,
-0.164995870598194, 0.0704836439566211, -0.0341351365684308,
0.253109207656731, 0.0986224720533927, 0.0788051554263656, 0.0348921204456385,
0.179396594950779, 0.183681460022634, 0.165162541829685, -0.131056895026594,
-0.043658366747374, -0.164919051972395, 0.000438005379084683,
-0.203731301273092, 0.156795003278361, 0.147134057978626, 0.0674666995224469,
0.0387098969653998, 0.0780676163151537, -0.020902824954741, 0.0863172379869736,
0.0254543096586073, 0.221894385139484, 0.605679699321179, 0.00510726652451246,
0.0122597366937914, -0.177909018628657, 0.112879199323714, 0.00722763837901889,
-0.226515396833189, -0.179516219763039, 0.0427914553668244, -0.301499857490821,
0.07077374887633, 0.0684478824696676, 0.109600278015023, -0.246938356324243,
0.0361050298642269, -0.273671368850468, 0.0361740234532937, -0.0379154357295922,
-0.140858482198813, -0.210557939874765, -0.183530140411959, -0.170246755949232,
-0.153695213839514, 0.197257212733899, 0.176328880404244, 0.267335680405474,
0.0108203611164604, 0.420718967278865, 0.0610513794750232, -0.0695142190235678,
-0.0211301631659223, -0.109598478540236, 0.0168455051311258,
0.0297211637285332, 0.0641419705460361, 0.196768178957862, -0.101223881022882,
0.156790420489079, -0.248396712859248, -0.0324362090799154, 0.116497846273851,
0.206631759837322, -0.26949738606709, 0.0523825016986687, 0.273991705979531,
-0.211717836239047, 0.343025575540513, 0.123225187210276, -0.130175293604759,
0.383814266477628, -0.0425846714760806, -0.120958839664608, -0.222129739206683,
0.0129931836329733, 0.144776049791581, -0.162517014653887, -0.0859042045925978,
-0.181054333724447, 0.0216585016380889, -0.0923250223535876,
-0.103640362754872, -0.223956862384563, 0.213416939275391, -0.21381164272001,
-0.235056303853187, -0.0630316125230271, -0.105100763021652,
-0.0360147146171049, -0.0709796908578211, 0.113552355258653,
-0.00724023466942227, -0.0325253318186826, -0.0460572605581843,
0.362655177721721, -0.102542875980008, -0.0801797191999012, 0.107198034030814,
0.183373945987406, 0.0309647928318137, -0.236970374446373, 0.0520664942182686,
-0.0402703424381589, 0.0504291632173252, -0.158130484290709,
-0.112283134517499, -0.300637215703608, -0.231320162708774, 0.0785169224262154,
0.23958352054759, 0.426485005426312, -0.0815406319434192, 0.00340039424643938,
-0.0375428689912434, 0.33351700972973, -0.0800149829632403, 0.140498383489292,
0.167368656035586, -0.0844117015281891, 0.0785085162751926, 0.0489486062002369,
-0.0561511863757523, 0.234387154088994, 0.0481142285230063, 0.0705808780719704,
0.221875388986379, -0.148031374880166, -0.0894717302805895, 0.231743942045981,
-0.210616362796393, -0.0812013243962482, -0.185771800435443,
-0.0438135376399639, 0.218374086144347, 0.121623152294345, -0.170058260806661
), .Dim = c(35L, 21L), .Dimnames = list(c("TU35976", "TU38792",
"TU35975", "TU4897", "TU9262", "TU18312", "TU38793", "TU30299",
"TU38794", "TU16422", "TU32068", "TU12311", "TU18325", "TU22117",
"TU20746", "TU19660", "TU41790", "TU43092", "TU18309", "TU19659",
"TU16277", "TU34798", "TU10953", "TU18317", "TU18307", "TU9186",
"TU17864", "TU15658", "TU14030", "TU15669", "TU9890", "TU33617",
"TU16279", "TU27042", "TU17167"), c("PC1", "PC2", "PC3", "PC4",
"PC5", "PC6", "PC7", "PC8", "PC9", "PC10", "PC11", "PC12", "PC13",
"PC14", "PC15", "PC16", "PC17", "PC18", "PC19", "PC20", "PC21"
))), center = c(TU35976 = 6.61238491831215, TU38792 = 10.2487864648085,
TU35975 = 7.41959195648095, TU4897 = 10.9311685931669, TU9262 = 11.0918083964267,
TU18312 = 11.7076238975148, TU38793 = 9.74398564779285, TU30299 = 12.892315214803,
TU38794 = 9.54139631961704, TU16422 = 12.3329627148834, TU32068 = 7.97610519995008,
TU12311 = 10.3617106799956, TU18325 = 7.59060411317826, TU22117 = 9.03799676912002,
TU20746 = 9.178102720165, TU19660 = 9.87502757868133, TU41790 = 8.9943466756264,
TU43092 = 10.9585816938684, TU18309 = 10.6298994150062, TU19659 = 10.3934726050624,
TU16277 = 8.32556725898017, TU34798 = 7.02103110645297, TU10953 = 7.01860382462013,
TU18317 = 7.36244270525821, TU18307 = 11.3778296252987, TU9186 = 10.7024172065364,
TU17864 = 8.2903853475227, TU15658 = 6.58819311715545, TU14030 = 8.134641008701,
TU15669 = 7.01449935590416, TU9890 = 8.13988320730633, TU33617 = 10.6809972208427,
TU16279 = 7.9614057227348, TU27042 = 7.34221912393268, TU17167 = 11.2592403834747
), scale = FALSE, x = structure(c(24.4897868475554, 12.4416010164283,
12.3200413661576, 6.91612686789472, 30.5655787046573, 33.5373707608136,
6.7521417841798, 31.2938791704631, 29.0635289584208, -23.3354523484558,
-15.6229701200559, 6.43010789907339, -24.9147092475567, -31.0094950116289,
-28.4137928961069, -21.5926413362042, -3.97585196983663, -20.8563124842282,
-0.593920033332157, -3.64572299218967, -19.8492949360489, -9.96693484414473,
-23.5449684594658, -18.996265716943, -23.7811581634436, 3.93738430048,
9.37353195322868, -19.423648038774, 7.87596964384287, -3.13661328005585,
-7.09016527838172, 3.71259770410569, 23.247076153704, -2.03249135782906,
-3.64325919770972, -6.94450353973124, 7.55306940052889, 31.346240781817,
-3.48931847017886, 29.4056029647248, 12.4861765573918, -6.8883231131661,
-6.8786118382072, 15.8325477899497, 15.3870141654559, -11.3237572934353,
-4.49560209730399, 1.115062768913, -20.2141393532659, 9.76193925462921,
-1.15731242202584, 6.22677163564855, 14.431909276018, -13.4812455995592,
-9.58096578254001, -10.0824168036361, -10.1417065317313, 6.11246229818054,
-4.84525793168699, 14.9352609583166, 16.1272133946181, -17.0836840355263,
9.35451814718846, -6.63382933222961, 5.79244304277846, 7.73028939318989,
-6.00113726185851, -8.57575275731579, -6.70414747438031, 11.0021295692077,
7.62205208005063, -6.03281844054887, -1.59444275830952, -7.70227835025931,
-6.4184560830441, -12.3696585650364, 1.42826497938592, 8.28315195406988,
-13.8464231033084, -6.95076374785849, 5.03179685038719, 16.6482451825695,
25.4557137741208, -6.16437895161065, -0.188023057317339, 3.51919999146276,
-3.06021456766217, 12.7300628490286, 0.340348842202129, -4.18342952961951,
9.91999265542052, -6.2673275997563, -6.19445525060979, -4.89432563084305,
13.3155996092707, -17.5672684031691, -11.0570104414685, -10.1512940397755,
-2.83422832939901, 11.2933772476208, 18.438427176494, -10.6101681685668,
2.08048057352476, 5.50745324586983, -0.13719717270706, -8.9884667602456,
-2.05928444313446, -0.305773393272364, 4.86554463593651, 3.76285186780856,
6.3828926870771, -9.3655246100571, 7.43651870615979, -1.67082019758757,
10.9895599029683, -13.9919213196391, -7.47010213808747, -5.18316808251628,
3.09935494483666, 0.0636557688816208, -1.54858827353552, 5.76445544260896,
-7.0244817763404, -6.91485718432494, 6.12343812013661, 16.0347161023267,
-6.04452306419056, 0.503323158087318, 1.11905480618314, -3.32310396567328,
3.1274151532053, -1.0139189902977, -2.77210893619864, 5.73300675328722,
-2.20016077060467, 6.88733384115661, 9.29615931028018, 7.89283806484175,
-0.0231784873612931, -13.6472361241707, 8.57047305874539, 9.82836851225084,
-12.0042402273719, -9.50514572716062, -7.91330903314343, 8.52172081325052,
-3.03276814511544, -8.49580406093405, -2.33842120098311, -0.32804627167184,
-3.21259723306479, 11.1229701308803, 0.598883041283415, 5.30460923308376,
1.00174438392002, -1.97289932239362, -10.2991918599681, -3.87835358264123,
-2.25659286308468, -0.166946108131337, -4.78949081440092, 3.19183354598423,
4.8383258847991, 0.890821633270318, 14.2375666218091, -6.1498396420083,
0.00274092818159516, 2.69868755607018, -6.72994548395586, -2.87563409646143,
6.39162216555694, -1.86430696261611, 1.49349337044142, -0.729377873822674,
3.91976580406387, -4.55169784656593, 4.01470477896383, 9.45815735561592,
3.87143750044459, 2.8485319829077, -3.54554995799273, 2.1592779575771,
-10.2265407210419, -2.57734842434015, 3.85728531842232, 3.91861430861892,
-6.08938026261039, 2.90775931244037, -5.65086822564584, -7.43442226552497,
-0.28834246695089, 4.46777893840091, 1.61867182502613, -1.24379496084697,
3.09330484415506, 4.6129428255061, 4.21327794340158, -2.73480960596815,
-1.34501116390485, -4.95417459779772, -1.50853746838879, 4.27315750243976,
5.82184171959861, -3.417633941994, 9.09574374372561, -2.51371486672449,
-5.3839624715871, 4.02549123440248, -4.25205519103745, -6.14575157593085,
3.39653326861517, 2.02195972923006, 0.737129651195039, -3.82060608464456,
-5.8916049884138, -1.48103866101421, -4.15462594283813, 2.85762680978641,
8.01848063520863, -0.33318325779775, -4.11692676832453, -3.14609004733945,
-3.87716607273847, 2.15863837910027, 5.89633795035078, 6.78408908767204,
5.31335299395838, 2.26271596099398, -6.04004269690954, 0.66943773718038,
-7.25501768327073, -1.85352333808763, 0.881321761644734, -3.37300816578494,
1.72336449307423, -1.12738959732561, -1.00315691989581, 1.22856047955961,
-5.49078567616185, 7.36566150881169, 1.5098504796481, -4.40047569826546,
1.27392003732061, -1.32431052177641, -3.49056795319117, -1.81317399652518,
4.21470358856118, -3.35477920804357, 0.76847614798403, 4.53289945581295,
1.08860795433414, 2.64380516830634, -3.5282442958868, 0.44695028657078,
6.40236529107031, 1.40990912221294, -0.646731159533322, -0.886819655598922,
-0.75355529529178, -3.22681748481868, 1.14443040171423, -1.01721629940821,
-1.91406441412825, 2.79795257160194, 1.42125707801834, -3.78513731974819,
5.22015934186783, -3.11038146921718, 4.2502371748613, -2.23573426796635,
1.2890326250481, -3.2044887828359, -0.0731034485321765, 1.00689042463065,
3.10770789167095, 0.653717111815948, 2.0046773828616, -5.05502285926648,
0.670269628945531, -0.0735066911646989, 0.567378074513561, -4.12216401435608,
-1.34361612453716, -2.31031886740573, 4.81653041064168, 0.0225272846463933,
-2.29390683891735, -3.34883149458521, 1.62567523985947, 1.46413728406672,
2.19861127383831, -2.2813563892377, 1.48246659900362, 1.20813467297598,
4.19016124903691, -3.8965814643401, 3.56184752348759, -1.3878373295987,
1.21458977444305, -0.654562661488364, 0.95632241225039, -1.55232118983605,
-2.44238882925152, 1.13264815277383, -1.78862294430386, 0.991256365490394,
-3.88951966260596, 0.202107196565979, 0.551090451321218, 2.66348016156913,
-0.854309938495967, -0.232716071672556, 0.851674118269924, -0.985653498045505,
1.36933618443016, -0.586856938122795, -4.1119301713715, 1.0178280457984,
0.338189042904915, -3.94726389810293, 1.81377558126379, 1.58372328030388,
1.72308989346244, 2.19193280022478, -1.61256462530986, 1.96795000825428,
-0.618527965915178, 0.0756192601567935, -0.659358025397941, 0.0429550715360451,
-0.381324238450748, -0.652905276182954, 0.170764900323459, -0.641706794178386,
-0.0221222785810836, 2.30873232738461, 1.27649223277723, 0.0239704300112933,
1.61258498499148, -2.05589657580369, 0.0265571467877399, 1.25029287962319,
-0.571431828720408, -1.73619374916576, -0.0382088144200337, -0.321239936961026,
-0.611330806517872, -2.02090803294801, 3.14396363436869, -1.10262030332252,
-1.04168768271449, 0.37429480503427, 0.104514606536275, -0.478767952336214,
-0.192546893385304, 2.1627661920923, 0.195395664072845, -0.782469508379251,
2.08546244595102, -0.0517755897340584, -2.39907743930418, -0.211807017194547,
1.71381905614487, 0.83072521628824, -1.60469551897193, 0.245800797687136,
-1.46819802906202, 0.849407708022543, 0.917523253251796, -1.60319562494994,
1.89391093257302, 0.378604048008997, 0.14292959825461, -0.126413703865446,
-1.54645185235878, -0.749700454486876, -0.0658077625468692, 1.55140944467166,
0.302322289091245, 0.458890126236888, -0.788003883776506, -1.8522098865949,
-0.303600234099745, -0.869244048674518, 2.09285131170145, 1.13594169034256,
0.00257581110260907, 1.28821287823238, -0.603719944771219, -0.0929817189825319,
0.646083878133434, -0.997976432253114, -0.0887921104756663, -0.252335674196718,
0.885171983511843, -0.455664633944459, -0.080944547822266, -1.11153038212543,
0.684953529364659, 0.00107091275654292, -0.0523472212020604,
-0.722665903230519, 0.112695429436832, -0.431332703576418, 2.28632079906497,
0.389976510624337, -0.672770348738164, -0.710872748519771, 1.11097033239239,
-0.117525841496901, -0.38659967304807, -0.311519238694312, -0.58769206070627,
0.880912821150173, -0.223180949008615, -0.0563287542495416, 0.677577799103385,
0.174786309548343, -0.486054451272286, -0.875421020334039, -1.99840144432528e-15,
-2.19269047363468e-15, -2.80331313717852e-15, 3.39658856596259e-15,
4.44089209850063e-16, 4.44089209850063e-16, -2.44249065417534e-15,
1.11022302462516e-15, -6.66133814775094e-16, -4.44089209850063e-15,
-3.5527136788005e-15, -2.33146835171283e-15, -6.43929354282591e-15,
2.44249065417534e-15, -1.4432899320127e-15, -2.19269047363468e-15,
1.0547118733939e-15, -1.11022302462516e-15, -9.99200722162641e-16,
1.33226762955019e-15, -1.55431223447522e-15), .Dim = c(21L, 21L
), .Dimnames = list(c("EDA_01", "EDA_02", "EDA_03", "EDA_04",
"EDA_05", "EDA_06", "EDA_07", "EDA_08", "EDA_09", "NDA_01", "NDA_02",
"NDA_03", "NDA_04", "NDA_05", "NDA_06", "NDA_08", "NDA_09", "NDA_10",
"NDA_11", "NDA_12", "NDA_13"), c("PC1", "PC2", "PC3", "PC4",
"PC5", "PC6", "PC7", "PC8", "PC9", "PC10", "PC11", "PC12", "PC13",
"PC14", "PC15", "PC16", "PC17", "PC18", "PC19", "PC20", "PC21"
)))), class = "prcomp")
Here is the minimal code to reproduce this problem:
library(factoextra)
library(tidyverse)
# create a named factor for coding the coloration of the variables in the plot
markers <- facto_summarize(pca,
element = "var",
result = "contrib",
axes = c(1, 2)) %>%
mutate(candidate = ifelse(name == "TU35976" | name == "TU18317" | name == "TU12311" | name == "TU3565" | name == "TU9890" | name == "TU18316", "1", "0")) %>%
select(name, candidate)
markers_vec <- as.factor(markers$candidate)
names(markers_vec) <- markers$name
(var_cluster_PC12 <- fviz_pca_var(pca,
col.var = markers_vec,
axes = c(1,2),
select.var = list(contrib = 200),
palette = c(
"grey90",
"black"),
label = "none",
) +
theme_bw()
)
Which produces this plot:
This image doesn't quite do it justice, so here is a plot of the full data set that shows how bad the overlap is:
Probably the simplest method, given your particular colour scheme, is to make the arrows all black but use the alpha channel to make the gray ones gray. This means the black arrows will still be completely black even if other arrows are drawn over the top:
fviz_pca_var(pca,
axes = c(1,2),
select.var = list(contrib = 200),
alpha.var = ifelse(markers_vec == 0, 0.2, 1),
label = "none") +
theme_bw()
I have a dataframe with observations of a temporal event. I can easily plot the observations based on their Xand Y values, either in base R or ggplot2:
plot(df$Y ~ df$X)
ggplot(df, aes(x=X, y=Y)) +
geom_point()
What I'd like to do though is add in the temporal dimension by coloring the observations differently depending on the time they were made. This could be done by defining thresholds for different time steps and assigning discrete colors to each time step. But what I'm interested in is a continuous coloring where the flow of time can be directly read off the flow of colors or color hues or color transparencies. Any help is greatly appreciated.
Reproducible data (slightly shorter than data underlying shown graphs):
structure(list(X = c(171.000358830368, 171.453956550099, 171.075127685269,
171.140924277581, 171.271520544141, 171.582558980529, 171.23762532992,
171.70218914837, 171.318375693213, 170.939546828382, 171.70218914837,
171.779948757467, 171.809856299427, 172.066064242221, 171.603494259901,
171.718139837416, 171.947430992445, 171.930483385334, 172.002261486039,
172.153793031971, 122.654817251517, 171.710164492893, 172.012230666692,
171.395138384244, 171.750041215507, 171.631407965731, 171.927492631138,
172.059085815763, 171.816834725885, 170.787018364384, 170.328436054326,
170.083194210252, 170.405198745358, 170.505887469957, 170.429124778926,
171.187779426653, 171.00634033876, 170.984408141322, 170.731190952725,
170.924593057402, 170.570687144205, 170.751129314032, 170.539782684179,
170.561714881617, 170.461026157017, 170.631499146191, 170.686329639785,
170.863781055416, 170.782033774058, 170.655425179759, 170.645455999106,
170.612557702949, 170.658415933955, 170.514859732545, 170.382269629855,
170.167932245806, 170.187870607113, 169.841940038438, 175.365863038509,
190.973612269532, 206.626222813496, 214.924568789419, 214.509850874236,
211.538038121446, 209.610995501136, 207.007045514458, 204.585531533737,
202.489012842319, 200.693563406635, 199.331773329375, 198.61997383072,
198.494362154486, 202.238786407918, 205.8695620019, 206.616253632843,
207.53541208909, 208.81146721273, 209.850255836818, 210.903998231887,
211.630751501523, 211.825150524265, 211.487195300113, 210.759445112412,
209.407624215806, 207.816542983517, 205.843642132201, 203.585622714198,
201.110275157949, 198.862224920598, 196.914247020916, 195.676074783759,
195.283289066014, 196.015026925976, 197.60710507633, 200.479226022586,
204.110001616568, 207.405812740595, 210.878078362188, 214.015379513825,
216.768867210305, 218.423751198776, 219.404718575074, 219.467524413191,
218.760709504862, 217.317172146244, 215.572565531892, 213.46408382369,
211.338654508377, 209.181323814973, 207.720838849244, 206.537497105681,
205.904454134187, 205.898472625795, 205.953303119389, 206.365030280376,
206.977137972497, 208.319986606516, 209.067675155523, 209.547192744954,
210.106463779612, 210.356690214013, 210.602928976153, 210.692651602034,
210.849167738293, 210.88405987058, 210.998705448095, 211.037585252643,
210.902004395756, 210.895025969299, 210.958828725481, 210.92692734739,
210.975776332592, 211.201079815359, 211.164193846942, 211.304759294155,
210.646793371028, 198.816366689593, 179.579835700717, 164.568243472769,
159.956500502489, 159.636489803513, 159.535801078913, 159.997374143168,
160.238628314981, 160.281495791791, 160.576583539133, 160.637395541118,
156.724492134644, 153.810500629644, 153.846389679996, 154.022844177562,
154.29300897327, 154.442546683072, 154.536256981214, 154.665856329709,
154.665856329709, 154.833338564686, 156.847113056681, 176.676810294436,
177.112463488991, 176.71868085318, 176.042770404877, 175.822451512436,
175.75964567432, 176.817375741649, 176.028813551962, 176.347827332872,
176.622976718907, 182.215687065486, 118.910392998086, 174.380907989949,
171.889609744655, 168.855988071813, 164.924143222097, 160.275514283399,
156.671655477181, 154.165403460906, 107.676123319726, 107.769833617868,
108.581324923058, 109.676937876871, 109.244275436511, 110.393721965853,
161.248506315174, 163.97707106002, 165.507340290323, 166.298893234206,
168.008607716271, 169.06434394747, 170.018394536004, 170.41118025375,
170.202824378093, 170.201827460028, 121.025853132745, 117.483803246579,
120.870333914552, 166.984772863163, 166.489304584687, 168.118268703459,
181.382263562859, 196.653054487796, 135.095357788943, 135.297732156208,
202.44614536551, 203.460011037964, 204.440978414263, 205.227546767819,
205.43490572541, 204.935449774673, 203.891676560258, 202.264706277617,
199.827241607851, 176.081650209426, 194.729002621683, 194.670184455827,
191.186952735516, 190.463190220076, 190.508051533017, 191.482040482858,
186.606114225261, 171.208714706025, 158.889798172571, 149.82880987666,
152.986049389604, 156.597883540345, 159.647455902232, 161.754940692369,
162.893421122992, 114.542894953814, 162.976165322415, 114.268742485845,
160.762010299286, 159.259654774813, 157.995562667957, 156.561994489993,
155.638848361484, 154.469463470836, 154.00390273432, 154.431580584353,
154.255126086787, 155.244068807608, 157.058459686534, 159.101144802424,
161.906472238301, 165.29399982434, 120.406767014167, 172.028181355737,
175.03588315888, 177.740521870158, 179.910812498411, 181.483949205524,
181.988389746588, 181.854802725832, 180.508963337618, 178.901931416283,
176.58609075049, 173.537515306668, 170.679351213327, 167.374567826712,
164.487493109477, 162.203553821774, 160.569605112675, 159.562717866678,
159.828894990125, 159.632502131252, 160.2495944137, 161.450880682439,
162.705003608642, 164.622077048298, 166.212161362521, 167.468278124855,
168.602770883216, 169.155063491416, 169.462114255542, 169.404293007752,
169.003531945484, 168.422328713389, 167.728473739909, 166.825265972708,
165.984864043623, 169.386348482576, 183.504702123976, 200.99862033463,
214.553715269111, 218.552353629205, 217.590327696148, 217.000152201465,
216.115885877505, 215.405083296914, 214.675339273083, 214.063231580962,
213.56477254829, 212.972603217475, 212.488101037718, 212.097309156103,
211.685581995116, 211.786270719716, 211.46426618461, 210.97178866033,
210.657759469747, 210.36167480434, 210.032691842776, 209.635918452769,
209.435537921635, 209.197274504018, 209.046739876151, 208.712772324261,
208.585166811897, 208.3319496233, 203.956476234505, 189.73444311431,
179.21097601654, 176.815381905519, 177.826256823777, 175.90718954799,
175.461567172782, 177.582011897768, 177.488301599626, 177.463378647992,
177.577027307441, 177.676719113976, 177.633851637166, 177.843204430888,
177.995732894886, 177.92196095805, 176.699739409939, 176.824354168107,
176.859246300394, 176.68877331122, 175.325986315895, 176.695751737677,
174.103764767783, 176.520294158177, 178.304777495142, 176.684785638959,
176.697745573808, 177.861148956064, 176.614004456319, 176.121526932039,
176.364774939983, 173.25239673998, 176.064702602315, 174.207444246579,
174.133672309744, 173.831606135945, 174.066878799366, 174.03597433934,
175.830426856959, 174.100774013587, 174.019026732229, 173.965193156701,
173.955223976047, 174.028995912883, 173.886436629538, 174.420784712563,
179.872929611928, 198.741597834692, 208.031877285631, 208.418681494985,
208.232257816765, 247.112062365176, 207.772678588642, 207.722832685375,
207.639091567886, 207.613171698187, 207.569307303312, 207.722832685375,
207.850438197739, 208.002966661736, 208.157488961865, 208.43463218403,
208.722741504914, 209.112536468464, 209.628940026312, 210.13038981318,
211.046557515231, 211.660659043483, 212.448224315105, 213.349438246175,
212.024534137333, 215.36321273817, 216.563502088844, 217.81264042472,
217.271313915238, 218.686937568026, 220.266055783531, 221.81426953901,
223.14814591044, 224.246749618449, 225.106092990776, 225.724182191289,
226.127934007753, 226.088057285139, 225.885682917874, 225.557696874376,
225.048271742986, 257.697338382997, 223.992535511786, 223.464168937154,
223.130201385264, 222.831125965661), Y = c(150.561649721083,
150.067708443465, 150.155166858818, 149.178216605502, 149.407794945803,
148.660423032789, 148.913853668186, 148.276798619991, 148.487493893341,
147.692417390134, 148.3284785927, 147.692417390134, 147.963737246853,
147.216365333839, 147.264069924031, 145.992941364529, 145.626212327424,
144.316323788391, 144.345145311632, 143.45465962804, 118.779460351005,
143.071035215242, 143.299619709914, 142.858352250634, 142.574112400738,
141.615051368744, 141.139993158078, 140.643070343573, 140.190870582374,
140.33597204421, 140.341935117984, 139.946384557638, 140.207765958067,
139.045960417756, 139.269575684283, 138.122677828407, 138.177339338002,
137.268964433088, 137.874216421154, 137.406115129891, 138.076960929472,
137.840425669768, 138.343311558047, 137.691348825417, 138.413874597706,
137.837444132881, 138.70109598449, 138.061059399408, 138.395985376384,
137.616810403241, 138.267779290242, 137.52736429663, 138.054102480005,
137.291822882555, 137.941797923927, 137.242130601105, 137.811604146527,
137.343502855264, 137.326607479571, 135.242513195539, 137.238155218589,
138.312502343547, 139.236778778526, 138.86110513076, 139.273551066799,
138.584816045896, 138.903840492808, 137.509475075308, 137.144733729462,
135.697694493625, 135.567500716225, 134.488184363121, 135.142134787009,
134.959267191271, 135.420411563131, 134.819134957581, 135.649989903432,
135.13120248509, 136.041565081262, 135.646014520916, 136.444072561011,
135.969014350344, 136.521592520073, 135.799066747784, 136.416244883398,
135.748380620704, 136.195611153758, 135.349848523472, 135.534703810467,
134.372898270156, 134.413645940945, 133.470486439016, 133.743793986993,
133.129597388266, 133.61658174648, 133.007354375898, 133.341286507245,
132.799640639435, 133.348243426648, 132.845357538369, 133.241901944344,
132.396139314057, 132.584969983569, 131.617958186543, 131.839585761812,
130.896426259883, 130.981896983977, 129.944322147292, 129.993020583113,
128.968365739605, 128.960414974573, 127.798609434262, 128.082849284158,
127.096954420181, 127.450763464109, 126.738176148109, 127.149628238519,
126.326724057699, 127.023409843635, 126.595062377532, 127.240068190759,
126.893216066235, 127.574000322106, 127.072108279456, 127.73997254215,
127.304668156644, 128.098750814222, 127.663446428716, 128.521135206551,
128.014273935757, 128.852085801011, 128.420756798021, 129.41460242703,
128.987248806556, 129.867796033858, 129.49510392298, 129.345033233,
128.14148617627, 130.371675767766, 130.815924763933, 131.900210345182,
131.529505925561, 132.387194703396, 131.814739621087, 132.715163760969,
132.186437886336, 133.071954341783, 132.563105379731, 133.595710988271,
133.373089567373, 134.173135298725, 133.956476951601, 134.797270353743,
134.562722785297, 135.403516187438, 135.163005545218, 135.992866645441,
131.572241287609, 129.463300862852, 122.968519677278, 123.474387102444,
123.527060920781, 124.246605156184, 123.732786965986, 117.399008772312,
123.862980743386, 124.421521986889, 123.70893467089, 124.808127936574,
126.778923818898, 124.377792779213, 122.935722771521, 124.517925012903,
124.076657553623, 124.838937151073, 124.239648236781, 124.81210331909,
125.153986215469, 125.120195464083, 123.901740722918, 123.046039636341,
120.966914580454, 120.263271875116, 121.53241274336, 121.521480441441,
120.071459668717, 121.39228050967, 120.229481123729, 120.093324272555,
118.693989626911, 118.729768069555, 117.589827133082, 118.195079121148,
108.618382640018, 113.465367772695, 107.992259893742, 117.911833116881,
117.31055651133, 118.451491293432, 122.158535489636, 128.127572337464,
125.310019979223, 125.468041434236, 128.644372064548, 129.146264107198,
128.288575329363, 128.90376577372, 128.18124000143, 128.778541224465,
128.358144523394, 129.135331805279, 128.537036736615, 128.94351959888,
128.079867747271, 127.956630889274, 126.952846803975, 126.741157684996,
125.517733715686, 125.563450614621, 124.14026367388, 123.462460954896,
122.454701487081, 125.26330923466, 123.721854664067, 123.747694650421,
122.82639975233, 123.206048782611, 122.267858508827, 115.131053046914,
122.130707812024, 114.510893374412, 122.343390776632, 123.249777990288,
122.697199820559, 123.432645586025, 123.527060920781, 123.925593018014,
123.003304274293, 124.011063742109, 123.917642252982, 125.041681659391,
124.907512499475, 125.786072035519, 125.602210594152, 126.27603793062,
114.507911837525, 126.00571191953, 125.382570710141, 125.88446275279,
125.144047759179, 125.545561393298, 124.766386420155, 124.993977069198,
124.193931337846, 124.197906720362, 123.290525661077, 123.403824062784,
122.544147593691, 122.783664390282, 121.98361865893, 122.607753713948,
122.271833891343, 122.822424369814, 122.699187511817, 122.789627464057,
122.228104683666, 122.501412231644, 121.437997408604, 121.568191186005,
120.996729949324, 121.368428214574, 120.744293159556, 121.303828248688,
120.7800716022, 121.339606691332, 120.856597715634, 121.619871158713,
121.218357524593, 121.974674048269, 121.564215803488, 122.684279827382,
121.975667893898, 124.308223585182, 126.514560881582, 128.363113751539,
128.091793894819, 128.702015111031, 128.509209059003, 129.155208717859,
128.840159653463, 129.56566696264, 129.235710213809, 130.013891341323,
129.748534558377, 130.711570972887, 130.380620378427, 129.980100589936,
129.586537720849, 130.563487974165, 130.419380357958, 131.614976649656,
131.552364375028, 132.790696028774, 132.718145297856, 133.984304629213,
134.010144615568, 135.165987082105, 135.426374636905, 136.790924685535,
135.130208639461, 129.957242140469, 126.592080840645, 127.358335820611,
127.887061695243, 129.049861081184, 128.986254960927, 130.73840480487,
130.834807830884, 132.108917927273, 132.045311807017, 133.056052811719,
132.901012893594, 133.848147778039, 133.514215646692, 134.505079738814,
134.604464301715, 135.360780825391, 135.037780995963, 135.768257533284,
134.735651924744, 135.809999049703, 134.344076746915, 135.788134445865,
127.311625076047, 135.801054439042, 135.322020845859, 127.889049386501,
135.320033154601, 135.615205306417, 126.665625417191, 134.175122989983,
126.25616101804, 134.430541316638, 133.848147778039, 134.240716801498,
133.809387798508, 134.271526015997, 126.20547489096, 134.244692184014,
133.780566275267, 134.294384465464, 133.818332409169, 134.256618331562,
133.891876985716, 134.534895107684, 127.426911169012, 142.477709374724,
146.033689035318, 147.092134630213, 146.619064110804, 170.36998695286,
146.552476453661, 147.216365333839, 146.709504063044, 147.474765197381,
147.046417731278, 147.842488080114, 147.290903756014, 148.023367984594,
147.353516030642, 147.901124972226, 147.171642280533, 147.729189678407,
147.072257717632, 147.535389780751, 147.15474690484, 146.832740921041,
146.568377983725, 146.13207975259, 144.739702026348, 145.216747928273,
144.640317463447, 144.137431575169, 143.147561328676, 141.941032735059,
141.073405500934, 140.496975036109, 140.057695268087, 139.778424646336,
139.654193942709, 139.203981872768, 138.923717405388, 138.393997685126,
137.79868415335, 137.209333695347, 136.706447807069, 165.401752653445,
135.932242062071, 135.629119145223, 135.168968618992, 134.679002723891
), starttime_ms = c(310849, 310865, 310882, 310899, 310915, 310932,
310949, 310965, 310982, 310999, 311015, 311032, 311049, 311065,
311082, 311099, 311115, 311132, 311149, 311165, 311182, 311199,
311215, 311232, 311249, 311265, 311282, 311299, 311315, 311332,
311349, 311365, 311382, 311399, 311415, 311432, 311449, 311465,
311482, 311499, 311515, 311532, 311549, 311565, 311582, 311599,
311615, 311632, 311649, 311665, 311682, 311699, 311715, 311732,
311749, 311765, 311782, 311799, 311815, 311832, 311849, 311865,
311882, 311899, 311915, 311932, 311949, 311965, 311982, 311999,
312015, 312032, 312049, 312065, 312082, 312099, 312115, 312132,
312149, 312165, 312182, 312199, 312215, 312232, 312249, 312265,
312282, 312299, 312315, 312332, 312349, 312365, 312382, 312399,
312415, 312432, 312449, 312465, 312482, 312499, 312515, 312532,
312549, 312565, 312582, 312599, 312615, 312632, 312649, 312665,
312682, 312699, 312715, 312732, 312749, 312765, 312782, 312799,
312815, 312832, 312849, 312865, 312882, 312899, 312915, 312932,
312949, 312965, 312982, 312999, 313015, 313032, 313049, 313065,
313082, 313099, 313115, 313132, 313149, 313165, 313182, 313199,
313215, 313232, 313249, 313265, 313282, 313299, 313315, 313332,
313349, 313365, 313382, 313399, 313415, 313432, 313449, 313465,
313603, 313619, 313636, 313665, 313682, 313699, 313715, 313732,
313749, 313765, 313782, 313799, 313815, 313832, 313849, 313865,
313882, 313899, 313915, 313932, 313949, 313965, 313982, 313999,
314015, 314032, 314049, 314065, 314082, 314099, 314115, 314132,
314149, 314165, 314182, 314199, 314215, 314232, 314249, 314265,
314282, 314299, 314315, 314332, 314349, 314365, 314382, 314399,
314415, 314432, 314449, 314465, 314482, 314499, 314515, 314532,
314549, 314565, 314582, 314599, 314615, 314632, 314649, 314665,
314682, 314699, 314715, 314732, 314749, 314765, 314782, 314799,
314815, 314832, 314849, 314865, 314882, 314899, 314915, 314932,
314949, 314965, 314982, 314999, 315015, 315032, 315049, 315065,
315082, 315099, 315115, 315132, 315149, 315165, 315182, 315199,
315215, 315232, 315249, 315265, 315282, 315299, 315315, 315332,
315349, 315365, 315382, 315399, 315415, 315432, 315449, 315465,
315482, 315499, 315515, 315532, 315549, 315565, 315582, 315599,
315615, 315632, 315649, 315665, 315682, 315699, 315715, 315732,
315749, 315765, 315782, 315799, 315815, 315832, 315849, 315865,
315882, 315899, 315915, 315932, 315949, 315965, 315982, 315999,
316015, 316032, 316049, 316065, 316082, 316099, 316115, 316132,
316149, 316165, 316182, 316199, 316215, 316232, 316249, 316265,
316282, 316299, 316315, 316332, 316349, 316365, 316382, 316399,
316415, 316432, 316449, 316465, 316482, 316499, 316515, 316532,
316549, 316565, 316582, 316599, 316615, 316632, 316649, 316665,
316682, 316699, 316715, 316732, 316749, 316765, 316782, 316799,
316815, 316832, 316849, 316865, 316882, 316899, 316915, 316932,
316949, 316965, 316982, 316999, 317015, 317032, 317049, 317065,
317082, 317099, 317115, 317132, 317149, 317165, 317182, 317199,
317215, 317236, 317253, 317269, 317286, 317303, 317319, 317336,
317353, 317369, 317386, 317403, 317419, 317436, 317453, 317469,
317486, 317503, 317519, 317536, 317553, 317569, 317586, 317603,
317619, 317636, 317653)), row.names = c(123318L, 123324L, 123333L,
123343L, 123354L, 123363L, 123370L, 123373L, 123376L, 123384L,
123391L, 123399L, 123406L, 123414L, 123424L, 123434L, 123438L,
123440L, 123444L, 123446L, 123449L, 123459L, 123468L, 123477L,
123483L, 123491L, 123502L, 123505L, 123509L, 123516L, 123523L,
123526L, 123532L, 123538L, 123548L, 123556L, 123565L, 123570L,
123574L, 123576L, 123580L, 123586L, 123593L, 123602L, 123611L,
123617L, 123623L, 123627L, 123631L, 123640L, 123648L, 123652L,
123663L, 123671L, 123674L, 123677L, 123684L, 123691L, 123696L,
123706L, 123714L, 123725L, 123733L, 123735L, 123739L, 123741L,
123745L, 123751L, 123756L, 123766L, 123775L, 123777L, 123786L,
123794L, 123800L, 123807L, 123813L, 123816L, 123819L, 123823L,
123833L, 123842L, 123850L, 123861L, 123869L, 123876L, 123884L,
123892L, 123897L, 123904L, 123910L, 123912L, 123917L, 123926L,
123933L, 123939L, 123945L, 123949L, 123953L, 123956L, 123962L,
123968L, 123975L, 123981L, 123988L, 123996L, 124005L, 124015L,
124020L, 124023L, 124029L, 124036L, 124045L, 124054L, 124063L,
124069L, 124079L, 124085L, 124094L, 124103L, 124112L, 124117L,
124124L, 124127L, 124130L, 124137L, 124146L, 124158L, 124169L,
124177L, 124184L, 124189L, 124195L, 124203L, 124213L, 124220L,
124227L, 124234L, 124244L, 124252L, 124257L, 124259L, 124262L,
124265L, 124271L, 124279L, 124291L, 124299L, 124304L, 124312L,
124321L, 124324L, 124328L, 124337L, 124343L, 124347L, 124354L,
124363L, 124415L, 124424L, 124426L, 124435L, 124440L, 124446L,
124454L, 124461L, 124467L, 124469L, 124476L, 124486L, 124493L,
124496L, 124504L, 124510L, 124519L, 124524L, 124531L, 124539L,
124550L, 124555L, 124558L, 124560L, 124564L, 124569L, 124579L,
124586L, 124595L, 124599L, 124608L, 124617L, 124622L, 124629L,
124634L, 124636L, 124640L, 124645L, 124653L, 124660L, 124669L,
124681L, 124694L, 124700L, 124709L, 124716L, 124722L, 124727L,
124730L, 124733L, 124738L, 124743L, 124748L, 124755L, 124766L,
124774L, 124782L, 124786L, 124793L, 124801L, 124808L, 124812L,
124821L, 124831L, 124844L, 124852L, 124861L, 124867L, 124873L,
124877L, 124883L, 124890L, 124901L, 124912L, 124917L, 124921L,
124926L, 124934L, 124943L, 124948L, 124956L, 124958L, 124962L,
124968L, 124977L, 124986L, 124995L, 125005L, 125013L, 125015L,
125018L, 125025L, 125036L, 125046L, 125051L, 125056L, 125063L,
125072L, 125076L, 125078L, 125082L, 125084L, 125091L, 125100L,
125111L, 125118L, 125128L, 125135L, 125143L, 125147L, 125155L,
125165L, 125168L, 125170L, 125178L, 125188L, 125195L, 125203L,
125212L, 125217L, 125220L, 125222L, 125226L, 125238L, 125245L,
125252L, 125256L, 125262L, 125269L, 125278L, 125285L, 125288L,
125292L, 125296L, 125307L, 125319L, 125322L, 125325L, 125329L,
125334L, 125342L, 125351L, 125360L, 125372L, 125379L, 125386L,
125389L, 125392L, 125397L, 125402L, 125411L, 125421L, 125428L,
125436L, 125446L, 125451L, 125457L, 125459L, 125465L, 125470L,
125477L, 125483L, 125490L, 125499L, 125508L, 125515L, 125523L,
125526L, 125532L, 125540L, 125548L, 125550L, 125553L, 125558L,
125564L, 125570L, 125580L, 125588L, 125594L, 125600L, 125604L,
125610L, 125617L, 125624L, 125633L, 125640L, 125647L, 125655L,
125661L, 125669L, 125676L, 125679L, 125686L, 125695L, 125706L,
125716L, 125725L, 125733L, 125737L, 125743L, 125755L, 125758L,
125762L, 125764L, 125769L, 125775L, 125783L, 125790L, 125800L,
125808L, 125815L, 125823L, 125828L, 125834L, 125838L, 125845L,
125852L, 125864L, 125871L, 125880L, 125889L, 125893L, 125896L,
125899L, 125902L, 125904L, 125912L, 125920L, 125931L, 125936L,
125944L, 125952L, 125961L, 125963L, 125969L, 125978L, 125984L,
125991L, 126002L, 126009L, 126018L, 126025L, 126028L), class = "data.frame")
You could do:
library(ggplot2)
ggplot(df, aes(x = X, y = Y, color = starttime_ms)) +
geom_point() +
scale_color_gradientn(colors = c("red", "gold", "forestgreen"), name = "time")
An alternative that might help keep track of the motion better would be to use geom_segment:
df$X2 <- dplyr::lag(df$X)
df$Y2 <- dplyr::lag(df$Y)
ggplot(df, aes(x = X, y = Y, color = starttime_ms)) +
geom_segment(aes(xend = X2, yend = Y2), size = 1) +
scale_color_gradientn(colors = c("red", "gold", "forestgreen"), name = "time")
I am trying to create a stratigraphic plot of geochemical element data which should be possible using package tidypaleo.
I want multiple plots of the different element data with Depth (cm) downcore set as the y axis. The data look as follows.
Image of data
I am using this code:
ggplot(wapITRAX, aes(x =BrTi , y = wapITRAX$Depth))+
labs(y = "Depth (cm)")+
geom_lineh()+
theme_classic()+
scale_y_reverse()
However, this only plots one element and I am trying to achieve a plot like this Image of plot
> dput(head(wapITRAX))
structure(list(Depth = 0:5, IncCoh = c(6.049230907, 5.975282432,
5.736199822, 5.658584418, 5.659008377, 5.597103404), BrTi =
c(50.50197628,
22.09236453, 23.48370927, 18.62638581, 14.36924414, 17.48777896
), AlIncCOh = c(16.69633736, 8.200449193, 23.70907643, 20.32310407,
28.62692352, 26.44224866), BrCl = c(8.04090623, 4.306048968,
3.417836951, 3.156895904, 2.787628518, 2.059316731), FeTi =
c(332.715415,
235.9371921, 372.726817, 390.7871397, 396.986099, 495.2624867
), CaTi = c(4.071146245, 3.27955665, 4.395989975, 3.677383592,
3.028670721, 4.523910733), ZrRb = structure(c(363L, 447L, 407L,
395L, 450L, 410L), .Label = c("#DIV/0!", "0.447638604",
"0.478169284",
"0.54554134", "0.548501778", "0.561420163", "0.579454254",
"0.579498861",
"0.580801291", "0.589758019", "0.590194076", "0.590277778",
"0.591357754",
"0.592870544", "0.593851133", "0.598519653", "0.599931082",
"0.600979737",
"0.601426307", "0.611710677", "0.617065868", "0.618499405",
"0.621310093",
"0.627720871", "0.63775246", "0.64005168", "0.643958869",
"0.644371941",
"0.645605974", "0.645661658", "0.646672915", "0.647348952",
"0.651578947",
"0.652401176", "0.656186383", "0.657906264", "0.658835905",
"0.662074554",
"0.662361624", "0.669589393", "0.67103429", "0.671371769",
"0.674335863",
"0.674781688", "0.676097561", "0.676639083", "0.677849462",
"0.680497925",
"0.680610514", "0.680725971", "0.683906537", "0.68855859",
"0.689067202",
"0.692353115", "0.692732291", "0.695411392", "0.696067091",
"0.696794872",
"0.699376436", "0.701762744", "0.702015197", "0.702432938",
"0.70361991",
"0.705235754", "0.705426357", "0.708084164", "0.708258528",
"0.708925221",
"0.715226656", "0.715263314", "0.717828827", "0.718975706",
"0.719799305",
"0.720363636", "0.72476489", "0.725426857", "0.725461098",
"0.726030739",
"0.7267645", "0.726998188", "0.727170554", "0.727533265",
"0.730362368",
"0.731182796", "0.735042735", "0.735849057", "0.736184046",
"0.737214792",
"0.738692342", "0.742223591", "0.742639327", "0.742714724",
"0.744126167",
"0.745790081", "0.746207701", "0.746606335", "0.747041077",
"0.749019608",
"0.751204307", "0.751326495", "0.752443737", "0.752994012",
"0.753398969",
"0.754587869", "0.755151515", "0.755466053", "0.756316411",
"0.757107679",
"0.759175941", "0.76070965", "0.761635833", "0.763760166",
"0.763861189",
"0.764076577", "0.765001936", "0.765509391", "0.765795207",
"0.76598579",
"0.767080745", "0.767225748", "0.767262192", "0.770291777",
"0.771194699",
"0.772540984", "0.774533358", "0.77486376", "0.777252364",
"0.778319123",
"0.778445883", "0.779076739", "0.779310345", "0.782186577",
"0.78238342",
"0.782541093", "0.785497178", "0.787594824", "0.788230584",
"0.788916736",
"0.789033068", "0.789052737", "0.789177605", "0.789596048",
"0.790439133",
"0.791078234", "0.792993631", "0.794899978", "0.795265235",
"0.795432921",
"0.798130469", "0.79861957", "0.798894446", "0.799410029",
"0.799679487",
"0.800096223", "0.800119868", "0.800280899", "0.800947004",
"0.80293448",
"0.804118993", "0.804120112", "0.805092452", "0.806161301",
"0.81631016",
"0.816337149", "0.81638756", "0.81827622", "0.818899466",
"0.81993865",
"0.821925431", "0.825329202", "0.825608077", "0.826135414",
"0.826325411",
"0.827108292", "0.829798658", "0.832460733", "0.832814584",
"0.833845295",
"0.833908046", "0.833980583", "0.83741705", "0.837604457",
"0.838292367",
"0.838860231", "0.839741935", "0.846485664", "0.846740374",
"0.84679304
1",
"0.84734599", "0.850027518", "0.850077882", "0.854863613",
"0.85565883",
"0.85591192", "0.857095047", "0.858967536", "0.85897779",
"0.859192457",
"0.859239235", "0.860408412", "0.860505166", "0.8609918",
"0.862821134",
"0.867779204", "0.872207328", "0.872918493", "0.873888706",
"0.874140666",
"0.877452229", "0.879856851", "0.88003663", "0.880155093",
"0.881013197",
"0.882190156", "0.882319978", "0.883029342", "0.883316008",
"0.884908053",
"0.885431692", "0.885849846", "0.889082969", "0.89763062",
"0.898384089",
"0.898527865", "0.9", "0.900719424", "0.904636261", "0.904913848",
"0.905223479", "0.908485273", "0.911212059", "0.916470588",
"0.917162698",
"0.917470525", "0.917952884", "0.918054256", "0.918781726",
"0.922268448",
"0.925987182", "0.926287744", "0.933958724", "0.939701616",
"0.940987438",
"0.941196465", "0.943041526", "0.946832078", "0.947419907",
"0.948080043",
"0.949085366", "0.950099404", "0.950502912", "0.950704225",
"0.953596288",
"0.954868709", "0.95505992", "0.955361596", "0.956422018",
"0.95671509",
"0.959064327", "0.96165309", "0.964232489", "0.965528048",
"0.965874467",
"0.968717195", "0.970270821", "0.971573209", "0.973415133",
"0.974608081",
"0.974986972", "0.978223496", "0.97935022", "0.980142566",
"0.981730052",
"0.982819606", "0.987483236", "0.989587207", "0.991002571",
"0.991932655",
"0.995306859", "0.997555012", "1.001706485", "1.002357873",
"1.005146199",
"1.010493827", "1.011544012", "1.013038906", "1.013469577",
"1.015357613",
"1.01541976", "1.015544041", "1.015937059", "1.016162826",
"1.019393939",
"1.019529957", "1.019994873", "1.028932491", "1.029518072",
"1.030101225",
"1.032310705", "1.03344968", "1.036769138", "1.037322515",
"1.041150092",
"1.041459782", "1.041961577", "1.042347697", "1.047411444",
"1.048390581",
"1.057811121", "1.060218978", "1.065395654", "1.06921167",
"1.071509648",
"1.073730469", "1.076014602", "1.077266637", "1.079983072",
"1.083003953",
"1.086196504", "1.086445657", "1.0867266", "1.0905", "1.092152628",
"1.094730942", "1.095692666", "1.1074142", "1.109534807",
"1.111154446",
"1.111188325", "1.113985094", "1.114045618", "1.128810226",
"1.129375525",
"1.132705479", "1.154833837", "1.155201819", "1.158940397",
"1.161348047",
"1.162953533", "1.175788329", "1.178036097", "1.179280397",
"1.18116463",
"1.198506534", "1.203134068", "1.207592892", "1.208610568",
"1.212525667",
"1.212841855", "1.214403157", "1.216756112", "1.228462377",
"1.23305986",
"1.241758242", "1.258212878", "1.265571914", "1.265707797",
"1.272666444",
"1.279325988", "1.289559543", "1.291191103", "1.296541003",
"1.305111821",
"1.315201411", "1.342263532", "1.344827586", "1.347944377",
"1.352831595",
"1.377717391", "1.37771934", "1.385440181", "1.403358682",
"1.411889597",
"1.412269549", "1.412544902", "1.433843384", "1.444088526",
"1.473706353",
"1.50172117", "1.508521601", "1.53030303", "1.532658694",
"1.537642783",
"1.54340949", "1.558746736", "1.595342983", "1.596774194",
"1.596969697",
"1.630653266", "1.702558635", "1.70468948", "1.718568102",
"1.740819711",
"1.760157274", "1.775510204", "1.804859813", "1.814081408",
"1.816513761",
"1.849050827", "1.870188679", "1.880195599", "1.887061404",
"1.91858679",
"1.929152149", "1.944140197", "1.946996466", "1.986547085",
"2.007683864",
"2.070983811", "2.099778271", "2.156359393", "2.159613059",
"2.163963964", "2.25951087", "2.261603376", "2.329896907",
"2.402555911", "2.414500684",
"2.444075305", "2.450268817", "2.484520124", "2.51119403",
"2.515064562",
"2.526086957", "2.554", "2.609715243", "2.61965812", "2.643854749",
"2.704166667", "2.883275261", "3.013186813", "3.02739726",
"3.206896552",
"3.320930233", "3.411627907", "3.688931298", "3.709677419",
"3.748267898",
"3.878865979", "3.936440678", "3.994230769", "33.15909091",
"4.095854922",
"4.29330254", "4.390957447", "4.514634146", "4.6367713",
"4.847665848",
"5.284023669", "5.387755102", "6.171339564", "6.183908046",
"6.36121673",
"6.847826087", "7.003496503", "7.193220339", "8.160550459",
"8.751879699"
), class = "factor"), MnFe = c(0.012176723, 0.010329834,
0.009460859,
0.004488071, 0.0033725, 0.003435313), MnIncCoh = c(169.4430276,
331.1977339, 490.5686845, 279.5752228, 272.3091921, 286.0408118
), CuRb = c(0.392971246, 1.484304933, 0.735426009, 0.491651206,
1.142857143, 0.4345898)), row.names = c(NA, 6L), class =
"data.frame")
Using your posted data. This should approximate the desired design.
First step, Transform the data from a wide format to a long format using the pivot_longer function from tidyr.
Then plot the data using "depth" as the independent variable and the parameters' values as the dependent variables.
Then use facet_wrap() to separate the plots. coord_flip() will make the independent variable (Depth) appear on the y-axis.
#fixed 1 column of data.
originaldata$ZrRb <- as.numeric(as.character(originaldata$ZrRb))
library(tidyr)
#Make wide
wapITRAX<-pivot_longer(originaldata, -1, names_to="parameter", values_to = "value")
library(ggplot2)
ggplot(wapITRAX, aes(x = Depth , y = value))+
labs(x = "Depth (cm)")+
geom_line() +
theme_classic() +
coord_flip() +
scale_x_reverse() +
facet_wrap(vars(parameter), nrow=1, scales = "free_x")
It is my data:
> dput(data)
structure(list(foldchange = c(-0.17853057272962, 3.60013440830337,
0.648944710423407, 1.38528656859267, 2.38882890772698, 1.91371568283765,
1.77591931363495, -1.51447851175922, 3.1416903855924, 1.51711016957237,
3.14707703341916, -1.44751697381751, 1.23658565660726, -0.512829478520189,
1.68928069854351, 2.07214007434345, 1.24799276690488, 6.25149659558487,
6.35918877435554, 5.86088034655694, 6.38890659730165, 5.05510489389194,
4.62060389613534, 3.75508710774868, 4.18575763169519, 5.31627264153051,
5.87091236649665, 6.71464565321037, 5.24000610137973, 4.25821377851955,
7.32277714374523, 3.1963295806222, 7.26249808789293, 4.44427454088613,
6.21495395454133, 8.74469985969472, 7.49982946564144, 4.45020943795387,
5.66199031471621, 5.29959827685333, 8.65819317196484, 5.86664903755707,
4.5740575604176, 8.24504501687473, 5.7916074097308, 4.18199181353134,
6.73956641707995, 4.60357435173805, 5.9205153184753, 3.65014593638562,
3.25607795403669, 5.56919529940933, 5.76811109641351, 6.10600807588152,
5.69234974521511, 5.5102283323841, 4.71232921328194, 6.55727667796477,
6.19995053763513, 4.64209842048131, 2.29238227264409, 6.79465189260383,
7.51968952300944, 7.81695579226993, 6.29926703626301, 4.64687557749141,
-2.44220257171186, 5.33199370895397, 5.18820654974805, 5.03498241997507,
6.29395095024283, 6.27602377186869, 6.78363927671209, 2.93759015053983,
6.65061604346668, 5.671080311536, 8.45199131823131, 5.79230415012306,
6.3270025568739, 7.0934690916107, 3.53800869528685, 3.08683779646569,
6.82111375813946, 3.02729078403818, 5.36024214796805, 4.04778690916444,
5.74765756930797, 4.10788604670319, 6.39978058654016, 5.7746717387066,
2.9247167920294, 7.54315906042106, 4.2742172444481, 6.61121261965006,
3.77012175922873, 4.94407566887151, 7.93185716981795, 7.05304621480995,
5.59261760605766, 5.42381827536197, 4.22645498896606, 5.15806113482742,
4.15403623593809, 6.40153592433128, 7.38902001442131, 2.72654942454391,
3.28741231093207, 3.79334363176751, 5.86527050546341, 5.10320299162235,
3.99883612485974, 6.45475273104195, 4.85567883821983, 6.55055641729645,
5.03746875764267, 7.27660375171087, 3.30817125205364, 5.23766518187252,
7.6588755830143, 3.53552741086444, 5.66455197986778, 8.40623211540503,
3.93151438658523, 4.26875667827774, 4.38704995079332, 6.75232207417316,
3.76563594385214, 6.08008097541859, 2.40905905886796, 5.50981339395085,
5.78780269825563, 0.2322416329745, 5.69410860233132, 4.94656296117567,
4.20594226169741, 4.50293112094816, 6.07430576125864, 3.67684848946483,
5.825851099141, 5.22439201628482, 7.72829018644622, 5.24910611944979,
4.01783420322782, 6.3888069709767, 5.26066649741256, 4.81678726754752,
6.5683773907454, 4.86957242886115, 6.76705114368644, 4.45769029291236,
7.77607596853254, 6.85213457577069, 1.40150885676552, 5.43409652313493,
3.21738153172066, 6.23015085020594, 5.50091556711613, 3.99948543388746,
5.85816098688073, 4.33775608630599, 5.91715214825299, 5.45674826103132,
7.66790792082782, 6.63325838131012, 3.89631178894691, 2.38526575667126,
1.58661549426288, 4.76626341270591, 6.73426272316295, 5.54006035262931,
4.07996836453406, 7.12087390022358, 6.96007461543701, 5.68202490906633,
6.58504044389069, 5.41036820315057, 6.61076809589319, 1.23772469006557,
4.1661166499875, 2.94059625298825, 8.38336956160413, 4.84906289871508,
4.93787691221829, 6.82379835301371, 5.82520798412864, 4.87582657907206,
5.36621724700676, 8.91922991774938, 3.49025109999629, 5.1232073414505,
4.27193651596412, 5.07417945071012, 3.61930149745523, 2.7469092502892,
6.67162003616042, 4.86698118654996, 7.53876919017093, 4.58878989189686,
5.78956520376246, 3.98567045767003, 8.14934433289609, 6.88879716040936,
7.00251456012974, 5.05095662412332, 6.39777439550296, 7.96268799093557,
4.82826575143863, 8.31032763539508, 4.74493707321909, 6.8929416113222,
6.84202549278968, 8.20121430968127, 2.91031632522241, 4.86190488550545,
5.5516465446887, 7.74996457744065, 2.25505738807845, 5.71069872298306,
1.97493599527532, 5.60445326341706, 6.39297603198736, 7.16298115056911,
4.52688105225386, 6.46061751569601, 4.78104064111529, 2.84526825975018,
6.5537923066898, 6.98258253798747, 5.0967396817644, 6.64593966293456,
4.8990397150507, 4.59878411928317, 5.55158425631398, 2.1065660739172,
2.40396884881286, 6.45421536580342, 5.98567305090568, 6.48593538806214,
7.41313242816247, 5.99024340460149, 5.63812101302136, -1.43793573368627
), all_pvalue = c(0.818887590433193, 0.00892139546812015, 0.434133425685163,
0.0536266013313456, 0.0450933986128537, 0.0367856407800243, 0.0941222253709068,
0.213526299326008, 0.00855327289085924, 0.0449444491492238, 0.00465098209958804,
0.0369667514121697, 0.0910501610463896, 0.378892060498093, 0.0595757077704777,
0.031626850730261, 0.0878934608628569, 0.0124722939899662, 0.0249040599008334,
0.0150448394759397, 0.0104024068916351, 0.0340577599419123, 0.0244947271472485,
0.122485247246688, 0.0211309039009709, 0.0201043058824927, 0.0152779076456381,
0.00938733157248341, 0.013807428170544, 0.01948348499862, 0.0119978631408916,
0.0997968486684091, 0.00504808432168479, 0.041226720485986, 0.0127407583225205,
0.00709338243276709, 0.00383605674079435, 0.0399854589187244,
0.0163222001450531, 0.00909800553027099, 4.81144191594885e-06,
0.0105231068132293, 0.0377428014886314, 0.00709059633291303,
0.00743882656849872, 0.0696633906261403, 0.00568508439144595,
0.0142206230933159, 0.0183625193075117, 0.0818974933908099, 0.0609423408539195,
0.00581572852382799, 0.00603085345605447, 0.00684099077236254,
0.0194814381299995, 0.0325641567545152, 0.0404062983698557, 0.00626431765905907,
0.00287654691487974, 0.0183318967557602, 0.202860908663261, 0.0033395497287839,
0.00804091896430431, 0.00524635934550195, 0.0100089274728679,
0.00470611875383887, 0.361169323059008, 3.90129727113067e-06,
0.0205225005371219, 0.0120500045076898, 0.00732055098038156,
0.0229916087025324, 0.00544774481324614, 0.108252753362848, 0.00911923198666818,
0.0122812312739145, 0.00343585528287351, 2.06257208918569e-05,
0.000904993210532763, 0.0015065294739414, 0.102118204143709,
0.118350948568527, 0.0136202759386966, 0.15251012082679, 0.0428316882385798,
0.0752744217284719, 0.00632986043900174, 0.0269567937932686,
0.00707537967267082, 0.0149356605279715, 0.163005190656644, 0.00929612911973378,
0.0196453775259569, 0.0133262667903121, 0.114115405959882, 0.0189109801950218,
0.00590387539250432, 0.00802692325541374, 0.0184346327727756,
0.00104714399950925, 0.0554774130259537, 0.0084910975380844,
0.0349856904843115, 0.0124547169142572, 8.60050865459788e-07,
0.1188515828269, 0.133720711339729, 0.0641704698591151, 0.0075124796175742,
3.94432669779951e-07, 0.0740343932996142, 0.00160775849728933,
6.61326355967731e-06, 0.0106610228625055, 0.00539717052083514,
0.00492255859958016, 0.122531121480312, 0.0159768620962635, 0.00365780586610517,
0.0135086464724098, 0.0256265865459836, 0.00391215568816396,
0.0271580638871089, 0.0512876942387616, 0.0135566028247977, 0.0169110062500104,
0.0848247460082605, 0.0158705161056627, 0.176345767878009, 0.0124719098589431,
0.0152388258990332, 0.866066700538701, 0.042979313425954, 0.0160125031962862,
0.0441592105265668, 0.0124108545467876, 1.80875437447348e-06,
0.0704671677844812, 0.0111351361909711, 1.62883074487165e-07,
0.000694779973318456, 0.0120079549431507, 0.0196183531679123,
0.0100186493225724, 0.0173294242221405, 0.0573614373022037, 1.30612522568528e-07,
0.023898721968545, 0.00665918028588502, 0.06907634259105, 0.00518730999717143,
0.0109717740506543, 0.394662670743417, 0.0165847750353483, 0.0593829446004973,
0.000113216713641592, 0.0222583368635018, 0.0694462386106761,
0.00437955933335859, 0.0631677226779205, 0.00649674049335009,
1.10699021652115e-07, 0.00543997929535101, 0.00292890932795308,
0.021574456097881, 0.17997075681454, 0.279770535484078, 0.050945514039484,
0.0102923410906512, 0.0103389721465925, 0.112227938286441, 0.00398090342551613,
0.0428676019413789, 0.00337313863923396, 0.0092822848687081,
0.00778850900332348, 0.00146877357609246, 0.448950342618974,
0.0801294647165026, 0.0797698230881222, 3.72112644308374e-05,
0.00728653989704988, 0.0320006421510141, 0.000420503611946363,
0.0207086037412542, 0.0448889114898146, 6.94555482668648e-07,
0.000135994029220134, 0.0982385638169219, 0.0130399195487442,
0.0297056523919866, 0.0208455457844783, 0.104595177359326, 0.100684824982166,
0.00110472016462074, 0.0664329210478157, 0.0140428240948167,
0.0468861767036331, 0.0102813794498838, 0.0693803856754811, 0.00566014993761021,
0.00275831473628789, 0.00234550137829788, 0.0116252548991317,
0.0143153463606759, 0.00995352784254985, 0.00238257560505637,
0.00182950280683248, 0.013195116994233, 0.000640921917352429,
0.000171523469251389, 0.0087256530793244, 0.0989386901919075,
0.0321103798387662, 0.0222773975090858, 0.007943310795799, 0.171053950985746,
0.012874269835152, 0.38817395138115, 0.00787508757030877, 0.0114487159712535,
0.0187831808209386, 0.0452479566115196, 0.00640720682677851,
0.0315482155790946, 0.125132833439637, 0.0115284490664364, 0.00538397509568388,
0.00517772137814985, 0.00176762574966497, 0.0273122011845722,
0.0191341545126795, 0.0235413416908084, 0.270326642321866, 0.179926862630332,
0.008734949388329, 0.0138918131322944, 0.00507817315892406, 0.000173786133243839,
0.00225184544628237, 0.0181059516295825, 0.541544494598043),
probename = c("Mark_1", "Mark_2", "Mark_3", "Mark_4", "Mark_5",
"Mark_6", "Mark_7", "Mark_8", "Mark_9", "Mark_10", "Mark_11",
"Mark_12", "Mark_13", "Mark_14", "Mark_15", "Mark_16", "Mark_17",
"Mark_18", "Mark_19", "Mark_20", "Mark_21", "Mark_22", "Mark_23",
"Mark_24", "Mark_25", "Mark_26", "Mark_27", "Mark_28", "Mark_29",
"Mark_30", "Mark_31", "Mark_32", "Mark_33", "Mark_34", "Mark_35",
"Mark_36", "Mark_37", "Mark_38", "Mark_39", "Mark_40", "Mark_41",
"Mark_42", "Mark_43", "Mark_44", "Mark_45", "Mark_46", "Mark_47",
"Mark_48", "Mark_49", "Mark_50", "Mark_51", "Mark_52", "Mark_53",
"Mark_54", "Mark_55", "Mark_56", "Mark_57", "Mark_58", "Mark_59",
"Mark_60", "Mark_61", "Mark_62", "Mark_63", "Mark_64", "Mark_65",
"Mark_66", "Mark_67", "Mark_68", "Mark_69", "Mark_70", "Mark_71",
"Mark_72", "Mark_73", "Mark_74", "Mark_75", "Mark_76", "Mark_77",
"Mark_78", "Mark_79", "Mark_80", "Mark_81", "Mark_82", "Mark_83",
"Mark_84", "Mark_85", "Mark_86", "Mark_87", "Mark_88", "Mark_89",
"Mark_90", "Mark_91", "Mark_92", "Mark_93", "Mark_94", "Mark_95",
"Mark_96", "Mark_97", "Mark_98", "Mark_99", "Mark_100", "Mark_101",
"Mark_102", "Mark_103", "Mark_104", "Mark_105", "Mark_106",
"Mark_107", "Mark_108", "Mark_109", "Mark_110", "Mark_111",
"Mark_112", "Mark_113", "Mark_114", "Mark_115", "Mark_116",
"Mark_117", "Mark_118", "Mark_119", "Mark_120", "Mark_121",
"Mark_122", "Mark_123", "Mark_124", "Mark_125", "Mark_126",
"Mark_127", "Mark_128", "Mark_129", "Mark_130", "Mark_131",
"Mark_132", "Mark_133", "Mark_134", "Mark_135", "Mark_136",
"Mark_137", "Mark_138", "Mark_139", "Mark_140", "Mark_141",
"Mark_142", "Mark_143", "Mark_144", "Mark_145", "Mark_146",
"Mark_147", "Mark_148", "Mark_149", "Mark_150", "Mark_151",
"Mark_152", "Mark_153", "Mark_154", "Mark_155", "Mark_156",
"Mark_157", "Mark_158", "Mark_159", "Mark_160", "Mark_161",
"Mark_162", "Mark_163", "Mark_164", "Mark_165", "Mark_166",
"Mark_167", "Mark_168", "Mark_169", "Mark_170", "Mark_171",
"Mark_172", "Mark_173", "Mark_174", "Mark_175", "Mark_176",
"Mark_177", "Mark_178", "Mark_179", "Mark_180", "Mark_181",
"Mark_182", "Mark_183", "Mark_184", "Mark_185", "Mark_186",
"Mark_187", "Mark_188", "Mark_189", "Mark_190", "Mark_191",
"Mark_192", "Mark_193", "Mark_194", "Mark_195", "Mark_196",
"Mark_197", "Mark_198", "Mark_199", "Mark_200", "Mark_201",
"Mark_202", "Mark_203", "Mark_204", "Mark_205", "Mark_206",
"Mark_207", "Mark_208", "Mark_209", "Mark_210", "Mark_211",
"Mark_212", "Mark_213", "Mark_214", "Mark_215", "Mark_216",
"Mark_217", "Mark_218", "Mark_219", "Mark_220", "Mark_221",
"Mark_222", "Mark_223", "Mark_224", "Mark_225", "Mark_226",
"Mark_227", "Mark_228", "Mark_229", "Mark_230", "Mark_231",
"Mark_232", "Mark_233", "Mark_234", "Mark_235", "Mark_236",
"Mark_237", "Mark_238", "Mark_239", "Mark_240", "Mark_241",
"Mark_242")), row.names = c("Mark_1", "Mark_2", "Mark_3",
"Mark_4", "Mark_5", "Mark_6", "Mark_7", "Mark_8", "Mark_9", "Mark_10",
"Mark_11", "Mark_12", "Mark_13", "Mark_14", "Mark_15", "Mark_16",
"Mark_17", "Mark_18", "Mark_19", "Mark_20", "Mark_21", "Mark_22",
"Mark_23", "Mark_24", "Mark_25", "Mark_26", "Mark_27", "Mark_28",
"Mark_29", "Mark_30", "Mark_31", "Mark_32", "Mark_33", "Mark_34",
"Mark_35", "Mark_36", "Mark_37", "Mark_38", "Mark_39", "Mark_40",
"Mark_41", "Mark_42", "Mark_43", "Mark_44", "Mark_45", "Mark_46",
"Mark_47", "Mark_48", "Mark_49", "Mark_50", "Mark_51", "Mark_52",
"Mark_53", "Mark_54", "Mark_55", "Mark_56", "Mark_57", "Mark_58",
"Mark_59", "Mark_60", "Mark_61", "Mark_62", "Mark_63", "Mark_64",
"Mark_65", "Mark_66", "Mark_67", "Mark_68", "Mark_69", "Mark_70",
"Mark_71", "Mark_72", "Mark_73", "Mark_74", "Mark_75", "Mark_76",
"Mark_77", "Mark_78", "Mark_79", "Mark_80", "Mark_81", "Mark_82",
"Mark_83", "Mark_84", "Mark_85", "Mark_86", "Mark_87", "Mark_88",
"Mark_89", "Mark_90", "Mark_91", "Mark_92", "Mark_93", "Mark_94",
"Mark_95", "Mark_96", "Mark_97", "Mark_98", "Mark_99", "Mark_100",
"Mark_101", "Mark_102", "Mark_103", "Mark_104", "Mark_105", "Mark_106",
"Mark_107", "Mark_108", "Mark_109", "Mark_110", "Mark_111", "Mark_112",
"Mark_113", "Mark_114", "Mark_115", "Mark_116", "Mark_117", "Mark_118",
"Mark_119", "Mark_120", "Mark_121", "Mark_122", "Mark_123", "Mark_124",
"Mark_125", "Mark_126", "Mark_127", "Mark_128", "Mark_129", "Mark_130",
"Mark_131", "Mark_132", "Mark_133", "Mark_134", "Mark_135", "Mark_136",
"Mark_137", "Mark_138", "Mark_139", "Mark_140", "Mark_141", "Mark_142",
"Mark_143", "Mark_144", "Mark_145", "Mark_146", "Mark_147", "Mark_148",
"Mark_149", "Mark_150", "Mark_151", "Mark_152", "Mark_153", "Mark_154",
"Mark_155", "Mark_156", "Mark_157", "Mark_158", "Mark_159", "Mark_160",
"Mark_161", "Mark_162", "Mark_163", "Mark_164", "Mark_165", "Mark_166",
"Mark_167", "Mark_168", "Mark_169", "Mark_170", "Mark_171", "Mark_172",
"Mark_173", "Mark_174", "Mark_175", "Mark_176", "Mark_177", "Mark_178",
"Mark_179", "Mark_180", "Mark_181", "Mark_182", "Mark_183", "Mark_184",
"Mark_185", "Mark_186", "Mark_187", "Mark_188", "Mark_189", "Mark_190",
"Mark_191", "Mark_192", "Mark_193", "Mark_194", "Mark_195", "Mark_196",
"Mark_197", "Mark_198", "Mark_199", "Mark_200", "Mark_201", "Mark_202",
"Mark_203", "Mark_204", "Mark_205", "Mark_206", "Mark_207", "Mark_208",
"Mark_209", "Mark_210", "Mark_211", "Mark_212", "Mark_213", "Mark_214",
"Mark_215", "Mark_216", "Mark_217", "Mark_218", "Mark_219", "Mark_220",
"Mark_221", "Mark_222", "Mark_223", "Mark_224", "Mark_225", "Mark_226",
"Mark_227", "Mark_228", "Mark_229", "Mark_230", "Mark_231", "Mark_232",
"Mark_233", "Mark_234", "Mark_235", "Mark_236", "Mark_237", "Mark_238",
"Mark_239", "Mark_240", "Mark_241", "Mark_242"), class = "data.frame")
I would like to create a nice graph (publication wise) to represent a data stored in this data frame. In general I would like to create a scatterplot or volcano plot with colors/shapes indicating what is important in my data.
I would like to achieve something like that:
Or:
As a filter cutoff we can start with: foldchange > 4 & all_pvalue < 0.05. I would like to also have a possibility to highlight (different color/shape and with a label) only couple of rows. Lets say as a starting point I would like to highlight Mark_23 and Mark_65.
Is it doable in R ? I have already tried something with volcano plot:
volcano = ggplot(data = data, aes(x = foldchange, y = -1*log10(all_pvalue)))
volcano + geom_point()
Can someone help me with going further ?
Here is an attempt to show what's possible
data$zones := interaction(abs(data$foldchange)>4,data$all_pvalue<0.05)
library(ggplot2)
library(ggrepel)
ggplot(data = data, aes(x = foldchange, y = -1*log10(all_pvalue),color = zones))+
geom_point()+
theme_light()+
geom_label_repel(aes(label=ifelse(probename %in% c("Mark_23","Mark_65"),as.character(probename),"")),
box.padding = 0.35,
point.padding = 0.5,
segment.color = 'grey50',show.legend = FALSE) +
geom_hline(yintercept = -1*log10(0.05),linetype = "dashed")+
geom_vline(xintercept = 4,linetype = "dashed")+
scale_color_manual(labels = c("bad", "also bad","good","not sure"), values = c("gray50","green","blue", "red"))+
labs(x = "fold change",
y = expression(log[10](p)),
color = "meaning")
there are many things here. geom_label_repel from library(ggrepel) allow to link your few points you want to show with a label. The easiest way of having the different colours is to create a variable that says in what zone you are. That is what I did with the zones variable, that you use in colour to have a changing color. You can manualy change it with the scale_color_manual function.
In the labs function I used the expression function that allows you to do superscripts. You can of course change the theme including, margins, axis, text size etc.
You can for example add theme(legend.position="top") to your plot to have the legend on top.