Skip to main content

plotting - Is it possible to define a new PlotTheme?


In version 10 one can define, say,


$PlotTheme = "Scientific"


Which changes the appearance of all plots. This is fantastic because on the one hand I can have nice plots very easily and globally (i.e. No need to tune every option in Plot[...]), and on the other hand when I exchange code with my collaborators they may not have Mathematica 10 and defining $PlotTheme is harmless for previous versions.


However, none of the provided plot themes fit my needs and I would need to combine some options (for example, the lines should on the one hand have colors, and on the other hand, has solid/dashed/dotted styles). Thus I wonder if it is possible to define a plot theme myself, combining a few built-in themes, and specify it globally using


$PlotTheme = "myStyle"

Note: This is related to Specifying non-conflicting PlotTheme options (which focuses on changing and tuning plot theme as PlotStyle, etc.), and @kguler has already provided a great answer.



Answer




There appears to be a mechanism for doing just that, though I have yet to map its capabilities.


As a basic example for the time being:


Themes`AddThemeRules["wizard",

  DefaultPlotStyle -> Thread@Directive[{Purple, Orange, Hue[0.6]}, Thick],
  LabelStyle -> 18,
  AxesStyle -> White,
  TicksStyle -> LightGray,
  Background -> Gray
]

Now:


Plot[{Sinc[x], Sinc[2 x], Sinc[3 x]}, {x, 0, 10}, PlotTheme -> "wizard"]


enter image description here


Hideous, I know. :o)


You can attach rules to specific plot functions using the second parameter, e.g. BarChart:


Themes`AddThemeRules["wizard", BarChart,
  ChartStyle -> {Pink, Gray, Brown}
];

Now:


BarChart[{{1, 2, 3}, {1, 3, 2}}, PlotTheme -> "wizard"]


enter image description here


These themes do not persist across a kernel restart so you can experiment freely I believe.
If you wish to make any changes persist you could use kernel/init.m.



I could not find a way to use AddThemeRules to make Themes that would combine with others in the way that the default ones will. I found that I needed to take things to a lower level and make assignments to this System function:


System`PlotThemeDump`resolvePlotTheme

This appears to be the true home of PlotThemes and one can look at its Definition to see everything, once it has been preloaded by Plot or some other means. To read a specific definition I (once again) recommend my step function. Let's check the definition of "ThickLines" for "Plot" (note that plot function names must be given as strings):


Themes`ThemeRules;  (* preload PlotThemes subsystem *)


Defer @@ step @ System`PlotThemeDump`resolvePlotTheme["ThickLines", "Plot"]


Themes`SetWeight[{"DefaultThickness" ->
 {AbsoluteThickness[3]}}, System`PlotThemeDump`$ComponentWeight]

(Defer replaces HoldForm to allow proper copy&paste.)


We can use this knowledge to create a new thickness Theme for Plot:


System`PlotThemeDump`resolvePlotTheme["Thick5", "Plot"] := 
  Themes`SetWeight[{"DefaultThickness" -> {AbsoluteThickness[5]}}, 

    System`PlotThemeDump`$ComponentWeight]

Now we can combine this with existing Themes just as we can the defaults:


Plot[{x^2 + x, x^2}, {x, -1, 1}, PlotTheme -> {"Detailed", "Thick5"}]

enter image description here


For more on the role of SetWeight please see:



Comments

Post a Comment

Popular posts from this blog

functions - Get leading series expansion term?

Given a function f[x] , I would like to have a function leadingSeries that returns just the leading term in the series around x=0 . For example: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x)] x and leadingSeries[(1/x + 2 + (1 - 1/x^3)/4)/(4 + x)] -(1/(16 x^3)) Is there such a function in Mathematica? Or maybe one can implement it efficiently? EDIT I finally went with the following implementation, based on Carl Woll 's answer: lds[ex_,x_]:=( (ex/.x->(x+O[x]^2))/.SeriesData[U_,Z_,L_List,Mi_,Ma_,De_]:>SeriesData[U,Z,{L[[1]]},Mi,Mi+1,De]//Quiet//Normal) The advantage is, that this one also properly works with functions whose leading term is a constant: lds[Exp[x],x] 1 Answer Update 1 Updated to eliminate SeriesData and to not return additional terms Perhaps you could use: leadingSeries[expr_, x_] := Normal[expr /. x->(x+O[x]^2) /. a_List :> Take[a, 1]] Then for your examples: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x), x] leadingSeries[Exp[x], x] leadingSeries[(1/x + 2 + (1 - 1/x...
Post a Comment
Read more

mathematical optimization - Minimizing using indices, error: Part::pkspec1: The expression cannot be used as a part specification

I want to use Minimize where the variables to minimize are indices pointing into an array. Here a MWE that hopefully shows what my problem is. vars = u@# & /@ Range[3]; cons = Flatten@ { Table[(u[j] != #) & /@ vars[[j + 1 ;; -1]], {j, 1, 3 - 1}], 1 vec1 = {1, 2, 3}; vec2 = {1, 2, 3}; Minimize[{Total@((vec1[[#]] - vec2[[u[#]]])^2 & /@ Range[1, 3]), cons}, vars, Integers] The error I get: Part::pkspec1: The expression u[1] cannot be used as a part specification. >> Answer Ok, it seems that one can get around Mathematica trying to evaluate vec2[[u[1]]] too early by using the function Indexed[vec2,u[1]] . The working MWE would then look like the following: vars = u@# & /@ Range[3]; cons = Flatten@{ Table[(u[j] != #) & /@ vars[[j + 1 ;; -1]], {j, 1, 3 - 1}], 1 vec1 = {1, 2, 3}; vec2 = {1, 2, 3}; NMinimize[ {Total@((vec1[[#]] - Indexed[vec2, u[#]])^2 & /@ R...
Post a Comment
Read more

How to remap graph properties?

Graph objects support both custom properties, which do not have special meanings, and standard properties, which may be used by some functions. When importing from formats such as GraphML, we usually get a result with custom properties. What is the simplest way to remap one property to another, e.g. to remap a custom property to a standard one so it can be used with various functions? Example: Let's get Zachary's karate club network with edge weights and vertex names from here: http://nexus.igraph.org/api/dataset_info?id=1&format=html g = Import[ "http://nexus.igraph.org/api/dataset?id=1&format=GraphML", {"ZIP", "karate.GraphML"}] I can remap "name" to VertexLabels and "weights" to EdgeWeight like this: sp[prop_][g_] := SetProperty[g, prop] g2 = g // sp[EdgeWeight -> (PropertyValue[{g, #}, "weight"] & /@ EdgeList[g])] // sp[VertexLabels -> (# -> PropertyValue[{g, #}, "name"]...
Post a Comment
Read more