恒星の特性を調べる
StarDataは,105 個を超える数の恒星の特性へのアクセスを提供する.3000個の恒星をランダムサンプルとして選ぶ.
In[1]:=
stars = StarData[{"RandomEntities", 3000}];
それぞれの実体について数多くの特性が提供されている.
In[2]:=
sun = \!\(\*
NamespaceBox["LinguisticAssistant",
DynamicModuleBox[{Typeset`query$$ = "Sun", Typeset`boxes$$ =
TemplateBox[{"\"Sun\"",
RowBox[{"Entity", "[",
RowBox[{"\"Star\"", ",", "\"Sun\""}], "]"}],
"\"Entity[\\\"Star\\\", \\\"Sun\\\"]\"", "\"star\""},
"Entity"],
Typeset`allassumptions$$ = {{
"type" -> "Clash", "word" -> "Sun",
"template" -> "Assuming \"${word}\" is ${desc1}. Use as \
${desc2} instead", "count" -> "4",
"Values" -> {{
"name" -> "Star", "desc" -> " referring to stars",
"input" -> "*C.Sun-_*Star-"}, {
"name" -> "CalendarEventName", "desc" -> "a weekday",
"input" -> "*C.Sun-_*CalendarEventName-"}, {
"name" -> "Word", "desc" -> "a word",
"input" -> "*C.Sun-_*Word-"}, {
"name" -> "Surname", "desc" -> "a surname",
"input" -> "*C.Sun-_*Surname-"}}}},
Typeset`assumptions$$ = {}, Typeset`open$$ = {1, 2},
Typeset`querystate$$ = {
"Online" -> True, "Allowed" -> True,
"mparse.jsp" -> 1.321761`6.572697926887039, "Messages" -> {}}},
DynamicBox[ToBoxes[
AlphaIntegration`LinguisticAssistantBoxes["", 4, Automatic,
Dynamic[Typeset`query$$],
Dynamic[Typeset`boxes$$],
Dynamic[Typeset`allassumptions$$],
Dynamic[Typeset`assumptions$$],
Dynamic[Typeset`open$$],
Dynamic[Typeset`querystate$$]], StandardForm],
ImageSizeCache->{114., {7., 15.}},
TrackedSymbols:>{
Typeset`query$$, Typeset`boxes$$, Typeset`allassumptions$$,
Typeset`assumptions$$, Typeset`open$$, Typeset`querystate$$}],
DynamicModuleValues:>{},
UndoTrackedVariables:>{Typeset`open$$}],
BaseStyle->{"Deploy"},
DeleteWithContents->True,
Editable->False,
SelectWithContents->True]\);
In[3]:=
StarData[sun, "Properties"] // Length
Out[3]=
ランダムに選ばれた恒星の有効温度を抽出し,そのうちでデータが利用できるものを選ぶ.
In[4]:=
temps = Flatten[
DeleteMissing[StarData[#, "EffectiveTemperature"]] & /@
Partition[stars, 100]];
温度の記述統計量をいくつか計算する.
In[5]:=
stats = {Min, Max, Mean, Median, StandardDeviation};
TableForm[{Map[#[temps] &, stats]}, TableHeadings -> {None, stats}]
Out[5]//TableForm=
恒星温度のヒストグラムをプロットする.
In[6]:=
h = Histogram[temps, 25, PDF, AxesLabel -> Automatic]
Out[6]=
ノンパラメトリック分布で温度の分布を推定する.
In[7]:=
skd = SmoothKernelDistribution[temps, "Oversmooth"]
Out[7]=
長い裾を持つパラメトリック分布で温度の分布を推定する.
In[8]:=
edist = EstimatedDistribution[temps, StableDistribution[a, b, c, d]]
Out[8]=
確率密度関数の推定値を比べ,太陽の温度に関連付ける.
完全なWolfram言語入力を表示する
Out[9]=