我正在为模拟数据编写一个感知器学习算法.然而,程序进入无限循环,重量往往非常大.我该怎么做才能调试我的程序?如果你能指出出了什么问题,那也值得赞赏.
我在这里做的是首先随机生成一些数据点,并根据线性目标函数为它们分配标签.然后使用perceptron学习来学习这种线性函数.如果我使用100个样本,则下面是标记数据.
此外,这是关于从数据学习一书的练习1.4.
import numpy as np
a = 1
b = 1
def target(x):
if x[1]>a*x[0]+b:
return 1
else:
return -1
def gen_y(X_sim):
return np.array([target(x) for x in X_sim])
def pcp(X,y):
w = np.zeros(2)
Z = np.hstack((X,np.array([y]).T))
while ~all(z[2]*np.dot(w,z[:2])>0 for z in Z): # some training sample is missclassified
i = np.where(y*np.dot(w,x)<0 for x in X)[0][0] # update the weight based on misclassified sample
print(i)
w = w + y[i]*X[i]
return w
if __name__ == '__main__':
X = np.random.multivariate_normal([1,1],np.diag([1,1]),20) …Run Code Online (Sandbox Code Playgroud) 以下不起作用:
然而,这完全适用于Jupiter Notebook.
如果我只是将其评论出来,则图表不会显示出来.(也许它不会反正出现)
import pandas as pd
import matplotlib
from numpy.random import randn
import numpy as np
import matplotlib.pyplot as plt
df = pd.read_csv('data/playgolf.csv', delimiter='|' )
print(df.head())
hs = df.hist(['Temperature','Humidity'], bins=5)
print(hs)
Run Code Online (Sandbox Code Playgroud) 为什么这不会返回每个社区(边界框)中的点数计数?
import geopandas as gpd
def radius(points_neighbour, points_center, new_field_name, r):
"""
:param points_neighbour:
:param points_center:
:param new_field_name: new field_name attached to points_center
:param r: radius around points_center
:return:
"""
sindex = points_neighbour.sindex
pts_in_neighbour = []
for i, pt_center in points_center.iterrows():
nearest_index = list(sindex.intersection((pt_center.LATITUDE-r, pt_center.LONGITUDE-r, pt_center.LATITUDE+r, pt_center.LONGITUDE+r)))
pts_in_this_neighbour = points_neighbour[nearest_index]
pts_in_neighbour.append(len(pts_in_this_neighbour))
points_center[new_field_name] = gpd.GeoSeries(pts_in_neighbour)
Run Code Online (Sandbox Code Playgroud)
每个循环都给出相同的结果.
第二个问题,我怎样才能找到第k个最近的邻居?
有关问题本身的更多信息:
我们正在以非常小的规模进行这项工作,例如华盛顿州,美国或加拿大不列颠哥伦比亚省
我们希望尽可能多地使用geopandas,因为它类似于pandas并支持空间索引:RTree
例如,sindex这里有方法最近,交叉点等.
如果您需要更多信息,请发表评论.这是GeoPandasBase类中的代码
@property
def sindex(self):
if not self._sindex_generated:
self._generate_sindex()
return self._sindex
Run Code Online (Sandbox Code Playgroud)
我试过理查德的例子,但它没有用
def radius(points_neighbour, points_center, new_field_name, r):
"""
:param points_neighbour:
:param points_center: …Run Code Online (Sandbox Code Playgroud) 如何在 jekyll markdown 中设置图像大小?

Run Code Online (Sandbox Code Playgroud)
我正在使用 jekyll 最小值。使用

Run Code Online (Sandbox Code Playgroud)
从在降价更改图像大小不起作用。如果可能,我也想知道如何旋转图像。
为什么^不在Python中进行平方?我知道取幂是**相反的,但究竟是什么^以及为什么不使用该算子呢?
例如2^2=0,3^2=1.
python ×3
geopandas ×1
geospatial ×1
gis ×1
github-pages ×1
html ×1
jekyll ×1
markdown ×1
matplotlib ×1
r-tree ×1