我是python的初学者,我正在使用python实现主成分分析(PCA),但我在计算平均值时遇到了问题。下面是我的代码:
import Image
import os
from PIL import Image
from numpy import *
import numpy as np
#import images
dirname = "C:\\Users\\Karim\\Downloads\\att_faces\\New folder"
X = [np.asarray(Image.open(os.path.join(dirname, fn))) for fn in os
对于下面的代码,我无法得到A、B和C的单个值的标准错误。我尝试了几种不同的方法,要么得到0,要么得到一个单数值。下面的方法是我尝试过的最新方法,但没有奏效。有人有什么建议吗?我是Python新手。
import numpy as np
import pandas as pd
import scipy
import pylab
import matplotlib
import matplotlib.pyplot as plt
import math
import sympy as sp
from scipy.optimize import curve_fit
x = np.array([600,
我让python读取文本文件。
import numpy as np
with open("Input2010_5a.txt", "r") as file:
for line in file:
date, long, lat, depth, temp, sal = line.split("\t")
line_data = []
line_data.append(float(date))
line_data.append(float(long))
line_data.append(float(lat))
我试图在python中编写一个EM算法来对不同类型的图像进行聚类。我对EM算法的理解如下:
因此,我在python中编写了相同的代码。这是我的密码:
import numpy as np
import sys
from scipy import stats
# μ: mean vector ∈ R^(self.m x self.n)
# Σ: covariance matrix ∈ R^(self.m x self.n x self.n)
# π: probabilities of each component
class gmm:
def __init__(self,n
我尝试了拟合OLS模型,它在没有稳健估计的情况下正确工作,但是我想改进我的回归,所以,像下面一样,我尝试用这个问题来实现这个问题,在注释中有其他的尝试来解决它。
我不知道是否正确地应用关键字,所以我接受任何帮助。
代码:
# Fit and summarize OLS model
sumrz = dict()
for i, ca in enumerate(ccaa):
x = sm.add_constant(data.dy[ca])
mod = sm.OLS(endog=data.du[ca], exog=x, hasconst=True, missing='drop
我的任务是对给定的二维5000个数据进行PCA和白化变换。
我对PCA的理解是用协方差矩阵的特征向量分析数据的主轴,并将主轴旋转到x轴!
这就是我所做的。
[BtEvector,BtEvalue]=eig(MYCov);% Eigen value and vector using built-in function
我首先计算了特征值和向量。结果是
BtEvalue=[4.027487815706757,0;0,8.903923357227459]
和
BtEvector=[0.033937679569230,-0.999423951036524;-0.999423951036524,-0.
我正在尝试使用python通过主成分分析(PCA)来实现人脸识别。我遵循本教程中的步骤:
下面是我的代码:
import os
from PIL import Image
import numpy as np
import glob
import numpy.linalg as linalg
#Step1: put database images into a 2D array
filenames = glob.glob('C:\\Users\\Karim\\Downloads\\att_faces\\New folder/*.pgm')
filenames.sort()
我有一个这个方程:$Var(X+Y) = Var(X) + 2Cov(X+Y)+ Var(Y)$。我已经用python实现了这些,但结果并不相同。有人能解释一下我的代码出了什么问题吗? X = [random.randint(1,20) for i in range(10)]
Y = [random.randint(1,20) for i in range(10)]
X_PLUS_Y = [X[i]+Y[i] for i in range(len(X))]
print(statistics.variance(X_PLUS_Y)) # 108.23333333333333
print(stat
我试图使用GaussianMixture模型进行图像分割,所以我使用了两个组件,协方差矩阵type=“满”,并尝试运行使用Spyder3.6,这是附带的anaconda。下面是代码:
from scipy.misc import imread, imshow
from sklearn.mixture import GaussianMixture as GMM
import graph_tool.all as gt
from graph_tool.all import *
X=imread('2.jpg')
old=X.shape
X=X.reshape(-1,3)
gmm=G
我目前正在创建一个最小方差投资组合,并决定使用PortfolioAnalytics包的函数optimize.portfolio。不幸的是,在提取权重时,它们都是NA,即使我的返回值中没有任何NA值,这将是(从我的角度来看)导致结果权重为NA的唯一原因。我的数据集由多个资产(+5000)组成,每个资产有60个观察值(每月)。
library(PortfolioAnalytics)
library(ROI)
#index returns is an xts object consisting of 3800 stock Ids(columns) and 60 observat